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15958 SURFACE HYDROLOGY: I-EXPLANATION OF PHENOMENA 
| KEY WORDS: Conservation; Erosion; Ground water; Hydrology; 
Infiltration; Mapping; Pollution; Rainfall; Runoff; Seepage; 
Topography; Unsaturated flow; Water resources = 
j ABSTRACT: Nonisotropic sloping soil leads to a horizontal flow component during i 
j | - rain and infiltration. Splashing of raindrops and flow in the transition zone between the — 
a 4 _ soil bulk and the air also lead to lateral flow. The result is concentration of moisture in 
- concave spots to the point of saturation and even seepage out of the soil. The same a 
| occurs at steep soil cuts. Such moisture concentration can explain the formation of — 
{ _ runoff during rainfall at low rates that do not exceed infiltration capacity, partial area 
contribution to runoff, net recharge of groundwater by low-rate rainfall, the formation 
| of rills and gullies in a relatively dry area, delay of runoff after a certain total amount 
| of rain irrespective of rain intensity, variations in soil formation in different parts of _ 
_ the landscape, poor distribution of irrigation water in nonplane fields, and leaching of — N 
| the soil surface by water that was assumed to run only above the surface. = 
¥ REFERENCE: Zaslavsky, Dan, and Sinai, Gideon, “Surface Hydrology: I— 
ig Explanation of Phenomena,” Journal of the . pean Division, ASCE, Vol. 107, No. 
_HY1, Proc. Paper 15958, 1981, 1, PP. 1-16 


_ 15959 SURFACE HY DROLOGY: I— DISTRIBUTION OF RAINDROPS ob 
_ KEY WORDS: Erosion; Ground water; Hydrology; Infiltration; Rainfall; cri 
Runoff; Sediments; Seepage; Topography; Water resources 
- ABSTRACT: -Nonuniformity of rainwater distribution « over a given soil area is due to 
two mechanisms. Raindrop splashes produce a net horizontal flow that is anf 
to the slope and to the rate of the rainfall. A change in slope causes excess 
concentration of rain splashes in concave spots and reduces the rain over convex 
points. The concavity can be measured geodetically. The size of the horizontal 
ts discharge was deduced theoretically and measured on a sloping plane in a rain = 
be simulator. Both indicated the same order of magnitude. The horizontal discharge over — 
short and steep slopes can reach a tenth of the total rain. The local excess rain can be 
1 up to half of the average rain. The other mechanism responsible for localized excess — 
a” rain is the slanting of rain over slope facing it, as distinct from a soil slope in the wm 
“shadow.” This mechanism can lead to local doubling of the rain Jo. 
a REFERENCE:  Zaslavsky, Dan, and Sinai, Gideon, “Surface Hydrology: II— 
a Distribution of Raindrops,” Journal of the Hydraulics Division, ASCE, Vol. a No. 
HY/1, Proc. Paper 15959, January, 1981 
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15960 190 HY pRowocr: CAUSES OF LATERAL 

KEY WORDS: Erosion; Ground water; Hydrology; Infiltration; 


Rainfall; Runoff; Seepage; Topography =| “ 
| Penetrating to a more permeable soil layer they bend downstream. This can occur in 
t unsaturated soil above a sloping water table or in a layer in which water accumulates 4 
3 - because of the presence of less permeable underlying layer. We considered the general s 4 
case of a nonuniform soil with rotational symmetry at an angle to the vertical, i.e., a . 
“ sloping soil. General proof is provided that the soil generally behaves as a anisotropic — q 
7 medium so that infiltration also invloves a horizontal net flow component. This net 
flow component is proportional to the slope, the vertical flow, and the coefficient of | 
- anisotropy. The existence of such a mechanism can account for the concentration of 
water in concave spots, possible seepage, and runoff and erosion phenomena at low _ 
rates of rainfall. It can also explain how a previous rainfall leads to earlier onset of 
purenence: Zaslavsky, Dan, and Sinai, Gideon, “Surface Hydrology: I1I—Causes of 
— Flow,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY1, Proc. — 
Paper 15960, January, 1981, pp. 37-52 
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: and the slope; (2)moisture concentration increases with the concavity of the soil 


KEY ‘WORDS: Ground water; Hydrology; Permeability; 

Pollution; Rainfall; Runoff; Seepage; Water pollution 
7 

ABSTRACT: In a nonuniform soil with rotational symmetry the 
along a coordinate n that forms an angle a with the vertical. During infiltration there 
_ will be a flow component in the parallel direction s downstream, which is proportional — 

to the rain itself and the slope. The horizontal flow is proportional to the vertical flow, 
_ the slope, and the coefficient of anisotropy. The last has been calculated for a fl : 
_ made up distinct layers. In the limiting case of very thin layers or near the saturation __ 
point, the vertical and horizontal conductivities are, respectively, the expected — 
arithmetic and harmonic means familiar from saturated flow. The coefficient of 

_ anisotropy increases with the rate of vertical flow to a power greater than unity. The — 
_ horizontal flow component can be many times the vertical one. It finds application in a3 
pollution and in the study of the water regime. 

REFERENCE: Zaslavsky, Dan, and Sinai, Gideon, ‘Surface 
- Sloping, Layered Soil,” Journal of the Hydraulics Division, ASCE, Vol. 107, No. mYi, 


Proc. Paper 15961, January, 1981, pp. 53-64 


15962 SURFACE HYDROLOGY: V—INS SURFACE ACE TRANSIENT FLOW ur wee? 


‘KEY WORDS: Infiltration; Irrigation; ; Permeability; 
_ ABSTRACT: Seantadacd computation of nonsteady flow in the transition layer at the 
_ soil surface was carried out with various geometric dimensions and rain regimes. The 
- main conclusions are: (1)Lateral flow in the transition layer is proportional to the ~f 


_ (3)saturation initially occurs at the most concave point within the transition layer; and 
_ (4)the saturation zone may spread to the soil surface and cause seepage runoff and 
- erosion. Both saturation and seepage occur after a certain amount of rain has fallen | 
and are little dependent on its intensity. At very low rates of rainfall, no seepage will 
' occur if the transition layer is sufficiently thin. Gullies and rills are not merely the 

result of seepage at concave spot—they themselves produce more and earlier seepage _ 

_ runoff and erosion; seepage of rain and runoff from concave spots will occur when the _ 

amount of rain is either greater or less than the infiltration capacity, 

_ REFERENCE: | Zaslavsky, Dan, and Sinai, Gideon, “Surface Hydrology: V—In- Surface | 

Flow,” Journal of the Hydraulics D Division, Vol. ol. 107, ™ HY1, 

Paper 15962, January, 1981, pp. 65-93 } 


_ surface and continues long after the rain stops occurring mainly in unsaturated soil; i 
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15986 MULTIPORT SLEEVE VALVE DEVELOPMENT AND — 
WORDS: Cavitation; Design; Energy dissipation; Flow regulators; 


ie Fluid flow; Hydraulics; Jets; Laboratory tests; Pipelines; Sleeve valves; 


- ABSTRACT: The history of the development and use of multiport sleeve valves is 
| reviewed. The concept of using the multiport sleeve valve as an energy dissipator as 


_{ well as a control valve for water pipelines is explored. It is expected that high-pressure — 


_ head flows (150 to 300 m) can be controlled with the use of the multiport sleeve valve. 

_ Head loss coefficients for two sleeve valve designs are presented based on laboratory 

tests of 200-mm valve. Physical and mathematical models are described which were 

used to develop a horizontal multiport sleeve valve for long pipelines. The multiport 

_ Sleeve valve is quite versatile with respect to hydraulic and mechanical design. rd “1 

| REFERENCE. Burgi, Philip H., Green, Edward O., and Thibault, Raoul E., 
“Multiport Sleeve Vaive Development and Application,” Journal of the Hydraulics 

_ Division, ASCE, Vol. 107, No. HY1, Proc. Paper 15986, January, 1981, pp. 95-111 +3 
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SURF. SURFACE E HyDROLOG :I— —EXPLANATION | 


Zaslavsky ' and Gideon Sinai” 


¢ The following is the first of a series of articles devoted to a new — 
‘in the study of surface hydrology. The main premise is that during rain and — 
flow in a nonsaturated soil, a horizontal flow component may exist which leads 
_ to a number of phenomena as yet unexplained or, at best, poorly explained. 2a 
‘The s senior writer was asked, as long ago as 1964, to review the engineering 
_ practice concerning erosion in water channels as well as its rationale. The resulting _ 
_ Study indicated that there is very often little correlation between high-water 


of water could be found. There had to another process by 
rain percolates into the soil and then rises again to form seepage forces. "2 ee 
In 1968, the senior writer was asked by J. Van Schilfgaarde of the Agricultural — 
Research Service in Beltsville to join the Hydrograph Laboratory for the purpose 
of reexamining the approach to surface hydrology in in ones and, in particular, 2 


we 
Serious doubts have been raised as. to o the of the elementary physics 


_ phenomena seem clearly to contradict existing theories, while the very definition 
7 of the soil surface (42) poses a fundamental is somewhat philosophical question. 
The main points in the new development can be summarized as follows. The 


— soil is ; layered, its surface being the interface between the bulk of the soil 


- _ and mathematics involved in the classical approach. A number of observed 4 


 ' Prof. of Soil and Water Engrg., Faculty of Agricultural Engrg., Soils and Fertilizers 
‘Diy., Technion—Israel Inst. of Tech., Technion City, Haifa, 32000, Israel, 
*Lect., Faculty of Agricultural Engrg., Soils and Fertilizers Div., " Technion—Israel — 
Inst. of Tech., Technion City, Haifa, 32000, Israel. 
_ Note.—Discussion open until June 1, 1981. _ Separate discussions should be submitted 
- for the individual papers in this symposium. To extend the closing date one month, 
 awritten request must be filed with the Manager of Technical and Professional Publications, -, { 
a ASCE. Manuscript was submitted for review for possible publication on March 23, 1979. 
| _ This paper is part of the Journal of the Hydraulics Division, Proceedings of the American 
Society of Civil Engineers, © ASCE, Vol. 107, No. _ HY1, January, 1981. ISSN 0044- el 


¥ 
q | 
i _ runoffs or surface shear forces and the formation or widening of rills and gullies 
The importance of seepage forces and piping mechanisms has been realize@ 
| 
| 
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. and the air. These nonuniformities behave on the whole as a nonisotropy. In bia 
sloping soil this is the cause of the horizontal flow component which, in turn, a a 
may cause local moisture concentration, runoff at low Tates of rain, and and seepage 

This article, the first in a series, presents some of the basic phenomena that 

- surface hydrology must explain. It establishes the basic concept of lateral flow, 
at fires i in the form of a hypotheses. A number of conclusions are drawn which 
not only purport to reformulate the science of surface hydrology, but have © 

= direct practical applications in irrigation and soil and water conservation. le 

_ The subsequent articles will elaborate the theory and consider the Partial 

a experimental validation of the different parts of the three suggested mechanisms. - a 
These mechanisms are: (1) The splashing of raindrops; (2) the flow in the 

a at the soil surface; and (3) flow in : soll. ; 


g neni in either of two ways: (1) The rate of rain exceeds the absorptive capacity 
of the soil at a given point, implying that at a lower rate of rain there will © 
£ be no runoff; and (2) a groundwater table or a perched water table is built 
up which eventually produces. seepage or outflows from the soil. Such an outflow © 
i... groundwater has always been considered strictly as a part of streamflows © 
or baseflows with a large time-lag, and certainly not occurring during one 
rainstorm. In fact, the time- lag towards the end of a runoff event, or between 
- such events, was the only accepted way to identify such a baseflow. :) =" ee 
These concepts are the basis of papers from the early 1° 1920s | onwards (such 
as Refs. 31, 32, 29, 36, 28, 22, 23, 24, 4, and 33) as well as many others — 
_ that followed. A typical statement, found in a comprehensive textbook (33) 
reads: ‘Storm runoff ordinarily is considered to be that portion of the total 
runoff . . . which reaches the point of measurements within a relatively short 
period of time . . . The time involved varies from a few minutes . sehind 
remaining flow of a stream consists of ground water discharge which lags = 
' the preceding rainfall by an appreciably longer period of time than storm runoff, 
the period ranging from a matter of days to a number of years.’’ Or, “The - 
a maximum rate at which a small plot can absorb rainfall is called its ‘infiltration Am 
capacity . . . most of the difference betweeen storm rainfall and storm runoff 
consists of water that has passed through the surface of the soil into the ground — 
No statistical organization of these models, such as in Refs. 34, 13, 
10, departs materially from the basic concept that runoff is constituted principally — 
of rainfall less infiltration. Stochastic models (such as 9, 11, 12, 21, 30, Zz ' 
8, 39), or deterministic models (27,25), | can be cited in this connection. Two. 
notions underlying all these treatments, whether they are stated | explicitly 


or not, , concerning the formation ob ot yo! 


1. The rate. of rai rain exceeds the infiltration capacity of the soil at a point, 


_ implying that at a lower rate of rain there would be no runoff. 
wal 2. There is a buildup of a table ¢ or of a water 


il 
wy 


one rainstorm. mulziom # diiw no Ust noesse yaim 


_ A physically and mathematically more rigorous treatment of saturated and 


unsaturated flow in the soil is represented by works such as Refs. 14-20, a 
3, and 1, but they too do not, in effect, deviate from the two concepts that 
_ water will run off either because of its inability to penetrate the soil or after 
first accumulating in the groundwater. Wherever the groundwater seeps out 
of the soil, runoff may be formed, but this is dealyed until long after the rain 
Actual observations of runoff occurring during a rainstorm within a delay 
as short as a few minutes, or even only a fraction of a minute, indicate that 
the foregoing concept cannot be general and at best draws but a partial picture. 
While it appears to be evident that rain which exceeds the infiltration capacity — 
of the ground will run off (on condition that there is no surface detention), 


print to streamflows and baseflows, s, after. a long delay, but not 


light rains much lower than the infiltration capacity. A number of hypotheses 2 ; 
_ have been advanced to account for this fact, yet they often hardly as much 


is deduced from runoff measurements rather than based on direct observation. 
While there is no doubt that less permeable areas may contribute more surface 
rainfall excess, it is also certain that ‘tunel, may @ occur without such areas * 
and at low-intensity rains. Vino" iW Of Ol igus! ihw aint wen 
_ The infiltration capacity of the soil has been presumed to be a unique figure. 
This simplistic concept, which has prevailed for so many years, does not allow 
for the dependence of runoff on antecedent moisture. In more sophisticated 
treatises (e.g., Refs. 14-20), a more realistic picture of unsaturated flow in 
a the soil is admitted. We shall not refer here to a number of articles that relate - 
the rate | of infiltration to the water | storage in the soil. The actual nature of 
the phenomenon can best be understood by consulting Braester (6). According — 
to this work, surface moisture gradually increases during the rain, the nor 
capacity of a uniform soil is simply its hydraulic conductivity, and the —_ 
& approaches saturation, after a long time, if the ri rate ¢ of rainfall | equals or surpasses — 


a certain ‘amount of | time is needed for surface s saturation on and flooding | to 
_ commence. This time will depend on the rain intensity, on a number of soil 
: properties in addition to the saturated hydraulic conductivity, and on antecedent — 
moisture. All papers relating the infiltration mechanism to rainfall and Lanesndll ; 
assume one dimensional vertical flow. In a nonuniform soil, a different definition — 
of the infiltration capacity isrequired. 
A high antecedent moisture alone cannot account for runoff phenomena. 

_ Analysis and observations of unsaturated flow have indicated that a very short 
’ time after the cessation of a rainshower the vertical soil column reaches a ad 
_ or less fixed moisture content known as the field capacity. Intervals between _ 

_ rainstorms as long 2 as a few days are insufficient to evaporate more than a 
few millimeters of water from the ‘soil, an nd often not even tha 


g | 
it still does not explain why sometimes quick runoffs are caused by relatively - 
| 
It presupposes that there are some parts of the soil surface which have a very 
low infiltration capacity, but the existence of such insufficiently permeable areas 
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evaporated can abe made | good bya few minutes ori rain. ean many y places; therefore, 
_ rainstorms during a rainy season fall on soil with a moisture content that varies 
within very narrow limits. Still, it is common experience that there i isa cumulative ee 
effect with gradually increasing runoffs. It is equally common experience in 
many regions that little or no runoff occurs before a few hundred millimeters — 
of rain have fallen. In the hills of Judea and Samaria, Israel, the relevant threshold _ 
is around 400 mm of rain, which is hardly affected by the time distribution. 4 
The intensity of a given storm has only a subsidiary effect on the runoff, a” 
constituting an addition to the accumulated total rain and antecedent moisture. 
The rise of groundwater, at any rate in Israel, cannot serve as an explanation. — 
7 These well known observed facts which ran counter to a preconceived physical - 
explanation, have merely spurred investigators to invent statistical tricks and a 
_ obfuscating factors. The worst part of these exercises has been that these Statistical a4 
abstractions were given names that pretend to grace true physical entities. The — 
| - dilemmas remain: How does runoff form when the rainfall does not exceed — 
; the infiltration capacity over the whole field? Why is there an increase in runoff 
At least two more ideas ‘should be mentioned ‘that attempt to explain ‘runoff. Hy 
- One is the formation of a surface crust (35) which has been shown to develop — 
s direct correlation with the cumulative rainfall. In experiments run by rain 
- simulators it was found that the infiltration capacity is reduced with increasing 
rainfall, e.g., some wind-blown loess soils of Israel can strat out with an infiltration | 


rain totalling some 200 mm, although there will be some recovery of surface : 
_ permeability after each drying period. The initial rate of infiltration for any 


_ new rain will be at least 10 mm/h-20 mm/h. Only the final rate of infiltration, - 
which is obtained after a few tens of millimeters, will be very low. bi minal 


indicate, however, that runoff starts much earlier. Other soils are much less 
sensitive to puddling by the rain. The crust formation can explain only part 
Be. other concept, which has already been mentioned, is that of a partial 
_ contribution or partial area. It states that small parts of the soil surface avs 
avery low infiltration capacity and thus contribute considerably to runoff while 
_ the other parts do not at all. This concept cannot either be proved or disproved. 
It is only another way of saying that there must be some reason for runoff 
_ despite the fact that the rainfall does not seem to exceed the infiltration capacity. = 
_ There can — be parts of the area where the rain exceeds the infiltratior 
In the following it will be demonstrated how some pa parts of the landscape 
contribute more to runoff than others. These parts, however, are not ome 
_ less permeable. Rather they are related mainly to the topographic configuration. - 
It would be erroneous to consider infiltration as one dimensional vertical flow. — 


4 This analysis follows largely the ideas set out in Refs. 37, 42, 43, and Sho 
_ Scalars and Vectors.—Traditionally, infiltration into the soil has been almost 
od synonymous with vertical flow. In reality, vertical flow is only one of three 
(but certainly two) flow components. horizontal-flow component cannot 


capacity of 30 mm /h-40 mm/h and end up with 3 mm /h-5 mm/h after cumulative z 


be added arithmetically to the infi iltration as if both were scalars. The commonly 


equation bom bowerem od pie q 


in which R = runoff; P = ant I |= infiltration can serve at 
as an overall scalar balance over an area in which the infiltration, J, is not 
actually measured at a point but is the difference between the measured 

precipitation and the outflow through a well-defined and measurable river or 
channel (assuming the boundaries of the drainage basin to be determinable by 
7 its topography alone). Thus Eq. | cannot be considered an equation for a specific | 
point, nor is it capable of predicting runoff, but it is to calculate net recharge — 
over a field, employing the measured values of the rain and the runoff. 
_ Errors in P and /.—Precipitation can be measured, albeit with limited accuracy 
(on the order of +20%). _ Infiltration capacity can be measured or i 


| iter in time or in space. A change by : as much as half a an an order | of a 
It is therefore unrealistic to expect much by w way of accuracy in the prediction 7 


of runoff, R, which is most commoly 5%-10% of the P (precipitation). Eq. 
AL or r any similar equation of differences, sophisticated as as it may look, cannot 


~ Rational Method. —The rational method expresses the runoff as a proportion 
of the rainfall. Representing a large volume of experience, it has the operational an 


_ recorded or estimated rain. So far it has lacked a convincing physical explanation - 
and is based almost exclusively on intuition. = 
- §.C.S. Method.—As described in the National Engineering handbook of the 
United States Department of Agriculture, the S.C.S. Method improves on ‘the 
Rational Method by allowing for an initial detention of rain and a a 
relationship between runoff and rain. It is interesting to note that on a large os 
proportion of very common soils, runoff does not begin before several inches" 
of rain have penetrated the soil. This is represented by the S.C.S. Method 
as being a function of the | total rainfall irrespective of the rainfall distribution 
7 Is it Possible to Measure Surface Runoff?—The answer to the question of 
; measuring runoff, R, depends on its definition. If by runoff is meant the outflow - 
-_ through a well-defined channel that drains a sizeable area, then such outflow 
is reasonably meaningful and measurable. However, at any given point in the 
field or when defined as overland sheet flow o or, as it is sometimes more ‘etme 


_ called, runoff supply, it defies unequivocal measurement as well as definition 
it is as difficult to define surface runoff as it is to define the soil surface 
itself. This problem of definition will be treated later in this report. It is reasonably » ll 
clear that at the soil surface there is a transition between the soil bulk and — 
the air. In the direction from the soil bulk outward, porosity as we well as hydraulic 
_ conductivity gradually increase. It may be hard to accept this very fundamental * 
argument about the transitive nature of the soil surface, but at least the practical 


result will strongly depend on the depth at which such interception will be 
common are is to a “deep enough” cutoff and let 


_ problem of intercepting the runoff for measurement will be appreciated. The = 


| 

| 


affects the quantity to be measured and undoubtedly tends to exaggerate the 


_ value of the apparent runoff. The alternative method is to provide a very thin 
horizontal threshold that is supposed to divide between the runoff and the flow 

within the soil. The question here is how thin the threshold is and at what 
elevation is it installed. The soil surface being irregular, it is legitimate to ask 


what is the size of the irregularities. to which the dividers should | conform in 


_ It would appear to be much more reasonable, and in fact feasible, to measure 


order to separate the soil bulk from the air. y 


the horizontal flux component or, still easier, the horizontal discharge (by the a 
vertical integration of the horizontal fluxes). This, in fact, is what is measured — 
near j a a vertical cutoff. Stagnation n near ‘such a a wall may ¢ cause > part of the water 


the two. The problem is not one of a mere technical limitation but is of a 


: place within the soil ¢ or + outside of it. _ There i is no means to distinguish tetvens 


q In summary, the notion of on-surface runoff as point values distributed over 


the soil area is fundamentally questionable and practically impossible to measure. i 
- Popular « as it is (from kindergarten onwards), the model of Eq. | still lacks " 

a real demonstration of relevancy to either a physical understanding, a consistent a 
mathematical formulation, or a practical measurement. 
Concentration oF Vi Water In Concave a 

of moisture in concave areas has been observed (38). This 

¥ oncentration occurs without any runoff in the usual sense of the word. By % 
‘concave”’ it is meant not only the bottom of a slope or of a valley but — Be 


~ transition from a steep to a moderate downward slope. Concentration of moisture _ 


without runoff in concave areas is a phenomenon that cannot be explained Si 
by any existing hydrological models. . Such a concentration has been observed 


J by the writers in a region of sand dunes in the Northern Sinai with 70 mm 7 
_rain/yr and an infiltration capacity of some 500 mm/h. It occurred both on 
_ steep and on moderate slopes. When travelling in the southern part of Israel 
_ following the rain, the green sheets of seasonal grass and shrubs can be seen J 
_ Covering concave parts of the landscape. Bedouin are wont to plant their barley ; 
only in concave . parts of the landscape. The accumulation of moisture has been = 
observed in areas where no surface runoff could possibly have taken place, 
where no water table was present, and where no highly impermeable layer 
In the sandy flowerbulb-growing area near Lisse (The Netherlands), it me 
been « observed that in a soil cross section under a concave surface the ‘sand 
«| 
was wetted to a considerably greater depth than under level surfaces (according 
to a personal communication from van der Valk and Knotterrus). This in turn = 
had its effect on moisture availability to plants and on wind erosion patterns. _ > 


_ Looking at fields under rain or after a rain, water often seems to appear 


» od Concave parts in fields often suffer from excessive wetness, bad trafficability, 


| 
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for a long time even wi there isnowatertable. q 
Any model attempting to explain surface hydrology should be compatible 
with this phenomenon of moisture accumulation in concave areas. It must ~~ . 
noted that surface runoff will not stop at a concavity caused by a transition _ 
from a slope of 12% to a slope of 10%. Every drop of water will either = 
vertically or become runoff with no accumulation at the point of concavity. — 
A new mechanism must be postulated to explain this phenomenon. —t*™*S 


‘Main Osservarion in Beer-Snesa Exreriment (37,38) 


2, at different points according to the formula 


7 The e elevation was measured at grid points | i i and j jas defined i in Fig. 1. i 


ij-t 


4 
of 
—Twe Grid for Measured Soil Surface Elevation 
- Moisture contents were measured at depths of 20 cm and 40 cm over an 
13 area of 70 m x 70 m. The field was planted substantially parallel with the © 
: ‘contour lines. The main slope v was 12%. Yields were also measured. The results | 
: a summarized in Fig. 2, the concavity being estimated by Vv? z as in Eq. 
2. The correlation between the moisture concents, C, two weeks after the rain, 
and t the curvature, was found atr= 90% tobe 
For the ‘one of brevity, detailed yields (which were correlated 
<a with the moisture) will not be cited here. Overall yields reached more 
than 2 tons/ha in the concave part (concavity equals 0.08 — 0.1 m~ ‘ in Fig. — 
nd and as little as 0. .2-0.3 — in some convex parts [concavity (- —0.09) 


runoff i in » the usual sense was observed. The soil was a loess-loam which - - 
A 


| 
7 
The curvature of the soil surface was measured geodetically through the elevation, _ 
4 
SS 
| 


was s plowed disked in the ordinary | water table, 
win water, or impermeable layer was found anywhere near the area dealt | 
A content measurements were made (salinity, fertility, clay content, etc.), but 
no visible trends that could be related to the variations in the moisture content — 
observation made in this experiment quantifies the ‘qualitative i impression 
- gained from observing many areas. It must be assumed that there were little 4 
_ or no differences in the rainfall distribution. The soil at the surface was uniform, _ 
fo that there were only small differences in the infiltration rate of the rain. 
_ The rain wetted the soil to a limited depth only, so that there was ; no net 


| of the groundwater. _ Where, then, did the differences in moisture 


of lateral flow within the soil. This iow would have to be such as to occur ' 


, FIG. 2.—Correlation of Moisture Content at 20 cm and 40 cm Depth 10 Days aft a 
] 


Rain with Soil Surface Curvature or Concavity 


7 in a nonsaturated soil and to lead towards accumulation in a concave sai 7 


- The concentration of water in concave pote of the landscape can iin 
oi some soil forming processes. The pedological ‘‘genetic’’ formation of a clay 
’ horizon is already very pronounced on flat land, but the layer so formed = 

4 is found to be thickest in concave parts (excluding hydromorphic alluvial bottom — 


land). On the high-lying convex parts of the slope, the development of B horizon 


_ Many theories have related this distinction to overland flow and erosion, 
but they a are, in fact, ‘unable to steep concave slopes 


sas") 

x 
q 
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can erosion carry away any B horizon from, under an A 
_ what way can it leave a distinct if faintly developed B horizon? Can it really ; > 
be that some tens or hundreds of years of B horizon development are followed 
In the opinion of the present writers, the development of a loamy or clayey 
B horizon on dune sand parent material can serve as a perfect model capable 
_ of unsettling existing theories and of offering some new insights. Observations : ; 
in Israel have clearly shown that the clay is imported into the sand by rainwater _ 
and settling dust. It is equally evident that water is the vehicle by which the 
clay enters the sand to accumulate in the B horizon. Dune sand has a hydraulic — 
conductivity of some 10° cm/day. Only an unusual rain-squall lasting a few 
minutes may have an intensity exceeding this hydraulic conductivity, the 
_ infiltration capacity of such short-term rain spurts actually being several times 
larger (5,6). It is very unlikely that any runoff and surface flow will occur — 
on sand dunes. How, then, can clay, carried by water, accumulate in concave | 
_ parts of the landscape without surface runoff, a high water table, or impermeable = ; 
_ The concept of lateral flow in the unsaturated soil has been introduced by 
, | a the writer (42). According to this concept, lateral flow is not caused by the | 
boundary conditions, such as the position of the water table, but by soil anisotropy 
o caused by soil layering. When the layers are horizontal, the main driving force rs 
a (gravity) is normal to the layers, so the flow is straight down. When the layers 
are at an angle to the horizon, the force of gravity points downstream from a 
‘the normal and will therefore cause a horizontal flux component pointing downhill. 
_It was thus reasonable to postulate (and was, in fact, later proven) that, 
g least under steady state, the average horizontal flux J, is proportional to — 
the average vertical flux, 7 , to the slope, tan a, and toa coefficient of anisotropy, 4 


"Tie single expression for lateral flow can have a number of possible 
explanations. In a concave landscape, the incoming slope is larger than the © 
. _ Outgoing one. Therefore, the incoming horizontal flux is higher than the outgoing = 
one, leading to moisture accumulation in low-lying (“‘concave”’) parts. In the _ 
_ mathematical terms for two dimensional 1al problems ems (z v vertical and x t horizontal) = 


Then, assuming as a first approximation that q, and U are not functions 

of x, the conservation law applied toEq.4yields 


a 

‘Manage to “‘catch’’ more water or clay than others to form a thick B horizon. 
U tan At q 
q « 
| 
=4, U 


in which = the me moisture content; and = the additional outgoing water 
3 caused by increased vertical flux. In short, any concave part of the landscape, 
_ (@°z/ax”) > 0, leads to higher moisture contents and a higher share of vertical 
infiltration. The more water there i is, the greater the development of the B 


horizon. Furthermore, on convex parts dons the landscape there is a shortage 4 


: moisture and a smaller vertical flow (a7 z/ax* < 0). In fact, as the B horizon 
develops, the anisotropy U does likewise, and g* grows progressively more 

“negative, The development of the B horizon in convex parts of the landscape 
then ceases of its own accord. This new explanation is interesting inasmuch 

as. it also interprets the fainter B horizon on straight and convex slopes - i 

_ genetically mature soil catenas. This explanation does not have to rely on the 
questionable crutch of erosion and runoff theories. 
_ The lateral flow component and the moisture accumulation on concave parts 

_ of the slope could explain the observations made in the Beer-Sheba experiment. 

If the lateral flow component and the moisture accumulation can be shown 

- to be of sufficient magnitude, they could explain the occurrence of saturation 

in parts of the field, leading to seepage out of the soil and the formation of 
overland flow. The ‘“‘partial area contribution’’ would then get a new meaning. 

Rain falling on areas with surface seepage will not infiltrate into the ground, p 
- but not because of the limited infiltration capacity, the value of which is totally 
irrelevant in this case. In addition, part of the area is waterlogged not because — 
of the ratio of rain to the infiltration | capacity. Tt is, ;, rather, _ because of some a 
mechanism leading to a local accumulation of water. The new concept, put 
forward by the present writers, postulates that this accumulation of water leads, _ 
in the first place, to locally increased groundwater recharge and may also lead 
to runoff. This is contrary to the traditional measuring of ‘‘partial area contribu- | 
tion,”’ * according to which there is less ‘groundwater recharge where runoff a. 

If the rain is able to penetrate into the ground, if only to a 0 a shallow. depth, 

: before turning into overland flow, it has an altogether different effect on leaching 
and pollution. In an appreciable number of observations of stormwater composi- — 
tion, the relatively high concentration of some solutes could not be explained 7 
unless the rainhadenteredthe soilfirst. || | 
_ If localized waterlogged areas can be formed by lateral flow, this would explain .? 
why a certain amount of total rain is needed and some time must elapse before 
significant runoff occurs. Analyzing Eq. 6 indeed shows that the rate of excess — 


accumulation, (aC)/(at) + q* » depends on the rate of rain, 


tz 
my aft 
mount rain. Waterlogged areas” can no longer absorb ‘Tain regardless 


od the rain first penetrates ‘into the ground and then seeps back out again 


‘in concave parts, this would explain field erosion by the formation of seepage — 
forces. Depressions in the soil surface are often a starting point of erosion, 
_ which can now be accounted for by local high accumulations of moisture and ; 


the concentration of streamlines, which in turn produce high seepage forces 
- Road cuts, steep river banks, and landslides truncate the soil layers. The 
7 tone flow seamen ae cu soil surface but cannot seep out as long as the 


i 
| 
i 
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soil i is ‘ius or the water pressure is negative. The streamlines then | bend 


downward and accumulate until sufficient pressure builds up for seepage — 
occur usually followed bad 
net groundwater ia in some areas of very limited rain. This is a no “a 
form of ‘‘water harvesting’’ whereby certain parts of the landscape are allotted a 
_ Several times more rainwater than the average. This runs counter to the a : 
4 ‘partial-area-contribution”’ theories, which state that the parts that contribute ‘ 
M to runoff are excluding most of “the rainwater by being impermeable. The lat — 
flow concept, on the other hand, clearly s shows that waterlogged areas can 
_ form, exhibiting a greater measure of deep percolation and supplying runoff 
Ithas been suggested (42) that all soils exception have amore permeable _ 
layer at their surface. This is due to the simple fact that the porosity approaches 
_ unity towards the air. This in itself produces a lateral flow component wherever — _ 
‘the soil surface : slopes. It has also been suggested (37,43,44) that while s splashing _ / be 
_ faindrops produce a real lateral flow component, very much as in Eq. 4, they © 
do not cause water to seep out of the soil. A detailed treatment of the mechanisms — 
_ leading to lateral flow will be the subject of future parts of this report. 
| Eq. 4 is proven to be physically sound, then it will have what may be > 
termed a fringe benefit in the accuracy of predicting runoff. Given errors in 
_measuring the rain, 7, will produce only the same relative errors in the horizontal 


_ flow component similar to the rational runoff equation. There is no amplification — ; 


fae the relative error because of the smaller value of the runoff compared with 7 


the rain and the infiltration as in Eq. 1. In fact, e- 4 does not : assume the — 


oF ano Ritts By Water Erosion (42) ome 
 Gullies 2 and rills often fo form in areas of little runoff. In fact, they develop eS 
- mostly by extending their upper end uphill, where the quantity of streamflow | 
_ is smallest. The most baffling observation is that this upward advance of erosion 4 
channels is often caused by undermining, which seems to be due to = ; 
- coming out of the soil and not pouring over it. Such undermining is followed 
by caving-in and then by) a gradual transport of the debris by overland flow. 
In many instances the present writers have observed subsoil channels that ae 
to start and end where there is neither overland flow nor a water table nor 
a distinct impermeable layer, aiqzono> 
explanation and other erosion phenomena depends on the acceptance 
of two processes, : (1) A mechanism | by which outcoming water ¢ erodes - 
the ‘soil; and (2) a poser hn a by which water comes out of the soil. ak goa 
_ The first mechanism is undoubtedly that of seepage forces. At t sharply co concave i a 
points, very high hydraulic gradients can arise due to the convergence of “. 
streamlines. The drag forces exerted by the outflowing water can then detach — 


soil particles, ‘overcoming even high cohesion (40, 26). 


i 
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the soil. When it does, it must do so at a positive pressure (at least somewhat 
- higher than atmospheric), while at the same time the soil must be saturated 
In common thinking saturated soil is related to either of two cases: (1) A ] 

_ high or perched water table situated above an eeneneatibine layer; and (2) overland — 

flow, which forms whenever the rain exceeds the soil’s infiltration capacity. _ 

_ The questions to be answered by the present writers were how water outflow 
could occur in the absence of a water table or a perched water table and in 
the event of the rain’s not exceeding the infiltration capacity. gus 
... has here been Suggested that there is a lateral flow component that can 

occur with any amount of rain and in unsaturated soil. This is the horizontal 4 

flow within the soil and close to the soil surface, which can cause car 


4 accumulation in concave parts of the landscape to such an extent that it may 
even reach saturation. Saturation, of course, can be followed by outflow from 
the soil, by erosion, and by Uoati at) 


7 badly gullied valley around Nahal Bohu in the Israeli Negev was ‘subject — 
to soil conservation and afforestation. In the preparatory stage, two aerial pF 
photographs taken 20 yr apart were compared. The tips of some of the gullies 
were advancing at an average rate of about 1 m/yr, the process invariably 
Z being tunnelling followed by a cave-in. The advancing gullies were in nearly 
AG every instance in concave parts of the landscape, where the topography is " 
_amphitheater-like. The measures proposed for counteracting further erosion were: 
(1) Underground drainage; (2) flow barriers; and (3) filters. Although there was = 
no obvious water table nor any other zone of saturated soil, the drains would 
nevertheless draw off water. If this model for the formation of rills and gullies 
_is correct, i.e., by outflowing water, then any increased infiltration of water — 
- into the soil augments gullying. Turning an area into irrigated land as well as 


certain ‘‘conservation’’ measures could thus 1 mean higher ground water recharge > a 
More Osservations: Surrace Hyorovocy of Straw 
After he had presented some of these ideas at an experimental station in . 
a Ohio, the first writer was shown an amphitheaterlike drainage basin with a 
‘spring at its mouth. There was no evident impermeable layer, nor did measure- 7 
ments indicate saturated flow around or below a very small area where water a 
was seeping out. Yet the seepage continued long after the rain had ceased. 
Since then, many such places, 6... -» agricultural fields and roadcuts, have been 
observed with evidence of erosion and seepage. 7020) gi 
‘The “straw roof’? example is probably the best means to unlodge some of 
_ the old-established concepts and point to better ones. The building of straw — 
or reed covered (‘‘thatched’’) roofs is a practice shared by primitive tribes 
near the equator with some very well-to-do private houseowners in Northern © 
7 Europe. An ‘‘expert’’ hydrologist, if asked for his opinion, would first measure a 


the infiltration capacity of the straw and would then pronounce it too high | 
} a roof. But pace the expert’s opinion, not a drop of rain gets through 
_ the straw into the building’s inside. Every drop of rain is shed clear off the 
roof, not a single drop runs along its surface as ‘‘overland flow.”” There is’ 
“ground water” or “‘perched water.”” How ow confusing! 4 


| 

—— 


_ The assertion that most soils are nonisotropic with respect to water flow — 
is hardly in need of proof, since this view has been a a commonplace, at any 
x 
rate, as far ar as saturated flow is concerned. The present writers extend it, without _ 
: proof, to nonsaturated soils, and this, in turn, leads to the conclusion that 
in any soil with sloping layers there will be a flow component parallel with 
the layers. It is postulated that the horizontal flow component will be proportional 
tothe vertical component andthe slope. 
If this postulate is correct, even if only approximately, then rainwater must 
accumulate in concave parts of the landscape and at soil cuts. Another mechanism 
_ of rain accumulation in concave parts is probably that involving raindrop splashing. 
_ Many phenomena can be related to the lateral flow mechanism and to the 
accumulation of water at certain points in an otherwise nonsaturated soil. Among © 
1. Major portions ns of the rainfall are absorbed by the s soil and o only later reemerge a a 
2. This mechanism could explain v why runoff follows a certain total amount _ 
_ of rain and not necessarily rainfall of a momentarily high intensity. ee erers 
«3. It is possible to conceive of runoff with rain nowhere exceeding the 
7 4, “Partial area contribution” to runoff is made not only by less permeable 
Soils but also by soils locally waterlogged due to lateral flo lo 
5. The same waterlogged areas also contribute to the net recharge of ground ; 
= even though the average rainfall would be unable to do so. pelt thas 


6. Local accumulations of water and water seepage out Of the soil onl 


7. The rain first infiltrates the soil and only later reemerges but with delays 
that are much shorter than those met with in the case of groundwater. Absorption _ 
and reemergence can place within one rainstorm or in practically any other 
8. The concepts described can explain various leaching phenomena that were 
formerly held unaccountable, 
= 9. Also capable of explanation is the formation of a soil catena with B 
horizon, most highly developed in a concave part and least developed in a 
convex part of the landscape. spabed bo 
The foregoing s statements, if | proven correct, should replace ¢ at least partly 
the notions presently current in surface hydrology, i.e., that runoff is aneeed 
ia by rain that cannot penetrate the soil or by water, from whatever source, 
the first reaches the ground water table. Further fundamental reservations can 


OMENA 
| 
|. ase, extreme, mdicates some Of the limitations of present-day _ 
= surface hydrology. No builder in his right mind would make a straw roof flat. ' 
il 
| 
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be oe concerning the traditional vi sian that runoff is equal to rainfall minus 


‘infiltration, reservations that are unconnected with the mechanism of lateral 


: ‘The practical implications of the new concept could be far reaching. Speaking — 


ig 


or improved applications. The whole subject of watershed hydrology, and erosion 
mechanisms is due for revision, together with commonly held views about soil 
and water conservation. The Beer-Sheba experiment 1 not t only well illustrates 


it also points to an ‘immediate practical slien. ‘Obviously, even the most 
uniform irrigator cannot avoid an eventual nonuniform water distribution, unless 

_ high local curvatures are eliminated by grading. This i is certainly a surprising — 

‘Future work could usefully make a detailed study of the mechenioms causing 

lateral flow. A part of such a study will be presented in a future ire installment 

of this report, but the extremely wide scope of this subject excludes any possibility - 
ofa complete coverage of the required theoretical studies, the computational © 
methods for specific problems, experimental support, and the formulation and | 
testing of practical applications. Rather, it is hoped that this series of articles — 
will bring about a more widespread effort with these aims in view. The next 

- article will propose a theory and preliminary experimentation with | raindrop : 
splashes. Next, we shall demonstrate in general the formation of lateral flow 

_ due to gradual changes in soil properties, whether at or near the surface 4 
deeper in a layered soil. Finally, we shall use numerical solutions to study 
: 4a in detail the process « of moisture accumulation in concave parts of tl of the landscape. . 
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ols are used in this paper: 
volume of moisture per unit volume of soil; 


h= "horizontal (water flux); <a 
columns in topographical grid; and wit 


excess s vertical flow. ; 
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a Dan Zaslavsky' and Gideon Sinai’ 45, 
¥ In a previous report (6) it was postulated that there pes! four mechanisms pe 
by which a horizontal flow component is associated with rain or with vertical 
- flux in the soil. It was further asserted that that component must increase — 
with the slope of the soil and could even be proportional to the slope. It certainly aN 
.- proportional to either the vertical flux or the rain. Such a lateral flow had 
to be postulated i in n order to to explain a number of phenomena See a 


while, for the specific case of vertical rain and symmetric outflight of the splashes, ' 
an expression is obtained for the horizontal flow. It is indeed seen to be 
_ proportional to the slope, the rainfall, and an experimental coefficient which 
f _ incorporates mass ratios, the kinetic energies of the splashes, and s some geometri- 
_A preliminary experiment actually ‘showed the downstream horizontal flow 
of splashes to be proportional to the slope and of an order of magnitude as 
 Slanting rain over a rough soil is another cause for nonuniform rain distribution. a 
% ‘Together with the concentration « of splashes i in concave spots, it causes variations © 
in the effective rain distribution. The variations and their dependence on the 
roughness of the soil surface will be assessed. 
a In their turn, the variations in the effective rain distribution may lead to 
"nonlinear and nonhomogeneous relations between the rainfall average and its — 
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= nonuniform distribution of the rain or, in other words, local ‘high 
concentrations, may explain in part the formation of runoff at low average — 
= of rain and a number of other phenomena. Newly adopted one 
such as variance of the ‘rain, “roughness « of the soil surface, the azimuth, and 
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Way Do Rainorop Propuce fide LATERAL Fiow? 


Raindrops that hit the soil break up into smaller droplets (‘ ‘splashes ), which — 
rebound (2,3,4). If the soil slopes, the splashes will travel farther downhill than 
uphill. The virtual center of gravity of the array of droplets resulting from 


es” 
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FIG. 1.—Typical Splash | Shape Caused by Raindrop Hitting Soil Covered with Water 


the original raindrop | will be found downhill from the first point of impact. 
‘This means that, given a certain distribution of rain intensity at some horizontal 
surface above the soil, the eventual “‘effective’’ rainfall distribution over the 
surface of the soil will be the result of a downhill translation. The horizontal 
- discharge amounts to the rate of rainfall multiplied by the horizontal shift of 
the center of gravity of the raindrops. An observer v watching the passing ng splashes 
will count more going downhill than uphill. The net difference amounts to a_ 
ey real net lateral flow, which can be demonstrated experimentally. At least — 
for moderate slopes we may state that the downhill translation increases with : 
the slope. Over a long and uniform slope the resultant horizontal flow ‘hom 


“not be 1 ‘noticable and the final rain distribution will remain unchanged, whereas 
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at the top of the slope, at its bottom, and at any point of duag in ‘the. ape, 
_ the effects of the horizontal flow will be felt, eS 
.. the following we shall present three more or less independent eatin, 
: (1) An experimental qvelention of the net downhill flow as a function — 


a raindrop hitting a water surface will produce a crown 1 of ‘Splashes (Fig. 4 


1), a phenomenon that has" been investigated b ee 
= splashes” arise at a fixed angle, | B, with the horizon. ‘The’ configuration of the 


o splashes and the distances travelled by each splash have been found to be 
Bb around the initial flight path; but the experiments have concerned 
only vertical flight. Splashes also occur on a nonsaturated soil ‘rere ; but bes 
we know no fort this case. ~ 


yaticred 


q on wo La) 
_ The present writers fir first adopted the. convenient assumption that the exit 

"ange B is uniform around the drop and that the exit velocity v, is more or. 
less uniform. It is noteworthy that neither the final solution nor the conclusions ’ 
ate sensitive to minor deviations from these assumptions. te 
one: _A single droplet having an initial velocity of v, , and departing at an angle, 
B will be assumed to describe a parabolic path (with no air resistance and vr 


over a flat gravity fi field). The components of its velocity are then - on 


’ 
4 


ie vertical ‘coordinate, positive » and upward aa = 0 at the initial 
horizontal coordinate in the ren of flight of a splash © 


_ hitting points); x= 
= gravity acceleration; 


at the initial hitting ‘Point of the raindrop); g 


i 
4 
which z 


and x,, x, = sumnatode the downstream and upstream distances from x = 


0 of the the points: of impact of (Fig. 2). If the of 
x cos*B (tanB +tana) . 


‘tana 


Consider now a three-dimensional picture (x, y horizontal coordinates and a 


2 upward vertical coordinate with the origin at the point of ‘The : 


4 


FIG. 3. of Splash. sh Range (r, yo on n Slope (Plane Vi View): — is for: Slope 
a = 20° exit angle 50° 

— of splashing, r,, of a mass, m,, is simply obtained from Eq. 4 and Eq. 5 by 

the slope (tana)to(tanacos®) = 

V2 cos 7B (tang + tana (6) 
f 

Each droplet at an angle 6, a on the average, a -Conjugat ea 


+ 7. The. 


the aver te at 


4 
| 
q | 
Ar, tana cos6,.......». 8) 
The distance component in the direction of the slope (X-axis) 1s: 
Ax, = 7,008 0, = — Vi cos’ 


"ghee m} represents the mass of a while m = = the mass of the | 
raindrop; and m’ = the total mass of splinters caused by the raindrop. The 
full expression for the shift, based on the sine assumption, is obtained 


many drops (in time and space) we Teach the limes, i.e., the 


1 
‘Note that the change in the center = gravity has been enti for the 
mass of droplets m’ and not for the mass of a single raindrop. The mee 
Be of Eq. 11 is easily found if B and vz are ¢ considered to be independent, | on G 
average, of @ and if m/ is a time average and independent 


fe =) tana = 2 sin 6, tan~ tan B, 


rey More generally, we obtain from Eq. 6, using, me definition of the er. 

2 cos 8 (tan B, + a cos @, ,)cos®, 


in which the summation over splashes in and space. It 
is divided into two terms. . The first (2/g) =m} v2. cos’ B, tan B, cos 8, would — 
Vanish if V, m and Ba are | independent of 8. If they are not, , the probability 


7 


average shift ¢ by 
| 
| 
d 
is then 
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B, 1 =m’ cos 29, mi cos’, sin? By 


16a an to 1/2 when ‘the m, distribution is independent 
6,. The term (Eq. is the weighted average of the vertical kinetic 
. aa the splashes. It expresses the maximum height they reach in their flight. 
The last term (Eq. 16c) is the weighted average of tan” 
; Pati in words, the shift of the center of gravity of the raindrop due to splashing a 
a is proportional to the specific kinetic energy of the splashes V2 /2g or, in other ee 
_ words, to the maximum possible height to which these splashes can jump. This _ 
kinetic energy is probably related in some manner to that of the originaldownward _ 
} velocity less some losses due to friction due to the pick-up of soil particles J 
and to the production of new water surfaces. Consideration of the momentum = ' 
“aeration requires the following two be (for a raindrop 


cos 8, sin 


i, i, having a a mass sm’. At the limes of many drops, a ‘symmetry y of the horizontal + 
_ momentum with respect to two orthogonal lines must be maintained. In addition, J 
7 the probability function of V, m, and B, with respect to any angle, 6, must § 
_ ibe continuous. Furthermore, the probability distribution function should = ; 
= around the slope ‘direction. There is a very small number of 
possibilities for the velocity and the mass distribution of the droplets | around — 


with a translation Ax’ The resultant lateral flow can be found by counting — 
the number of drops passing through a vertical control surface. This is obviously 
- the amount of splinters per unit time, multiplied by the distance upstream (Ax’ ; 


horizontal flow, Q,, is then simply 


| 
% 
| 
| 
7 (185) 
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which Az’ = the time and area average of ‘the: splashes; P= 
a the rate of rain; =m? / am, = € = the ratio between the splashing mass and 
the mass rain, and = = = the of 


SZ the rain, and on the type of soil, but it also depends on the total amount — . 
of rain and its intensity, since they determine the wetness conditions at the 
soil surface. The ratio between the mass of ‘ i and the mass of 


In the expression for Q. in Eq. let us introduce ‘the explicit values 


introducing a new parameter, ws the maximum height of the flight trajectory a 


in which V,, and B, = respectively the s = speed and the of exit- velocity 


vector, Eq. 19 will become, with cos’ 6, = 


a _ The parameter 5 may be estimated on the st strength of indirect een 7 
By looking at the mud stains on walls, the height 8 can be observed qualitatively. 4 

A more accurate value may be found by measuring the drop density in the — 
air or with the aid of blotting paper stained with methyl blue. From the =_ 
‘Stained blue the mass of the drops at every height can be estimated. on 9m 

- To sum up, the horizontal flow in the direction of the slope is _ ina 


wy aly oinz. 
22 is ‘similar to in the first part of this report 
(6). Due to the extremely short delay and small storage inherent in the pn 
phenomenon, eq. 22 may be considered to represent a quasi-steady state. The 
4g . parameters 5 and U are due to the rain energy and soil condition. Their form 7 
recalls a layer thickness and coefficient of anisotropy. salt vhesgll 
-Preuminary Measurements of Laterat wedesiqe 10. 
quantitative experiment was made with impermeable and relatively 
es on various slopes. Rain was portent by a raintower 17 = J 


a | 
| 


_ high, in which the ° average flight of the e drops were vertical. The rain was reasonably a : 
- jloanily distributed in time and space, and the flight velocity closely approxi- 


mated the final velocities encountered in the _ air. Measurements were made ‘i 


downstream edges of a wide send: Miuentictinn: a considerable number of © 4 
fell beyond the sides of the slope. Thus the absolute values of lateral a 


4 


& measured are rather low and could in real life be at least 20%-30% higher. : 


From the measurements taken, actual mass moments could be deduced and 

_ changes in centers of gravity calculated. However, a simpler estimate of the 

_ net lateral flow could be obtained by subtracting the net downstream splash © 
‘the late from the net upstream splash discharge. The ‘difference is simply — 


lateral @.- In the results are as a function the 


ransiati 


Slope tan x O 0035 O105 0194 O287 
4. Q. due t to Splashing Rain of Intensity over an imperme- 
able, Smooth Surface of Slope 


© There are two conclusions that can be drawn from the measurements: w 
A considerable lateral flow occurs due to raindrop splashing; and (2) within — 
_ the range of our experiments, the lateral flow due to raindrop splashing i rapsnecnndl 
monotonically with the slope, while with slopes of up to some 20% it increases — 
An idea of the order of magnitude of AX can be gained by observing ‘the 
i height of raindrop splashes on vertical walls. Typically it reaches 30 cm. In y 
ss 22 this can be taken to be representative of 8. A typical value for B is 
45°, so. that tan? B = 1 — 1.4. Let us also assume that 2e = 2; we then 
have U U=14— 4- 2, so that = — 


| 
- 
| — 
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ag tana to 601 tana, in centimeters . 

experiment reported in Fig. 4 (up to a t of 20%) a rather 


_ the associated mass thrown up by the splashing eet is larger, so that 
€ > I, or both. The measured and the estimated figures are very similar. The — 
actual figure will probably vary with the roughness, the aggregate strength of 
the soil, and the soil moisture content. Similarly, it may change with the intensity 
of the rainfall, not only due to changes in the moisture regime at the soil surface, 
_ but also through the increase in the specific kinetic energy of the raindrops — 
ee with increased rain intensity (brought about by an actual increase 


Possiste Ratio Between Laterat Frow and Rain 
ah 


a uniform ‘slope the contribution of t the | lateral to 
may be since it is a constant that does not depend on the : slope’ 


4 lateral and the rain discharge i is pa) 


‘The parameter Ax in the sdiealile experiment was found to be neat. 
tan (in centimeters). The ratio ‘thus depends on the steepness and the 
length of the slope. Assuming tills, gullies, and depressions to occur every | 
2m, with a ridge in between, a typical value of the slope length from a ridge 
to the nearest gully would then be L = | m. ‘Further, assuming an inclination 
; not exceeding 10%, Ax/L can reach 6. 6%, i.e., , 6.6% of the rain flows towards 7 
the depression due to splashes alone. _ This. is a considerable amount of of lateral 
The term (tan /L has an important physical meaning, the | 
of rain in concave places being obviously proportional to this ratio. It is the 
__ slope times the drainage density of the landscape. Geometrically, it is the curvature <a 
of the landscape. In order to simplify understanding, it is useful to ce , 


a esa of a soil surface ‘sinusoidal in | shape, the elevation, 2, , of which is 


The parameter A/2 = the amplitude of the sine wave; 
length . The first gee of z is the slope. Its maximum value is (7A /2L), 
average value, 5 = The curvature is arri at the of the 


4 
| 
4 
| 


i terms of the average slope 5, it is s(n? 5/2L 
Excess Rain in Concave 
In the foregoing it “that the lateral flow due to raindrop 
"splashing may be responsible for the accumulation of rainwater in concave | 
locations. A more rigorous expression for e excess rain may be obtained be a 
Eq. 22. _ Consider the slope, tan a, to be: a vector senting from two c components — 
in the x and y directions, with z = the elevation. Then 


the horizontal flow in the 3 x, y plane, is also a vector parallel with ta 


aie > assumption, of course, that ‘@en rainfall P as P we the coefficients, U and 
| (Eq. 22), are independent of the coordinates ‘x and y. If they are not, .- 
_term should be added to Eq. 30, i.e., grad (z) + grad (US). The total —, 
of rainwater landing on a soil would, ¢ to 30, 
: which V*z = the curvature which can be measured geodetically [see 


_ The orders of magnitude of the ‘components of the term in brackets of Eq. 
31 are of some interest. It has been shown that the term, U- 8, can reach 
a value of 66 cm (at least in the aforementioned experiment). For slopes of 
ok m length and elevation differences of only 0.1 m, the curvature is the order - 
of 5 x 107° cm~'. Thus the excess rain at such concave points can be 33/100. 
‘The effective total can then be nearly 1.5 times the original rain. os. ~ 
a _ A precipitation higher in certain locations than the average rainfall constitutes 4 ‘4 
an extremely significant topographic effect from both the hydrological and the 
_ ‘The curvature at some points in the field can be very high. ‘Theoretically, 
the curvature at the intersection “of two plane slopes (second derivative 
elevation) tends to infinity. This would appear to mean that an infinitely large 
_ effective amount of rainwater should be expected there; a problem that hes 
Pd yet to be studied in detail. What is certain is that there is a scale effect, which | 
_ depends on the typical distances of splash flight and on the minimum wavelength ~ 
be into account in in describing the shape of the soil, assuming that 


The excess rainwat¢ be obtained through the conservation 
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a place of infinite curvature owe ast ta in any way contradict the physical reality, 

‘since it i is confined to an infinitestimal area. The writers have conducted 

experiments with a V-shaped, symmetrical slope, and runoff and erosion appeared 

almost immediately at the bottom edge. However, this phenomenon may also — 

be “pes by ‘Secondary splashing and by the flow in the surface transition | 

— ‘The problem of raindrop a action on a the 1 soil s structure will not be fully treated 

an here. For the present analysis it is sufficient to note that the raindrop splashes — 
“carry with them soil material the amounts of which are measurable. Typical 
figures obtained by the writers showed that, in a Loess soil, there was 1%— 10% 


by weight of splashed soil in the rainwater. _ Typically, treatment of the soil * 


A 
ai A oi ton of? 29 
bas FIG. '5.—Boundaries of Surface Drainage Basin 4 


- reduced this figure * 0. .25%-0. 5%, in 1 some cases it at least st reduced the fi flight 
distance of splashes. Each millimeter of rain in the previously calculated example, ‘ 
_ of 1 m slope with fluctuations in elevation of 10 em, produced 70 g of raindrops _ 
flowing towards the | depression. The 1%-10% splashed soil means 0.7 g-7 3 
of soil/m’ splashed towards the depression which is a considerable amount 
‘It can explain the accumulation of splashed material in soil depressions. One — 
can see such splashed material on almost every soil after any rain. It can explain 
the appreciable erosion caused once there is an actual runoff coming out “dl 

_a depression and capable of carrying away the ¢ splashed sol. disch 
gy A typical annual rainfall of 500 mm can move 0.3 mm-3.5 mm of soil downhill. 
in the preceding example and ones more than 10 times this amount. doistw a 


q q 
= 
| 


gs Eq. 31 may be integrated over an entire field. For rainwater to be conserved, 
condition must be fulfilled: te vistaiborman 
4 aban ag ba 
| 
‘in which We= mig finite loss of, rain | across 
the boundaries that becomes relatively 1 negligible for pd enough areas. ene ; 
first integration, Eq. 33 reads (Fig. 5) 


a2 


in which the term P - + UB [(@z/ax), — (az/ax), ] = the net addition of rainwater — 
to a strip of unit width in the y due to splashing in the: x 


2 


J 


- 6. into (5); Rain Angle B, Soil ‘Slope a 
The order term defines the net gain in the y direction. A drainage basin is 
defined by a boundary along which the slope normal to the boundary is zero. 
Therefore, over the entire surface of a drainage ope ys 0. There must 


_ The accumulation of ‘rainwater j in concave places is an important cause in 
_ the nonuniformity of the final distribution of the precipitation. There is another : 
t form which has been published in a report before (5). Consider a slanting | rain _ 


(Fig. 6) falling ; at an angle B to the vertical over soil sloping at an angle a, 
the azimuth of the slope veing y and that of the rain, w. The effective rain 
(based on simple geometrical considerations) is then 


| 

| 
| | 
a | 
| 

| 
| 
| 
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It can be shown that 1 7 tan a can easily be +1, leading to a virtual doubling» 
‘ of the effective 1 rain on one slope and i its reduction | to zero on the other. T he 


, &.g., y and 
w have ‘the sat same peters i.e., , the: two-dimensional case where the slope and 


the rainflight have the same aspect. Eq. ds then becomes = ~—— 
= P(l + tan a tan B).. ob GB 


4 When a = B = 45°, Pio sih is either 2P or zero. For a vertical wall a finite 


‘Tate | rain on an infinitesimal point (tan a - > It would be 

‘interesting to combine the two mechanisms of concentration in concave areas — 

(Eq. 31) and slanting rain (Eq. 35). 

_ A further complication of the rain behavior at ground level may be expected & 

as a result of variations in the intensity of the rain, P, high upinthe air, = 

Cuances in 5 U To Be Expecteo Due To Cuances in Rain Data wad § 

_ Formulating a law for the lateral flow due to splashing would be ; simple enough 


if the anisotropy U and the jump "na 3 were ain of the rain. However, 


ot bas gait 


} 
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Bic: 
-— is probable that the ratio of the mass splashes to that of the ‘total rainfall 7 
as well as the recoverable kinetic energy both change with the rate of oo 
rainfall. Investigations (1,2) indicate that larger raindrops are formed as | the = 
intensity | of the storm increases. Such larger drops certainly attain higher final 
velocities. It has been observed (2) that, with a given soil, the specific kinetic — 
_ energy of the splashes is proportional to the original kinetic energy of the nin 
drops. The product U- 3 is likely to be proportional to the recoverable kinetic y 
0) energy of the splashes, and it can therefore be expected that the dependence 
in U- 5(P), is sit similar to that shown in The horizontal discharge, 


ba Q,, will therefore be related to the rain intensity, P, by an expression waned 


| 
| 
| 
4 


corrects Eq. 22 in which U-d was considered to be a constant. 
xperimental proof of Eq. 38 has been attempted, but it is reasonable to assume — 
_ This in turn indicates that the average, UB, would be the sum of their values 
at average rain intensity, P, and a correction term proportional to the variance 
of the rain over the area. This means that - fluctuations in rain intensity | lead 
to an additional contribution of the lateral flow pro ortional to the rain’s variance. , 
7 
It also means that short and strong bursts of rain, even localized ones, can ~ 
_ produce heavy lateral flows out of proportion to the rainfall. Higher intensities | at 
will be associated with more extreme concentrations in concave places. im A A atl 


ation, Runorr, ano Grounp Water RECHARGE 


en 


the main conclusion to be drawn from the foregoing analysis is that ev 
ander rain waiform in the air, the effective rain can at some points on the 


Furthermore, a rain nonuniform in time and space is liable to increase “the 
amplitudes of the fluctuations in the precipitation over the field, 
_ Let us assume that at least in some cases the local precipitation can eventually 

exceed the infiltration capacity and that the excess will form runoff. The same 
model may, incidentally, apply to ground water recharge, which is the difference 4 
between infiltration and evaporation. The mathematical form of the 

= — will conform to a difference equation 
R= P- > 

7 in which R = some net effect (runoff, recharge, or even crop yield). A common 


error is to write the average, R between the the averages 


T, thus 


“the; averages P and Thaving been measured and then calculating R. More correctly 


Earlier on it was demonstrated that P ee C4n fluctuate widely between zero and i 


a value several times P. If has an expected value, 


and /, respectively, T and then as a first approximation (neglecting 


1 = the correlation between flictuations in P and I. In most cases 
_ such a correlation exists, e.g., when / is surface retention and R is either runoff 
r water recharge, the correlation is positive. Surface crust causes rain to 
accumulate, especially i in concave places. Thus, when / is infiltration, a negative 


It may happen that R = 0, ie., that the average rain equals th the average a 
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infiltration (plus retention, ete.).. Runoff may occur due to a local 


concentration of rain as expressed by the statistical terms of Eq. 42. Similarly, — 
very slight rains may effect a net water recharge in some concave places _ : 
- the landscape due to lateral flow. Eq. 42 has been cezlculated (5) assuming “a 
; an normal distribution for | both P and J. It can be shown that ¢ the he result i is bes “ae 


—— ——efe .. (4a) 
Pp and I must not necessarily have a normal diettibalion. nor, for hat matter, 
-. they, if only because P = 0, J = 0. However, other distributions, while 
"mathematically more complicated, will at any rate qualitatively lead to similar 
conclusions. It is interesting to study the shape of E of Eq. different nt ales 


of Pando. ForR=0Eq.43areducesto 


q 
_ Fluctuations of P can be of the same order of magnitude as the average ff 
value of P, and | even larger. The value of o, can, therefore, be about P/2- 
or larger, while the runoff or net recharge can still take up a significant part 
of the total rainfall (1/4 and more). This is entirely due to local —— 
= of rainwater thanks either to slanting rain over a rough surface or to lateral 
flow caused by splashes. dw ais) lo sft 
_ At very high values of P (and of R) the second term in Eq. 43a vanishes f 
and R 4; increases proportionally to R. Although this is true only if the = 
deviation, o, 1 Oo, Temains unchanged, it has been found that the fluctuations in J 
- effective P increase with the average value of P (Eqs. 31 and 35). Therefore, 
_ the ratio (R/o,) in Eq. 43 changes only as long as J remains a constant. Thus © 
probably all three terms in Eq. 43 increase in a similar manner. Contrary to - 
what has been suggested before (5), the term denoting local concentration of 
_ precipitation does not become negligible with increased rain intensity. ‘Iti is even 
~ possible that the effect of the fluctuation increases as a result, e.g., in view 
_ Amost interesting conclusion presents itself here. The simplest | linear equation 
(Eq. 29) has led, upon averaging, to a nonlinear relation, i.e., Eq. 43. This 
‘is due to a combination of fluctuations and a threshold J. Eq. 43 secails 
ae relations between rain and runoff that are far from linear with — 
R SP tiga dip change Aa Segre, hes be 


In most hydrological models the slope of a field is considered to be an important 


- feature. In the foregoing it has been proved that the local curvature is is just: 
as important a feature, at least as far as raindrop splashing is concerned. 7 pass. 
_ The specific mode of averaging depends on the explicit formulation of he , 


end function. As an example, the oldest and most common hydrological model _ 


q | 


is exemplified i in n Eq. 40 and so corrected is in Eq. 42. In Eq. 43, a 


Statistical distribution is assumed. The variance of P has been related to variations q : 
3 Let us calculate the variance of P, taking o? from Eq. 31, and ‘assuming — 
U and 8 to be constant concentration only in concave spots) ~~ b issvon 


2 
op, = P?U 


dy 


Another variance is due rain (Eq. 35). 


Note that y changes from one - point to th the other as the aspect of the ‘slope 
changes. This means that the roughness of the: field, as as defined in. Eq. 46, 
is not a contant topographical parameter but may vary with changes in the 
aspect of the rain. Still more complicated cases can be constructed by negro al 4 ‘a 
7 the two mechanisms of slanting rain and concentration in concave spots. So 
3 far, only integration over space, as in Eqs. 45 and 46, has been carried out. . 
Fluctuations of the rain intensity P in space : and time and the dependence of 
4 3 on P would require the calculation of yet another variance. ae 
It becomes important at this point to consider the need for measuring different 
topographic and rain parameters that have always been neglected in the past. 
¢ 
ax parameters are bound to have a decisive effect on both the local and — 
the temporal concentrations of rain, which can cause ground water recharge, | 
runoff, or at least different degrees « of wetness of different parts of the soil. ar 
_ Asan illustration let us imagine a two-dimensional sinusoidal landscape. With 
7 the elevation z fluctuating around the average Z at an amplitude A ,/2 and half 
cycle x= L, aU a given soil wave), we have 
\ 
Sl 


the curvature (dispensing with the index as one cycle only being considered) 


_ The roughness or variance according to Eq. 45 is found from om tre ol. sales 


w(4 


PU sin 


| 


5/2 
As mentioned earlier, the term A’ 2/L 2 denotes the maximum curvature 


is the e drainage density. } For L ~im m and A only 0.1 m, ‘UB can be 0.50 50 
m (as. measured). The standard deviation will be around 0.2P. 
_ A second example is that of slanting rain (Eq. 46). Assuming 1 = = I and 


rain fall at an angle of 45° 


with (A/L) = 0.1 the value is about 0.18 P. Much higher values are obtained 
_ when the relief of the soil is steeper and the angle of the slanting rain shallower. 
.. A rough estimate of the final variance of the precipitation on the soil surface + q 
= be gained by adding up the separate variances. The total standard deviation _ 
can thus easily amount to 30% of the average rain intensity even if the rain a 
- originally perfectly uniform in time and space at some height above the soil 


Summary ano CoNcLusiONS 
il theory has been stool relating the horizontal flow to rain splashes. 4 
Relying on certain assumptions that flow was found to be Proportional to the 4 
slope and to the rainfall itself as well as to a given fraction of the kinetic — 
energy of the splashes. A rough estimate of the incalculable parameters showed 7 
them to be of the same order as the measured ones. Measuring the horizontal — 
shift of the center of gravity due to splashing is sctalanely simple, and te 
— validate the Proportionality | to ) the Slope up to an inclination of 20%. 


Local concentrations of rain in splashes occur in concave parts of the landscape, 

- concentrations of rainwater of up to 50% of the average precipitation over a 
___ the entire area having been shown to be possible. The measure of the concentra- _ 
_ tions is proportional to the curvature of the soil surface, but the behavior at | 
4 very high curvatures, i.e., a sharp change i in slope, has still to be investigated. a 

Rain concentration there is probably moderate. 

ier _ Another cause for higher local effective rain is the interaction between rain 

- slanting (the angle at which the rain falls) and surface roughness. This interaction - 

can easily cause a doubling of the rain over some parts of the surface. The 

combination of concentration in concave spots and over some slopes facing” 

the slanting rain produces on the whole a rain variance of a very high value, | 
and even if the most traditional view, i.e., that runoff is rain minus infiltration, _ 
is adopted, it is seen ‘that the rain can, locally and temporarily, be much higher z 


while ation at the soil surfaceis 
| 


than the average. The rain can thus produce waterlogging and runoff even ae 
its average does not exceed the soil’s infiltration capacity. Ins a similar _ 
local ground water recharge can « 
infiltration does not exceed evaporation. A suitable way ‘to average “these 
- parameters has been described. It is very interesting to note that, taken on - 
average, a basically linear relationship produces a strongly nonlinear one much — ; 
like the relations found between rain and runoff. 
he The conclusions based on the foregoing analyses do not only make for a 7 
better understanding of effective ‘Tain distribution. They clearly point to the 
necessity to measure hitherto neglected parameters that affect surface hydrology. —_ 
They are the variances of the rain, infiltration, and whatever other parameters _ 
eater into the equation. Furthermore, the character and the slant of the rain 
must be related to one type of soil roughness parameter which involves the 
a 
square of the slope and its aspect to ‘the rain. The other type of "Toughness: 
parameter is independent of direction and is found to be the square of the — 
: local soil surface curvature. Both types of roughness depend on the details — 
in which the enrveying of a soil surface and the calculation are performed. — 
a ‘They are not unique. Just how the roughness is to be calculated and how it 
r _ affects the hydrology of a watershed should be the subject of future studies. — 
tt may be relevant to future agrotechnical practices, to the characterization 
of watersheds, and to the choice of prototype watersheds to be studied. ee 
The calculation of horizontal flow with the aid of rain splashes enables the 
discharge of soil eroded by splashes 
Inc., 
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following s symbols are used in this paper: 
twice amplitude of soil surface elevation; on 
infiltration, surface retention or any entity subtractable | from original 
= mass; indy, com sdy asve bas 


= 
® 


rate of 


net ef fect—rain minus infiltration, minus 
“coefficient of 
velocity of raindrops, of rain splashes; 
volume of rain splashing across border; 
horizontal coordinate, distance of point of impact with ith subscript 
horizontal coordinate (normal to x and to Loe 
elevation or vertical coordinate, positive upward; 
= soil surface slope angle (inclination); 
= at between flight of outgoing splash and horizon, or r angle between 
raindrop flight and vertical; 
azimuth angle of slope; 
maximum height of flight of rain splashes; J 
ratio between mass of splashes ant 


o= exit velocity of splash: tarefylly Wy some 


= upward or upstream; hur coworkees 

effective; and Nest by of the 

an rec waiter citer ty Dow whe andl 
the channel by of wales fiom farther .- 


= of total splashes belonging one and 
= excess (of rainor moisture). Morel, 
Expand fren (be upaand apd outward” 
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tad) on it 8) le tsiweH vd. benoiigen 

— the first part of this report (15) we stipulated the existence of certain 
which rain and vertical infiltration do not penetrate ‘straight 


down into the soil. They induce a horizontal flow component which is proportional — a 


to the slope, the rainfall, and some coefficient of anisotropy. . Such mechanisms _ .- 


could explain a number of phenomena i in hydrology, erosion , and soi soil- formation — 


~ have been made in the past and recorded fatefully by some investigators. The - i. 
Seaeiee by J. D. Hewlett and his coworkers (8,9) should be mentioned — 
in particular. They reported two fundamental observations. One is: “There is a 
a doubt that the largest component of stormflow from forests and ‘most 
-wildlands begins as s subsurface flow. Sometimes. the term. interflow is used 
synonymously with subsurface stormflow, but other times it is used to — 
any nonvertical subsurface flow above the water table.’ The second observation 
is that: ‘‘Channel precipitation is augmented first by Snenelen of the surface 4 
_ and subsurface areas contributing water either by flow through the soil pores — 
into the ‘channel or by displacement o! of stored water from further upslope . le 
=| The area adjacent to the stream channel is the main source of stormflow because 
ue the lower slope continually receives soils and ground water from upslope; thus 
the soil near the stream is wetter at the beginning of a storm, causing | 
Source area to expand from the channel upward and outward’’ (see Ref. 8, | 


4 ‘Prof. of Soil and Water Engrg., Faculty of Agricultural Engrg., Soils and Fertilizers — 
Diy., Technion—Israel Inst. of Tech., Technion City, Haifa, 32000, Israel, 
Lect., Faculty of Agricultural Engrg. Soils and Fertilizers Technion—Israel 
Inst. of Tech., Technion City, Haifa, 32000, Israel. 
Note.—Discussion open until June 1, 1981. Separate should be submitted 
: for the individual papers in this symposium. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
" ASCE. Manuscript was submitted for review for possible publication on March 23, 1979. — 
This paper is part of the Journal of the a aa Division, yy of -— American 
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p. 92-95). They named this phenomenon the “Variable Source Area. an —_ ‘a 
_ It is the purpose of the present report to explain these undeniable phenomena, _- 
‘and to formulate them inarigorousmanner, 
a In the second part of the report (16) the mechanism of raindrop splashes 
_ was shown to be in | good ac accord with the aforementioned stipulation. ‘It produces 
a real horizontal flow component which is s proportional to the rainfall, the slope, 
a simulated thickness of the splashings layer, and a certain coefficient. The 
thickness and the coefficient appear to increase with the rainfall intensity. The a 
Pe of raindrop splashes can explain how rain accumulates in concave 
spots. Slanting rainfall hitting slopes that face it is another mechanism that © 
can explain local rain excess. These local rain concentrations can explain why a 
flow in soil depressions can commence long before the average rainfall exceeds 2 a 
the average infiltration capacity. However, this mechanism of rain redistribution ) 
- fails to explain several phenomena that call for horizontal flow within the soil. 
_ The splashing of Taindrops cannot account for soil moisture accumulation over 


a 
- greater distances, such as encountered in the Beer-Sheba experiment (15) or | 


_ mentioned by Hewlett et al. (8,9). It is also. a very brief phenomenon that 
- ceases a fraction of a second after the rain stops. In contrast, some processes’ 


: of accumulation take days or weeks. Erosion due to seepage forces must involve 


rain that first enters the soil and eventually seeps out. 
In the present part of the report we shall co consider in principle he mechanism 
that can produce a downstream flow c component more or less parallel to the 
_ soil surface and within the soil. In most general terms, the arguments supporting | 
_ such a mechanism are very simple. Let us consider a soil which is not uniform, 
its nonuniformity taking the form of parallel layers. It may be said that this — 
soil’s properties have a rotational symmetry designated by some vector di, 
and the soil is on the whole ; anisotropic, one o of the principal axes: being _ 


If the force ‘motivating the flow. is directed aw: away from X,, the flux. will also 


be directed away, but not at the same angle. Consider the soil layers as being 


in a slope forming an angle a with the horizon; A, is therefore at an angle 
of a to the vertical. The force of gravity thus forms an angle a with A, and . 
points downstream. Infiltration will therefore not proceed vertically downwards a > 
will have a horizontal component downstream. 


In the following we shall deal in greater detail with the mechanism of horizontal 


fl flow within the soil. Its order of magnitude will be assessed and in some cases | 


ea we shall follow the | phenomenon through the tae of sti streamlines 


i very important specific types of nonuniformities i in ‘the soil are analyzed: I: 
(1) Cyclic layered soil; and (2) soil surface air transition layer. Here they — 
considered only qualitatively. They will be the subject o of 2 a detailed 1 quantitative © 


7 When a streamline moves from one medium of conductivity K, to another 
‘ of conductivity K,, it will form different angles, ya and Yas" , with the line normal — 
to the interface (1) so that (Fig. 1) 
we te interface (1) so | ) 


q 
| 
| 
| 
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ceva 


oF Fig. 1 can be considered an example representing two soil layers whose interfac 
_ forms 2 an angle a with the horizon. On passing from a less permeable to a 
‘more permeable layer, the ‘streamline will be diverted from a normal direction 


toward a diagonal direction, thus producing a parallel component. If, as in 


anit ES 


FIG. 1 ‘of Streamlines on Crossing Intertayer Surface: Hydraulic Con- 
this case, initially y, = a, the | ratio between th the horizontal and vertical | fluxes 


as can be: by simple trigonometric considerations. 
small slopes a), (tan’ a) K,/K, << 1, so that" 
a. 
or orl/2s sin 2a ~ ~ sin a and s ‘sin? << 1, sothat we 


| 
| 


3: In general cases of unsaturated K, 
and K, being the weighted averages of the horizontal and vertical couductivides, 

J __ If q, happens to be the rate of infiltration J, which is equal to some steady = 
_ rainfall, then K, = J. By substituting J in Eq. 2 or 3, we obtain the horizontal _ 

flow component as a function of the steady rainfall or the steady rate of infiltration 
and the slope. ° This is in full conformity | with our previous stipulation (15). ae ity 
In a medium with a gradually \ varying conductivity, the change in the flow 
‘direction will be governed by exactly the same rules. The angles of the streamline 7 
are then related to the direction of the vector, grad K, and a change in tan ‘ 
_ y is expressed by its scalar product with an elementary baal length, ds, along — 


and therefore on integration at any point 


By simple inspection of Fig. 1 and Eg. 5, tea arrive at the following conclusions om 


is regard to a uniform slope with more or less vertical entry of the water: . 
1. If the hydraulic conductivity increases with depth, the streamlines will oa 
Carve downstream to form a horizontal flow component by virtue of the fact 
- the right-hand side (r.h. ‘s) of Eq. Sa is positive in a downward direction. ag 

Note that K can increase not only as a result of a change in soil but also 
gs 2. If this conductivity decreases with depth, the streamline may turn upstream 

_ but not beyond a direction normal to the soil surface, i.e., tan y can vanish © 

but cannot change sign. Thus, there will always be a flow component parallel 

- either to the soil surface or to the interface between the layers. This component 
can decrease to q, sin a during steady flow in a uniform soil and to zero _ 
in a highly impermeable part of the soil, i.e., K << q,; tan y — 0. The parallel 

_ 3. The fact that the parallel flow component is ‘Proportional to to the slope 


the concentration of moisture in concave spots. 


ais seems mathematically simpler to use the components q, and q, in Fig. ; 
-* respectively parallel and normal to the interfaces. From Eq. 1 and some — 


in the case where is vertical at the point of entry and y, "= a ‘is the : 


40 

= 
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dn is positive ve downward), and re simple 1 result is 


: From Eqs. 6 and 7 one can review the foregoing conclusions more simply. 
; 4 The parameter g, > 0 only when K > 0 or a — 0. The parallel flow component 
exactly proportional to the slope. 
as The simplest case of curving streamlines is that of steady accretion at a 
_ phreatic surface due to rain (or loss therefrom due to evaporation). If the soil — 
above the phreatic surface is sufficiently thick and uniform, the streamlines — 
will enter it vertically. Moreover, at the point of entry, K, = qi | (Fig. 1) (for 
proof see Eq. 14 with z—> Saturated hydraulic conductivity K ,, a is observed 
_ near the phreatic surface. The slope is that of the water table a = ¥,. The 
‘ ‘streamlines will enter the water table at an an 


at an a ngle to the ‘normal, which, fo 
to Eq. 7, is | 


4% 


or at an _ to 0 the vertical i in the xz coordinates, ite: according to Eq. 


= 


— 
| qua — 


+ — tan’? a 


. 

4 It is interesting to note that in the case of smaller vertical fluxes the horizontal _ 4 
component at the water table will be relatively larger. One can actually calculate — “ 
the horizontal projection of a streamline as it descends from the soil surface. 
: If one assumes a saturation surface at elevation zero, one can derive the equation - 
the pressure head from Darcy’s equation ‘Sher 


-Telationship between the conductivity K and the pressure head we integrate 
Il, from Egs. 10 and 11 we assume 


| q 
| 


1981 


in a = and K = the hydraulic 
s Introducing this expression into Eq. 11 and integrating, we obtain the distribution 


in which Zo = the elevation of face at b= 
.. 6 in which K,, is replaced by any K bees Eq. 


+ exp [— cosa tana a dn “ay, 


— exp —a(n cosa — (16) 


‘at and case of a large | term n cos a — Z, Eq. 16 is approximated b 


by: 


“As =ntana+ 


acos a 


The first term is simply the deviation of the vertical from the normal. 
second term is the deviation of the streamline from the vertical. Thus the horizontal _ 
component of the streamline over a sufficiently great depthis 


It is independent of the depth. In many drainage problems the slope is between | = 
0.01 and 0. 1, and the — conductivity Prange — 10 mm/day-100 

at mm/day-3 | y. The coefficient — 
‘acan havea a typical value of 0. 01 cm ~'. Itis evident that the helpeuieh projection : 

_ &x can range from a few centimeters to many meters as the water moves 
from the soil surface to the water table. This result is particularly striking when | 
one considers certain shallow water tables of waterways and drains. Even in 
the case of relatively flat land, the groundwater slope increases towards the se i 
drain and the lateral flow component within the unsaturated soil becomes 

_ This special case of nonuniform condvetivity and lateral flow i is due to the 


ay 


| 
{ | 
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LATERAL FLOW 


(Refs. 2-7). However, instead of the numerical and thereby 
, ue the insight, we prefer the simple analytic derivation that provides a — 


to the rae s (Fig. 2); m and s form an angle a with the z and x coordinates, 

respectively. The flux qo is taken to be known somewhere in the medium. : 

Of special inte interest is 4 which is vertical in the evaporation) 
orin the negative z ‘direction (as infiltration). 


2". % bebaves an 


bee | ou 


a cyclic me dium and ju steady Naw, 
.—Upright and Sloping Coordinate Systems and R Related Fic Flow 


_ In order to solve an actual problem, we require information on the hydraulic — 


_ properties of the medium and certain boundary conditions. The most general _ 
case is of very little interest if its solution is too complicated. Here we we 
| interested mainly in two cases: (1) ' With uniformity along s, e. g., aW/as = 
 0,aK/aS = 0, etc.; and (2) with a uniform as in item | and a steady- 
While these are not the most general one, “they are sufficient to lead us > 
to some of the more important conclusions or at least provide us with a ‘good < 
‘starting point, lending an insight into hydrologic processes. ba fot wl? 77 


=> @ 

| 

bes nent 

Z 

According to Darcy’s law and assumption about unifor 


in the n, s dir n, s directions are beg 
(18) 
q, = —K\ — + cosa 
in which K(n, Vv) is the hydraulic conductivity that can vary ry with 1 the Pressure 
head, W, and the location along n, explicitly, but not inthe s direction. __ 
_ The fluxes expressed in terms of the x, z z system are derived | from Eqs. | 
18 and 19 by simple geometric transformation tas 


The ‘medium is assumed to be : isotropic | at every point. Therefore, qs and 7 
q,can also be obtained directly from the Darcy equation (the n and s directions 
are the principal axes and do not raise any am of isotropy): 1) ie 


= 

@q 
q z 


‘- 22 and 23 are in good accord with Eqs. 20 and 21 if we ononarll 
the gradient vectors in the uniform slope condition from the n, s system to- 
At every point it is assumed that the force and flux are parallel, so that 
the ratio of the flux components is exactly equal to the ratio of the force — 
components. ‘(simply divide Eq. 18 by Eq. 19 using the same K in both or 
divide Eq. 20 by Eq. 21 or Eq. 22 by Eq. 23. This is the essence of assuming» 
isotropy ata point, 
% The first question to be answered is whether or not the medium as a whole 
- behaves isotropically or not. In order to reach a general conclusion, let us 
consider path lines (that are identical to the streamlines only under steady flow). | 
2. ratios of path components in the s and n directions are found by dividing — 


Eq. 18 by Eq. 19, and integrating over n and dividing by § dn for averaging: — — s 


Note t that (q,/4,) = (ds/dn), along a path line so that the numerator in Eq. . 


24 is simply | ds. The entire fraction expresses the ratio of parallel to normal 
paths of flow. The fraction (q,/q,,) is the angle of the flux vector at a a point, 
while the r.h.s of 24 the of the flux vecto or. 


g 
j 
| 
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‘The rhs. is s obtained by finding expressions for q, and qn Eqs. 


-< 19. It is clear that if aV/dn varies with n, the overall slope of the water 
path also varies with n. The overall ratio of the parallel force F. to the normal 
force , is found by integrating the ‘parallel gradient “components se separately 
over n and dividing them by the normal component integrated overn: 
_ Clearly the ratios in Eqs. 24 and 25 become identical only in the following 
two cases: If aW/an = 0; or if sin a = 0, i.e., horizontal soil | layering. It 
is evident that if the 1 ratios of the overall flux components differ from those j 
: 4 of the overall force components, the soil generally behaves as an | anisotropic 
_ The ratio between Eq. 24andEq.25is al 


naa 


is smaller than or equal ‘to the arithmetic mean in the first integral. In effect, 
= means that on the average, the co conductivity i in direction s is greater than > 
c that in direction n. The equality to unity in Eq. 26 is obtained only when the - 
- numerator equals (§dn)*. This occurs only if a¥/an = 0, i.e., in a uniform 
« with a uniform flow regime. It is interesting to note that such anisotropy 


‘in the p preceding i isa a special case. 
G In a cyclic medium and in steady flow, ¥ fluctuates around an average value 4 
‘ a and aW/an fluctuates around zero. Thus, over a complete cycle, the first integral 


Eq. 26 26 Teduces to COS a and nd Eq. 26 rec reduces 


lig Care must be taken in the case of a horizontal : soil s since Egs. 26 and 27 
_ were obtained by division of Eqs. 24 and 25, both of which vanish under such 


In summary, in cases where aV/an #0 and a # 0, ami flux 


= from the force. If the force is on the average vertical, like gravity, there 


_is always a net horizontal component downstream. 
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_ We have demonstrated that streamlines curve downstream on passing from 


a less permeable layer to a more permeable one. Alternatively, within a uniform 

‘soil, the streamlines bend downstream w when the pressure increases with depth. z 

We have also proved that in any nonuniform sloping soil there will be a net 

horizontal flow component. Furthermore, we have furnished evidence supporting — 

our stipulation (15) that under conditions of vertical infiltration the horizontal — Fr 
flux is proportional to the vertical flux and to the slope. The total horizontal 
flow is therefore also dependent on the overall thickness of of the anisotropic 
medium in which the vertical flow takes place. M298 ui oats 
_ The simplest case is that of a saturated, layered soil. We have also dealt 
with another simple case, i.e., a uniform unsaturated soil with a water table. 
Hydraulic conductivity increases equeees the water table and there will be a 
lateral, downstream flow component. 

_ In otherwise uniform soil, another change in K can come about at the wetting | 
front. There, the pr pressure head and K drop to almost zero with depth. Therefore, _ 
streamlines that enter the soil vertically by infiltration will tend towards a direction _ 

- normal to the wetting front. If the front is at a slope, there will be a net 
horizontal component upstream, e.g., if aW/an is very high in Eq. 24 as at 

Unsaturated flow in a layered soil can show irregular behavior. In an unsaturated — 

) layer overlying a less permeable layer, the pressure head and K will increase 
with depth. The opposite will occur in the underlying layer. Thus the streamlines _ 

_ may describe a zigzag course moving downstream within the more Siimiidite 
ye and upstream within the less permeable one. However, the net effect 

_ is movement downstream, as has been proved in Eqs. 26 or 27. Livabteragiinoil 7 

During evaporation of moisture from the soil there will also be a a gradual 

_ increase in the pressure head with depth and the net effect in both uniform 

and layered soils will be of a downstream, horizontal flow component. 
‘In the following we shall consider qualitatively a two-layer problem and a a 
surface transition problem. In forthcoming parts of this report they will be 

i Let us consider a a permeable layer of saturated conductivity K ,, ,, and thickness. 

) (Fig. 3). The rain is of intensity 
oF which is ; smaller | than the hydraulic conductivity, | of the top layer. 

- Thus the flow in this layer will be unsaturated and under negative pressure. 

_ If D, is sufficiently large, the flow regime at the top will asymptotically approach , 
oe q, and the hydraulic gradient will approach unity (see Eq. 14 with z 

0). As the flow approaches the interface between the layers, the pressure 


increase se towards saturation, vIn Fig. two © adjacent vertical 
_ sections are shown. Compare the pressure head curve in the upper one with 
that in the lower one. At two points along a horizontal line C-D the elevation 
is the same. The pressure head Y is higher in the upper section (curve 1) 
than in the lower one one (curve 2). _ Therefore th there must be a flow compere 


j 
q 
| 
| 


in the horizontal direction, i.e., 
Let us study this problem in the light of the streamline equations 3 oud a 
4. At the top of layer 1 (Fig. 2), K, > q, < K,,. At the bottom of this 4 = 
b: layer K, = K’,,i.e., it increases and may approach K., sai» Owing to the presence = 
Of the less permeable layer and as a - result of the unsaturated flow, there is a 
a build-up of moisture, pressure, and conductivity above the interface. The — 


flow direction in the top layer will then change from vertical(g,=0)to 


thus evident 
~ that if ‘there is a "more ‘permeable layer at ‘the soil surface and if the flow 
is unsaturated, the streamlines will bend downstream. It can be ned that 


par 


2 


g + every soil, without exception, has a more permeable layer at its surface. Therefore, 
. : in every sloping soil subjected to prolonged rainfall, the streamlines will bend 


downstream. In other words, at the surface of every sloping soil there will 
be a horizontal flow component downstream. This is true both under completely 
| unsaturated conditions and in the presence of a water table. Furthermore, one 


1. The horizontal flow component is proportional to the vertical flow at all — 
depths. The same also holds for the total horizontal discharge. Implicitly, the _ 
horizontal flux also increases due to the coefficient of anisotropy ((K"/q,) 
= 1] which i in ) turn increases with the rate of vertical flow. ‘Thus, , the horizontal 


= 2. If there is a change in | slope, the concave, there 


4 
q 
~ 

né 
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is also a concentration of moisture. This is han to the fact that horizontal — 
inflow exceeds the horizontal outflow. A very high saturated hydraulic conducti- — 

_ vity (K,,) usually implies a rapid reduction of conductivity due to suction, — 

i.e., a higher value of the coefficient a in Eq. 12. This means that a relatively — 

small soil depth (z — z, in Eq. 14) is sufficient to make the conductivity equal — 
g the rate of infiltration and to cause the flow to become vertically downward. 

_ Thus the very high conductivity at every soil surface validates our ee 
that the rain initially penetrates the soil in a | vertical direction. The two- layer 
problem will be treated quantitatively in a later pa part of this report. Te > 


bf. In the first part of this report (15) we expressed reservations 7 ahett the possibility — 4 
of measuring surface flow because there is no ) explicit definition of what actually P 
constitutes the soil surface. In Ref. 14 it is argued more generally that in soil 
physics meaningful entities are averages over time and space and produce .* 
continuum. The same reasoning should be applied to the soil surface, e.g., _ 

_ porosity can be measured over a finite sampling volume or over an a 

a the air well above the ground it will be 100%. Somewhere in the soil it may 

be 50%. At any intermediate point it changes s gradually. Hydraulic conductivity — 

4 changes in a similar way from some finite value within the soil bulk to a very 
_ high value in the air, e.g., the change in K,, in Eq. 12. In a similar way the ; 
air entry value \, will be reduced to zero at. the transition from the soil bulk 
to the air (/, — 0 in Eq. 12). Finally, the rate of K reduction under suction 
will increase [coefficient a in Eq. 12 (See also Ref. 10)]. The air can therefore 
be considered as some limiting form of the porous medium itself (very simil 

to a very coarse gravel). The pressure of the raindrops is always atmospheric. | 
The rate of flow well above the surface is in fact the rate of rainfall, and 

_ the unsaturated conductivity is conveniently equal to this rate. The surface 

; transition layer is not merely a mathematical artifact. It may be very thin in 
- some uniform, smooth sands, but can obtain a thickness of several decimeters — 


_ The soil surface is defined by a transition of properties. Its direction is defined _ 
6 a surface normal to the property gradient. One of the many such amrnees 
could be chosen in accordance with some convenient conservation requirement, 
eg. g., that the total porosity -Temains unchanged (14). The said property may > 
change abruptly at this representative surface. We may conclude that a clear-cut 
“soil surface can replace the transition only as a mathematical artifact. 
F. important consequence of the aforementioned is that every rainfall, 
_ irrespective of its intensity, penetrates completely into the soil (at least into 
its surface transition zone). If we are consistent in our views, » there is no si such © 
thing as surface runoff over a . large surface because the rain can never exceed 
the hydraulic conductivity of the uppermost part of the soil surface, i.e., the | 
air itself. Saturation due to rain of high intensity, when such occurs, will always 
be manifest initially within the surface transition zone. 
‘The splashing of raindrops has been proved (16) to cause a horizontal a 


component proportional to the slope, explicitly proportional to the rate of rainfall 
and implicitly proportional to some fractional power of the rain. Evidently, — 


s steady flow — the surface transition layer is related in aap the same 


LATERAL 


g way to the rain. The 1 response. time of ‘the : splashing raindrops is measured 


_ in fractions of a second. It can therefore be considered quasi-steady (behaving _ 
 Meatically to the rain reas Flowi in the transition layer can be e delayed, depending 


_ which is expected to be normal to the surface. Therefore, the flow will tend 
‘to be normal to the soil surface (¢, — 0) (Fig. 1) and even to have a slight 
upstream flow component (q, < 0) (Fig. 1). Well behind the wetting front the | 
main force will be gravity and a horizontal downstream flow (q, > 0) (Fig. — 
- The flow in the surface transition layer can sometimes be observed as s tiny 
trickles of water or shiny soil surfaces. Very small concave spots where —_ 
concentrates can be considered in detail as such or be averaged out as parts 
_ of a thicker transition layer. The exact limit depends on the arbitrarily selected 
scale of observation. Calculation of a steady flow through a transition layer — ei 
. possible if the hydraulic parameters are known. Nonsteady flow will be treated me: 
in detail in the last part of this report. 
SiGNiFicANce oF Frows in Surrace Transition LaveR od? 
a us first consider a well-cultivated, heavy soil, with a good stable structure. , 
‘The hydraulic conductivity at the top (not necessarily near the air) can reach 
‘ 10-' cm/s or even | cm/s. Such values have actually been measured by the | 
senior writer in a drainage research field at Hazorea, Israel. In the subsoil ¥ : 
i the hydraulic conductivity may drop to 10°° cm/s or less. This means =“ 
the anisotropy coefficient (K,/K,) -— 1 can range from almost zero at an 
extremely low rate of rainfall up to the order of 10°. One can actually 
serve lateral flows of a few meters or tens of meters over a vertical infiltration 
_ of a few decimeters. After rain or irrigation water puddles form in slightly — 
_ concave spots. This example is obviously extreme and almost trivial. The 
observations reported earlier north of Beer-Sheva (15) are far less trivial and © 
‘fit the foregoing analysis. Even a change of two orders of magnitude in 1 the 
ia hydraulic conductivity within the surface transition layer and a slope 


can produce a horizontal flow equal to” the vertical one. The — in 


We have no direct measurements of field anisotropy. The of such 
per and corm of certain other aspects s should be the subject of = 

3 The two processes, raindrop splashing and flow in the surface transition layer, 
can combine to produce lateral flow and moisture concentration. Under laboratory — 
conditions one can test the flow within the transition layer separately by applying — - _ 
moisture lacking the high kinetic energy of the raindrops. In nature it will be 
difficult to distinguish between the two. In an earlier part of this report (16) = ff 
we mentioned an experiment involving a a V-shaped soil slope. Runoff started 
F almost immediately at the sharp bottom e edge. Raindrop splashing alone cannot — 

account for this. Since in this case there was no layer or pervasive 

_ saturation of the soil directly by the average rain, the phenomenon may be 
related to flow in the surface transition — ‘This deduction is also supported 


— 


by the fact that, the runoff a at ‘the ne 


~~ Apparently, the water flow in the surface transition zone can han quite ont a 
in cases where raindrop splashing is unpronounced. It can respond to more 
: extreme curvatures of the soil surface (greater than | m~') where it can produce 
runoff and erosion with very small amounts of rain and within a relatively 
short time. At the other extreme, it can lead to the accumulation of moisture _ 
at relatively moderate curvatures (smaller than 0.1 m"'), as in the Beer-Sheva 
experiment where the effect of raindrop splashing was negligible. 
The thickness of the transition zone is dependent not only on geometry s 
also on effects such as the transfer of momentum in liquid flow (12). Most — p 
soils have a looser structure at their surface with some aggregates, root holes, - 
and other irregularities. Freshly exposed soil cuts will develop such pool 
transition zones in the course of time. Surface lateral flow and erosive mechanisms 
_ Some surface transition layers may be very similar to a straw roof. Examples — 
are litter-covered forest soils and possibly even certain grass-covered areas where 
the old growth may have a marked orientation parallel to the soil surface. __ : 
The concept of surface transition may find application it in the solution 


face. Accordingly, the “flow medium will be « as consisting of “highly 
permeable layers at the surface. The concept of a seepage face thus becomes — 
redundant. The streamlines simply bend downstream in the transition zone, 
diverging from the normal on entering more permeable layers. In what has 
‘ been termed the seepage face, the flow is more or less parallel to the surfac 7 
but occurs within the soil. Such a comprehension of seepage is significant not 
only in contributing to the physical consistency of our flow analysis. It can — 
explain how certain small, geometrical, surface irregularities can cause local 
outflows and erosion due to seepage forces. 


__ Itis significant that a downstream horizontal flow component occurs whenever _ 
= pressure head and K decrease with elevation. This is certainly the case 
during drainage and the drying of the soil surface by evaporation. The presence = 

| of a surface transition layer heightens this effect. Thus lateral flow and moisture 

concentra s of the landscape will continue long after i 

‘fain has stopped. This process will keep the concave parts wetter for long — 
periods and shorten the time interval preceding the onset of runoff when the | 
next rain falls. It appears that lateral flow takes on universal importance during © 
saturated and unsaturated flow, prolonged rainfall, and the drainage and frying — 
of the surface. It occurs at seepage faces, in natural, exposed soil, and possibly 
4 The derivations were made for a soil which is isotropic at any point. The 7 
nonisotropic behavior is only on the average. This derivation does not exclude _ 
= where the soil is intrinsically eee aaa The mathematics is only more 


and the phenomena more extreme. well vet 


BE _ 
’ edge involved liquefaction of the a 
4s ana © OU a can occur only if water seeps out of the 4a 
7 soil. Raindrop splashes may be considered to remain outside the soil. However, 
| 
| 
4 
| 
= 


pore 


- Famer have been emphasizing the importance of leveling fields. The term a 


% leveling appears to be a semantic error bound up with a misinterpretation of is 
the mechanism. In practice, the leveling of fields entails the provision of plane — . 
_ Surfaces (also having gentler slopes). It is possible that leveling is not as important : 
_ as “‘planing.’’ Very small scale, sharp, local curvatures, such as furrows iokh 
- peanaifornisies in the moisture distribution that are averaged out by the soil 
_ itself and by the plants. Nevertheless, sharp-edged furrows invariably lead to” 
erosion and the rapid development of runoff. In the light of the previous analysis, — 
- some thought should be ‘given to the shape of furrows. Nonuniformities in 5S 
_ the moisture distribution can also be caused by more moderate curvatures. ; 
_ Moreover, they can occur on : such a ‘scale that cannot be evened out by the , 
root volume of a single plant or by lateral redistribution of water at greater 
_ soil depth. Thus a plane surface configuration of a field is a prerequisite for 


even distribution of rain or irrigation water, 
Fae We have completed the description of the four mechanisms leading to lateral | 
| flow components: (1) Splashing of raindrops; (2) flow within a surface transition — 
a layer; | (3) flow in a layered soil; and (4) bending of streamlines near a water 
table. The last three are merely different forms of soil nonuniformity. However, - 
they are here arranged according to depth andtime lage 
In all four phenomena the horizontal flow is found to be proportional od 
the vertical flow. All four have been shown to possess possible significance. — 
“ They could thus explain m most of the phenomena m mentioned in the first part 
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conductivity; 
_ water flux i in velocity dimensions; 
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as = slope of soil surface or soil layers; 
yz angie between incoming streamline and normal to interf: 
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= pressure head; and 
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Hypro.ocy: IvV— —Fiow IN SLOPING, ‘a 


as one » of the causes of a lateral flow ‘component induced by infiltration. It 
has been generally proved that as soon as the soil is not uniform and it is 
sloping, there must be a downstream flow component parallel to the slope. 
In the following we shall treat specifically a soil consisting of distinct layers. 
_ If the hydraulic characteristics of each layer are known, it is possible to determine — 
the flow and calculate the total horizontal - discharge. In fact, the description = 
of the soil as being made up of distinct layers also enables us to — 
pecific case of a cyclic two-layered soil and steady-state flow, the 
soil actually behaves as an anisotropic medium. The average parallel and normal 
_ conductivities can be calculated. The following is based to a large + on 


_ Consider a nonuniform soil with properties changing in | a direction n. T 
parallel direction is s (Fig. 1); n and s form the 


4 


We assume first: (1) Uniformity along s S, €.g., avsas = = 0, aK/aS = 0, 
_etc., VW being the pressure head and K the hydraulic conductivity; (2)a uniform. 
“dese as just mentioned and a steady-state flow; and (3) the medium is cae im 


at any point. Any anisotropy is, on the average, induced. = 
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; a ___In the three previous parts of this report soil laverine has been mentioned _ 
| 


ccording to s law Fig. 1), the faxes in in the n and s directio 


<4 
on 


in which K(n,V) = baad hydraulic conductivity that can vary explicitly with 


at 


at | lator od? bar wath 
FIG. 1.—Upright and Sloping Coordinate eeume and Related Flux Components 
The ‘flexes expressed in terms of the x, system are. are obtained from 


"ole to point isotropy, g, and q, can also be derived | directly and simply — 
g the Darcy equation (the m and s directions are the principal axes and 


qx 


to Eqs. 3 4 if. we 


Sare: 
| 
| 
- 
| 
| 
radient ve tors 


in Eqs. 12 and 13. 


control surface of Fig. 2. 
ave Rave to" Sos 


—P cos cos a dA= = 


P being “the intensity of horizontal above 


Ss 3 
These three equations (two in Eq. 8 


_ According to Eqs. 8, 9, 3, and 4, ‘the flux components in the x wy site — 


) = -(P cos” a + K sin? a) 


- si ina 


A from the n, s system to the x i 
accordance with the c 
same fluxes cross the soil surface and any other surface parallel to it. On a 
ovide the basis 
| 
a 
— av 
| 
= 


—=s—tana = 


in which z = the elevation of some surface n = constant 
By substituting the left-hand side of into inte, 10 we ot obtain yet anothe 


can be derived from 11 in one of two forms 


“Ths the case, Eq. 15 becomes: 
=(q, + K) tana =(K P) cos? 
Note that during a rainfall gq. has a negative value directed downwards. If — 
the pressure decreases with depth, the gradient i is larger than when q, is — 
upstream, and if it increases with depth, q,is positive. 
3, 10, 


92 

or Q,=$q, = tan a cos” K dz — P 

Bei convenience we define the a averages: a 

=—Sq.dz; . 


7 


tana = q, U tana D; cu 

in which G, = the average downward vertical flow. ‘From Eq. 19 one obtains 


an alternative expression: for the ne average horizontal flow with steady rainfall — 


,D= pus tan @; 
. Eq. 214 gives ‘the basic ratio between the hostnentel and vertical components 


of flow, which applies both in steady and nonsteady flow. Eq. 22 provides — 
a a similar ratio between the rain and the average horizontal flux, which applies — 
i in the case of of steady fl flow. Note the magnitude of is ‘not ot necessarily 


56 
q 
4 
j 

| 


is on ‘the x- or White angles a 
are sufficiently small, cos’ a — 1 and Eq. 22 is simplified. Anal ore 
The horizontal flow component is the value we actually need. It is related 
z= 4s to a unique coordinate system and can be integrated over a map. We shall — 
therefore have to continue using it although the employment of n, s coordinates 
_ The total discharge in int aad s direction i is Q,. By integrating Eq. 9 we obtain: a 


both in n steady and flow. a? 
_ The average parallel flux component g, is 


and 22 again feature the basic that in 
® first part of this report (1) and was later proved for splashing raindrops 
2). This relationship has also been demonstrated for the bending of streamlines _ 
$e far the uniformity of the slope has been limited to the requirement that | 
“aW/as= * 0 or that this term is negligible in the calculation of the first approximation : 


of the vertical ‘flow regime. Eqs. 21 and 22 are not limited to a two- -layered Do 
- eyclic soil, but apply to any soil with a varying conductivity along the n-axis, 
whether due to soil properties or to the boundary conditions. Eq. 21 is valid i 
for both nonsteady and steady flows, while Eq. 22 is only valid for the latter. 
In specific cases we have | to assume relations between K and Vv and find nd the 


AN ano Orvers of Macnituve 
Net Recharge of Deep Groundwater. —tThe net recharge is introduced into the 


_ sell in the form of seasonal pulses. ‘However, these seasonal pulses are damped 
and the flow. becomes virtually steady after several wavelengths. In the thane 
a equations P must now be considered not as a rain but as the net groundwater — 
recharge. If the net recharge is P as and the moisture content is C, then 
4 the wavelength is about P/C, e.g., if the moisture content is 0.3 and the a 
net recharge of 
300 mm/yr is equivalent to about 10" cm/s. 
‘Let us now consider a series of soil layers” oe from cy or sock with 
a conductivity approximating the rate of net recharge, i. €., 10™ 
permeable ones in which moisture accumulates. = 
The average ratio between components parallel and normal to the layers is 7 : 
given by Eq. 25. The parameter K, will be affected by the layers with high 
conductivity. _ The highest conductivity values will be obtained for coarse-grained — 
‘ soil i in the presence of perched water. Clearly the ratio can attain several orders 


= 
‘ 
1 
a 
: dy-state flow we obtain from Eq. 8 the very simple expression , 
| 
[ 
\ 
| 
i > 


of magnitude. The same can be estimated from Eqs. 21 or 22. The coefficient — 
of anisotropy is obviously not a constant. It varies from zero at low flow rates 
Flow Behind Wetting Front.—Sufficiently far behind the wetting front - 
flow can be almost steady. Such flow is established within a relatively short _ 
time, especially if there is an impeding layer at the soil surface. In such a 


Two-Laver Propiem 
_ Consider two parallel soil layers 
Their thicknesses are D, and D,, respectively, and their saturated conductivities _ 


f 


ap 


FIG. 3.—Two-Layer Pressure Head Distribution: (a) State A—No Saturation, State 
- B—tnitial Saturation; (b) State C—Partial Saturation; (c) State D—Complete Satura- _ 


. 7 case the formulas incorporating rain P can be used as a cineca 
| 


K, and K waz» Tespectively . Each layer is isotropic and uniform. The 
consis of a cyclic arrangement of these two layers repeating itself many times. mes. 


“cos dn: = 
& fas 11 we similarly obtain along the z coordinate 


‘The boundaries of integration are shown in Fig. 3 which palit two layers 


« four flow stages. The conditions are presented in Tablel1, 


The actual values of ¥, and WV, can from an equation 
as a function of the rain P. 


" TABLE 1 —Conditions of afta 


Stage of flo flow n- ‘boundery 


unsaturated 
appearsin 


Partial saturation so 
‘ -n, is saturate 7 
completely saturated 1 —(D, + D, ) v 


The solution is readily obtained by introducing a relation K() for the two ; 


layers. Then WV(n) or ¥(z) is found as well as K(n) or K(z). From these, 
turn, one can calculate 0, from Eq. 19 or from Eq. 23. 


As an example, let us consider the experimental relation for the ith layers _ 
with K,, as the saturated conductivity, and a, and WV, as the acunadl 


K,= Ki exp a,(W + Yan for V+ 


cos a(n — n,) = — 


a 
| 
J 
e| n-boundary W-value Figure 
4 

| 
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ee 29, when wt written for two layers of thickness, D, and D,, gives ns ieied 
exp (— D, cos a) = 
We cannot | express VW, and W, explicitly in an “piper manner, although 
these are the two unknowns in Eqs. 
Note that the boundary conditions in the foregoing ta table assume transition — 
to at In this transaction occurs when 


Pressune DistRIBUTION STREAMLINES AND Averace 
the case of a cyclic repetition of the t layers, certain general 
deductions can be made. > 
a us consider Eq. 26 and calculate K, _ according to its definition in Eq. 
24. Eq. 29 is most convenient for this purpose. If we — K, as a variable 


~ 


[exp(-a,cose wt 


[exp (—a, cosa D, 
a,cosa 
_ We now prefer to have K . and K adjacent to each other on either side — 


of a single interface betwen the two The condition at interface 
is that V is the same on both sides, e.g., V,. 
_ Consider first a case in which = K,.. this | the 
q K(¥) curves of the separate layers intersect each other. In this special case 


the layers seem to ene is to zero. Then: 


oa) 


on this value, - the layer which has the higher saturated conductivity will 
more permeable. If it is less than P,, the layer which has the lower saturated 
cnithestiiie will be the more permeable one. The transition value of rain P 
is very small. On the whole, the phenomenon it 


al 
] | 
: with the conditions of Eq. 33, all the terms in the brackets vanish, the integral j 


: ‘The streamlines zigzag from layer to layer. In the more permeable layer the 
_ K build-up with depth and the streamlines diverge from the normal in a downstream se j 
direction. In the less permeable layer the streamlines converge toward the normal. 
At the interface the streamlines undergo abrupt changes in direction. — a 
Ata gh rate of rainfall when K,, = K,,, and K,, = K,.2; Eq. 32 becomes 


K 


an s indeterminate term. For this « case, and for decreasing D, and D, values, _ 
it is suggested that the exponentials be expanded by a series as follows: w ) 


exp (na, cos a D,) — cos - 


Gn tana 


‘ 


coefficiemt 


“Cha © 
FIG. 4, —Anisotropy tropy Coefficient Uas Function of Rain Intensity. Pin Cyclic Two- Layer 


P 
ofa,D,ora,D,cosa,weobtain 
| 
| 
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mean of K. It is the exact mean in pronto flow and the approximate one 
r unsaturated flow when the layers’ thicknesses D, and D, approach zero. 
_ The coefficient of anisotropy | ( K./G q.) — 1 has been calculated as an example 
with Y, = 0. The abcissa in Fig. 4 gives the rate of rainfall P in terms relative : 
to K,,,. Two specific f points are given: P/ et = 2 28 where saturation appear. ‘ 
at the surface; and P/ K,,, = 0.63 where P = P, and there is no lateral flow. _ 
* The method used for a two-layer problem can also be employed for a transition 
zone at the soil surface if the latter is represented by a sequence of layers, 
each uniform and isotropic and varying slightly from its neighbors. Some of 
the details of such a study are interesting, e.g., if the s saturated a al 
conductivity decreases monotonically with depth, there must be a maximum | 
in the pressure head distribution within the surface transition zone and above | 
_ this maximum the streamlines bend downstream; below it they bend back towards a 
the normal. As the rain gradually increases, saturation will appear first within 
te surface transition zone at this point of maximum pressure. The site of 


ye upward as the rain intensity increases. 
Wis is “instructive that the coefficient of anisotropy increases with» the rate 
_of rain to a power higher than unity (at least in the two-layer problem). Therefore, A 


“the lateral — increases with the rain intensity to a power higher thantwo. 
“= 


Water ConcentraTion AND Erriux 


Lateral flow due to ‘splashing | of raindrops and to thin su: surface transition layers 
- follows the rain after a short time lag and may be considered to be quasisteady. 
7 Its rate was shown to be proportional to the soil surface slope. In the panes. 
-~part of the report any layering is considered which does not have to be parallel | 
to the ‘soil ‘surface. Therefore, the moisture will not ‘necessarily concen 
only in concave parts of the landscape but also in concave parts of the lay § 
~- There are many cases where the soil layers are more or less parallel to the 
soil surface. Sand dunes serve as an example. On observation one phot 
. such repeated layering, sometimes a few millimeters thick, due to segregation — 
_ Of the particles by size and slight chemical and physical stabilization of freshly — 7 
sand layers. ‘Where the soil layers are parallel to the surface, the 


moisture concentrations due to splashes, the surface transition zone, , and the 
“a - Moisture concentration is also due to other causes. The most common one a 
‘ is the truncation layers due to an excavation or an erosion cut. In an unsaturated — 
soil the diagonal streamlines approach the soil surface but cannot emerge. They 
then bend downward and moisture accumulates until saturation is reached and 
seepage from the soil commences. Geological faults have been known to give 
rise to the outflow of water as springs. This was hard to explain | when no 
saturated water table could be observed in the vicinity of such faults. _ The 
aforementioned unsaturated lateral flow offers an 


4 In pollution studies one wishes to follow the actual flow y path of the water. 


| 

| 


— Itcan can be prec predicted that infiltrating p< polluting water can move laterally over large 
- distances before reaching the groundwater table. At one observation site, more | 
oor, m of interbedded clay lenses and sands were sampled above the _ 
groundwater table. It was estimated that the spreading diameter of a point source 
- of polluting water would continue to increase as long as some impounded water 
- _ remained within the sandy layers. In this case the coefficient of anisotropy - 
| is of the order of the ratio between the conductivity of the sand and that 
' of the clay conductivity. The slope was recorded as a percentage andthe anisotropy _ 
_ in thousands. Thus, the pollutant was expected to travel a horizontal distance y 
of several hundred meters to several kilometers before reaching the water table. 
_ Observations tended to confirm the qualitative predictions. Surface pollutants — 
may readily move near the surface in depressions and rills and gullies where cc 
the onset of salinization is ‘early and water seeps out rather than in. ip lowueel 
7 _ The following may be considered an interesting application of lateral flow _ 
in regard to pollution. A convex cover of soil layers, if properly constructed, ig 
can act as a “‘straw roof.’’ Though highly permeable, it can prevent the péuitisting 
_ Water from leaching through a source of pollution. Very often engineers restrict 
a - their search to impermeable materials when attempting to prevent flow through - 


a structural element. Here, we suggest another possible way of obtaining such > 


_ Layered soil usually behaves as an anisotropic ‘medium. If the layers are 

at a slope, the force of gravity operates at an angle a to the principal coordinates _ 
‘ n and s. The coefficient of anisotropy relates the parallel component to the — 
7 normal component by sin a, and the parallel component to _the rain by tan 
q a. The average horizontal hoe component is related to the rain at steady flow 

by ((K,/P) - - lj cos* a a, , and to the vertical flow component by (K,/G, )- 

. The last parameter increases with the rain to a power greater than unity. 
The relationship between the horizontal and vertical flow components is remini- 
scent of the rational formula for runoff. However, runoff is generally replaced — 
by horizontal flow, and the coefficient increases with the rainitself.§ = 
_ The horizontal flow leads to moisture concentration. However, the latter does 
“not necessarily bring about saturation and seepage f from the soil. These depena 
on the thickness of the profile involved. 

An interesting outcome of this analysis" is that j it provides. a hint about the 

- possible form of the conductivity tensor in unsaturated flow. =~ 
In certain special cases the we sagea,. 4 can be calculated from the properties — 
of the individual layers. However, in the future it will be necessary to adopt 
‘posi method od by which horizontal flow can be estimated from the anisotropy. 
The latter may may be obtained from watershed measurements using some overall 
-Telationship between the horizontal and vertical flow components. 
Also of fundamental interest is the fact that the conductivity pressure head — a 
curves of two different soils intersect at some point. This crossing has an effect 
on the overall behavior of the ee soil. Some models that assume parallel 


= applications relate to problems of pollution. The pollutant may “ 
spread over a significant distance horizontally while into the ground. 
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| 
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interest in solving engineering problems. 
_ The next part of this report will deal with the computation of nonsteady 
Zaslavsky, D D., and Sinai, G., ‘Surface Hydrology: 1—Explanation of Phenomena,’ 


a Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY1, Proc. Paper 15958, 
Zaslavsky, D., and Sinai, G., ‘‘Surface Hydrology: II— Distribution of Raindrops,” 


Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY1, Proc. Paper 15959, 


Zaslavsky, D., and Sinai, G., “Surface Hydrology: of Lateral Flow,” 
_ Journal of the Hydraulics Division, ASCE, Vol. 107, No. HY1, Proc. Paper 15960, 
4. Zaslavsky, D., and Sinai, G., “‘Lateral Flow Due to Rain in Unsaturated Soils,’ Second - 
A nnual Report, Project A 10-SWC-77, Agricultural Engineer- 


aa sy mbols are used in this paper: cboon 


= thickness of individual layer; 
P = sate of rain («LT ) or rate of net recharge; 
q = flux vector or flux compensiat 
= = coc coordinate parallel to Ba ft 
= angle between soil layers and horizon; and 
v= head for ¥ > 0, suction for < 


"relating to interface between 


= relating to point where unsaturated conductiviti 
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IntropuCTION H t 
_ The final | part of the report on Surface Hydrology deals with a more problematic 7 


aspect, i.e., the numerical computation of the nonsteady flow through the surface 

_ transition layer. We have previously established that this transition layer has 
iy properties grading from those in the soil bulk to those in the air. ‘Sometimes _ 

| : can be fairly thick, as in a plough layer. In some well cultivated, heavy, od 
clay soils, the hydraulic conductivity can m change by six orders of 


layer i is greatly reduced, ‘its thickness being of of the ; grain diameter. 

: . should be noted that a change in porosity alone, e.g., from 30%-100%, 

can alter the hydraulic conductivity by one order of magnitude or more. 

i _ In an earlier part of this report (8) we presented a general proof that whenever > 
_ there exists a sloping transition layer the streamlines of the rain will enter ‘ 
the soil a and then bend downstream. During steady r rainfall the horizontal flow . 

will be proportional to the rain and the slope. Does | this a also apply in the case 
_ To date, the hydraulic properties of the transition layés have never been 
actually measured. We regard it as being made up of a discrete number of 
_ layers changing stepwise. Some of the parameters assigned to these layers are ; 
taken from actual soils; others are our proposals based on certain rules. tod 
_ are still at the stage where we required knowledge of the nature of flow phenomena _ od 
in principle. The exact numerical result is of minor significance. The numerical © 
: solutions have been used as a method of experimentation to study the influence | 


of the slope, the curvature, the rain regime, and the thickness of the transition be A 


‘Prof. of Soil and Water Engrg., Faculty of Agricultural Engrg., Soils and Fertilizers 
Technion—Israel Inst. of Tech., Technion City, Haifa, 32000, Israel. 
“Leet., Faculty of Agricultural Engrg., Soils and Fertilizers Div., Technion—Israel 
Inst. of Tech., Technion City, Haifa, 32000, Israel. 
Note. —Discussion open until June 1, 1981. Separate discussions should be submitted 
for the individual papers in this symposium. To extend the closing date one month, a 
x a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on March 23, 1979. 
_ This paper is part of the Journal of the Hydraulics Division, Proceedings of the American — 
Society of Civil Engineers, © ASCE, Vol. 107, No. HY1, January, 1981. ISSN oul 
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layer < on the flow regime. Some major results and conclusions from the numerical | 
inthe on are given here. Detailed solutions and some more results can be found 

in the work by (5). 


i, A mmneited soil surface is depicted as ‘in Fig. ‘1. The t transition layer is a 


underlain by a thick, uniform soil (H). _ The banat of the soil surface | is ” 


‘mee i n = ci sublayer depth which is 10 m in all cases. The amplitude A 
rw in the different cases from 2 x 10~* m-2 m. However, in most cases, a 
=2 m. These extreme values | were - chosen to ‘produce a given range - 


; pa The half-cycle L has the values 40 m, 20 m, 2 m, and 0.2 m (see 


_ The boundary conditions are as follows’ Fig. |): 1) The vi vertical line below 


“the hilltop is a symmetry line and thus a streamline; (2) the vertical line below 
4 the valley b bottom i is a symmetry line and thus a streamline; G3) the initial condition 


is a hydrostatic state throughout the profile, with the phreatic surface ¥ = 


basd asdi bev | 


sree as? ~~ 


a 
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FIG, 4 .—Geometry and Boundary Conditions of 
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It can an be _argued that a state of field capacity could | be , substituted for 
true state of However, the former depends on the history 


Main Geometric meters 
Data, in meters — 


“Intermittent rain. 


of its attainment. “This substitution should therefore be scans until a 
ime a as computers a are able to digest a long-term regime. 
_ The transition zone at the soil surface is represented by cight layers 5 = 


if 
1.—List of Solved Cases and Main Results 4 
1 Geometric in T. Time | 7. | Time 
Ratios” t | in. steps 
num- Li centi- in | steps, | centi- nm 
ber | L | hour |A/L “meters hours | TS,,, hours | 
388 2 | 2 0.45 | 10 |0.1 | 0.005 | 149 | 149 | 107 a3 | 21.3 
338 2 045 | 1.0 |0.1 | 0.005 | 13.8 | 138 85 | 19.1 | 19.1 | 102 
539 045 | 1.0 | 0.005 | 13.8 | 13.8 18.75 | 18.75 | 1 
— 0.45 | 0.01 | 0.1 | 0.005 | 15.5 | 1550 163 | sooo 155 
| 2 «| 045 | 40 | 0.1 | 0.005 | 17.2 58 | 60 7.3 09 
460 ~ 0.45 0.1 | 0.005 31.5 | ibood) 
| 2 | O15 | 10 | 0005 | 7.1 7.7 
752 F | 20. 7 0.15 | 0.01 | 0.1 | 0.005 8.0 | 800 1200 | 50 
536 20 0.05 | 1.0 | 0005 | 26 2 | 34 | 62 | 24 
33 |2 | 2 |o 1.0 |01 | 0005 | 097] 67 | 8 | 67 | 8 
535 L | 40 0.45 | 0.01 |0.3 | — | — | 
| 2 | 02 | 045 | 1.0 | 005 | 104 | 104 33 | 24.4 
s2iK| 2 | 02 | 045 | 0.01 | 005 
| 0.02 | 045 | 1.0 | 0.01} 0.005 | 16 | 6 | 6 | 38 
3791] 2 | 002 | 045 | oor ooo | — | — | 
378 A} 2 | 0.02 | 0.15 | 1.0 | 0.01 | 0.005 73 | 751 2B 14 
38H} 2 0.02 | 0.15 | 0.01 | 0.01 | 0.005 | — 
2 | 02 | 015 | 10 | 01 | 0.05 56 56] 84 
J 2 | 02 | O15 | | O05 | — | — mio 
0.05 | 1.0 | 0.1 0.05 | 24 44 ” 5 
«| 0.02 | 0.05 | 1.0 | O01] 0005 | 32] 6 | mw] 52) 8 | 
BY} 0.02 | 0.05 | 0.01 | 0.01 | 0.005 — — | 
$24 | 2 | 0.02 | 0.005] 05 | 0.01 | 0.005 | 2 82] 4 46 
$25 | 0.2} 0.02 | O15 | 1.0 | Os | 4s is | 4s 
525 P| 0.2] 0.02 O15 oor | os 
526 0.2 |2x 10] 0.15 | 1.0 | 10° |5 x 107 88 | 25 Fre = 
527 | 02] | 0.025] 1.0 25 | 42 | 42 2 4.6 
529 | 0.2 |2x 0.025 | 05 [10° {5x10} arg 4 
0.2 0.1 0.005} 05 |05 | 25 | 2 | 86] 42] 2 | 86 
0.2} 0.02 | 0.005 | 05 as | | ee] 
531 U| 0.2 | 0.02 | 0.005] 001 | 01 | J 
$32, | 02 |2x 0.005) 05 23 | 78] 46 | 46 
$32 W| 0.2 |2x 0.005] 0.01 | 10° |Sx 10°] — | — 
| 4 
. rid points representing the separate layers were chosen in such 


= 


| 


a that they comme uniform thickness of each in a direction normal 


the soil surface. Thus the successive sublayers do- not exactly | follow Eq. 
_ Each layer is considered to be uniform and isotropic. There is a “stepwise 
i ~ in the hydraulic properties of the various layers. The detailed values 


of the porosity, the retention curves, and the saturated a conductivities © 
presented i in 1 Figs. 3 and 4. 


1 


a” 


y the slope i in 1 Fig. 4, increases as as the saturated conductivity i increases. “ea a 7 

_ 3, The intersection of the conductivity curves for the layer pairs 1-2, 2- 3, : 

3-4, etc.,followsa monotonic order, 
4. The moisture content at saturation is the same as the porosity and therefore | a 

it epee 100% in the : air above the top — and it decreases with depth: 


5S. The rate of moisture decrease due to ‘suction, as expressed by | the 
in Fig. 3, increases with porosity. — 


- wt The intersection of the moisture curves for the layer pairs 1-2, 2-3, 344, 


7 
| 
| 
4 
| (5). ose for layers |, 2, and were selected artificially in conformance wit - 
| 


ee 
pot bur is rather sped as an 


which the ith flux. = the vobamettic moisture content; 


av Moisture Content, &, (Volumetric, cm3/cm3) 


; = FIG. 3. ee Retention Curves of Eight Sublayers Constituting Surface Transition 


‘ in which K,, = the hydraulic conductivity capable of being an anisotropic tensor; _ 

@ = the hydraulic head (6 = z + ); z = the elevation; and | = the | a 


a 

| PIN _ ore 
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= 
Inserting Eq. 3 in Eq. 2, one explicitly obtains samusneratl 

& q y 


= — | K,— + K W)K;, 1z- =0 .......--@ 
in which the hydraulic conductivity K,, = a product of two terms; KW) = 


Kio) 
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earlier dealt with anisotropic media but did not provide computations. In the 


5 
2 
= 
— 
git & 
| 
FIG. 4.—Hydraulic Conductivity as Function of Moisture Content K (6) in Eight 
Sublayers Constituting Surface Transition Layer = walks 
in reality, although the numerical algorithm used here is capable of solving 
flow problems for anisotropic media, the problems run here all involved isotropic 
In the problem under consideration, the anisotropic behavior is merely 
7 the result of the overall average behavior of the layered soil. Models proposed y 


following computation, ‘the anisotropy is not assumed to be a fixed value at »” 
a point but is rather obtained as an overall behavior. The anisotropy changes “4 
with the moisture regime itself. is contrary to the inference from Eq. 
4, where ‘the entire change with moisture is linked to the isotropic term K*(). 
Therefore we suspect that the assertion in Eq. 4 that was never of say measured — 


; the previous parts of this repor (6-9), leads us to » believe that any anisotropy 
in the soil’s hydraulic conductivity will be found be moisture 

N 


ye. The solution was 5 derived by the finite- element method. It was based on the 
combination of a variational principle (Ritz Method) and a weighted residue 
(Galerkin Method) : as described by Neuman et al. (2-4). This method is applied 7 
by minimizing a proper functional, allowing the coefficient of (@/at) to vanish ¥ 
and K*() = 1 at saturation (Eq. 4). At any given moment the solution e 

(x, z,¢) may be described by a complete set of f functions ¢, (x, z) and time- -depen- 


In a finite sum of N gives. 


vectors forming this sp ace: 


MIP = 
4 being the column over is defined. For finite N 


min, and i in our case _ 


The solution can be obtained by feeding ‘the coefficients F . for or , according 
tothefollowing 


This Galerkin method ‘provides: a solution (x, z) at any given moment nt t. The 
.|6h—c[—E dp /att must therefore be determined independently. . For more details 
the interested reader is referred to the original articles (2, socom ote Ne 
‘The grid points for solution are presented in Fig. 5, the upper part of which 
shows the details of the transition layer. The calculation involved the | following | = 


‘program calculating the grid coordinates for each problem, 


i 
a 
} 
hich is the projection of 
| 
| 


_ the boundary condition 
A as The main multilayer calculation in terms of the pressure head, moisture _ 
content, total head, influx, and outflow. It was based on a previously written _ 
& 
Legend 


A Schematic finite element mesh 


ma Vertical Crossection 
_. Scale Distortion 
Transition 


Surace 


FIG. §.—Details of Numerical Calculation Mesh Nodes 


= 


program (3) that was modified and adopted to the present problems. ick adi of 
3. Data processing and plotting as explained in the following. 


Estimate of Sovution Quauity 


The semi-implicit method of calculation with a system of linear equation 
makes possible convergence of the solution. The choice of the time step Ar : 


or outflow of water is expected to introduce a gradual, monotonic change within © 
the flow domain. Excessively large Ar steps cause fluctuations that gradually 
decay toward the true value. A gradual, stepwise change in the rain over a_ 


a 
| 
Je 
| 4 _ small number of time steps can eliminate this problem. We shall not go into 
details of the analysis. The time step was finally obtained by a trial-and-error 7 


4 and after th 


bar) 


in which At, is then ‘compared with At, “until At, < At,,. ‘Obviously, there 


were sensitive points in the field that determined the necessary time steps. __ 
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FIG. Depth Distribution at Different Times: Problem 438 with | 
Transition Layer and Low Rate of Reintall: 
_ The density of grid points must also be related to’ the ra rates of change, the 
- typical time constant being (Ax/K) ~ At. The density actually employed was 
much lower than that dictated by the foregoing criterion. A series of tests — 
_ was carried out to determine the influence of the location of the grid points — 
on the solution. The greatest numerical error was found near the wetting front. 
_ Typical error errors in head were of about 4 cm, being equivalent to a relative error _ 
the rain flux g locally exceeds the hydraulic conductivity K, the time 


_ takes into account the dynamics of the flow at that ime: 

= = 


therefore those involving a thin transition layer (small Az) yond high rain intensity. 
In fact, they called for a relatively ‘greater number of time steps in order . 
_ The part of the problem that presented the | greatest difficulties in computation 
was the appearance of seepage at the soil surface, 
In order to estimate the sensitivity of the Solution to the grid density it was x 
increased horizontally from columns-21 columns, a numerical error 
of only 0.25% at the wetting front (not more than 3 cm of water). The entire = 


was accurate to cm of water. the accuracy 


hr. 
with No Transition Layer and Low Rate of Rainfall 
can be aceaai by doubling the number of grid nodes, involving an approximately 
4 fold increase in the cost for executing the program in an IBM 370/168 computer. 
a grid of 297 nodes was only used. 
_ The accuracy tests were run on several problems involving different geometri 
_ Finally, the moisture content profile at different times (Fig. 6) ¢ and the surface a 
moisture at different times (Fig. 7) were plotted and found to be in good accord = 
with classical one-dimensional solutions of infiltration. The run presented in = . 
_ Figs. 6 and 7 is from problem 438 with no transition layer and a very low 
rate of rainfall (0.01 cm/h). This value can be related | to the conductivity 
¢ “the soil (10°° cm/s). and the final moisture content can be found from Fig. 
4. The results are quite similar to those fee one- “imencioanh infiltration (1). i ‘ 


ot expect a perfect fit, if only because there is some lateral movement _ 


; q of moisture even in a soil lacking a transition layer. This also explains — 


= are two different moisture curves in Fig. 7, one for the hilltop (x = 


0) and one for the vailey bottom (x= 20), 
conclusion, there is no way of ascertaining a shadow of doubt 


what is presently required of the solution is far less than perfect numerical — 
accuracy; it should rather serve as an indication of correct trends. The numerical 
- solution is used as an experimental tool. The conclusions from these numerical 
_ experiments should comply with a number of earlier analytical conjectures and 

_a large voiume of observations (6-9). Tor 
The following ai are lists of problems solved and s ‘some e of the ne 


ano Notes Asout Resutts (Taste 1) 


Column a —Different are given for each wavelength. The 
; - large group for L = 20 has only one amplitude,ie.,2m. 

_ Column 3.—The thickness of the transition layer 8, includes all eight sublayers. _ 
In the first group for L = 20 m, it is either 0.45 m or 0 m (the latter indicating 7 

the absence of a transition layer). 


Column 4. .—Rain intensity P. In the ‘first group “for L= 20 m 

values are 4 cm/h, 1 cm/h, 0.01 cm/h, and (a), the last indicating intermittent 
tain (problem 460). Note that the hydraulic conductivity of the subsoil is 10 = 
cm/s = 0.036 cm/h. Thus the lowest intensity is 3.6 times smaller and the 

highest 270 times greater than the soil conductivity. In some problems the initial 

: rain was increased gradually stepwise due to difficulties in convergance of the — 
numerical solution. Therefore, P gives only the intensity of the main ‘rainfall. — 
It can easily be checked that the product of time T and intensity P does not 
Columns 5, 6. —Average ‘slope A/L and curvature A /L?. ‘They should Id 
an important influence on the flow regime. However, the only results recorded ; 
— are the time and the total rain depth up to the onset | Of saturation 1 and 


Columns 7-9. —Total rain depth D D,,., time T,,,, and the number of calculation 
time steps 7S,,, up to the first appearance of saturation within the profile a 
- all cases involving transition layers, moisture concentration was observed | 
on the concave part. The reason for the absence of saturation will be pare by 
Columns 10-12. —Total rain depth D,,.,, time and the 
~ calculation time steps TS,.., for the spreading of saturation up to the soil surface _ 
and seepage in problems involving transition bs dame In me with no transition 


Int the case of thinner transition layers there i is a toward a decrease 


5 tion of some 60 different problems with various changes in the parameters, 
| he ion eems to ale ormahble he recy tc are reasonable thermore 
sth L is 20 
q 
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result. It becomes manifest at any tain intensity, from low (problem 522) ——o ao 
| intermittent (problem 460) to the highest (problem 815). The rainfall ranges 
_ from 0.01 cm/h-4 cm/h, while the conductivity of the subsoil is 0.036 cm/h. a 
q 7 In soils lacking a transition layer, saturation appeared only when the surface | 
_ was flooded by high-intensity rain (problem 533, D, at = Dycep): It never appeared 


_ with the low rainfall (problem 438). This is exactly as predicted by the classical 


transition’: thicknesses (is cm in problems 752 A and 752. 
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FIG. 8.—Cumulative Rain Depth Necessary for Onset of Saturation or Seepage in 
>. and 5 cm in ‘problem 536), the times or total rains ns necessary ‘for attainment 
of saturation are intermediate between those for thick transition layer (45 cm) : : 
and those in the case where there is no transition layer (problem 533), 
Some of the results for the total rain necessary for saturation are presented bi 7” 


ae saturation over an wide range of rain intensities (problems 522, 
2. The thickness of transition layer has the main effect. For high-rain intensity, — 

_ it is proportional to the total rain depth noasallary for the onset of saturation. — 
3. In the case of thinner layers and very low rates of rain, no saturation — 


J 

q 

. 
mem 

or 
. 
4 
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‘oth lem similar conclusion can abe tute about the occurrence of seepage (saturation 
_ reaching the soil surface). There are two notable trends: (1) At lower rates 
of rainfall the total rain necessary for attainment of seepage is larger; and q 
the effect of a thick transition layer is to increase considerably the total 
rain prior to ‘seepage, provided that its intensity is sufficiently high. If the = 
= is of low intensity, very thin layers will never attain a state of seepage. at) 
= results are quite reasonable. There is a lateral movement of moisture 
that tends to increase the total water increment in the concave valley. However, 
saturation and seepage will actually « occur depends on the interplay 
— several mechanisms. The added moisture spreads over a given soil depth. 
a This depth increases over the course of time. For a given total rain the wetting _ 
depth is greater at a low rate of rainfall. Thus at low rates of rainfall there _ 
i is more time for lateral motion but also more time for deeper moisture distribution. _ 


uckness of wetted 


: = - Thickness o Main rain inten- 

Problem Half cyte, | Amplitude, = ot sity, in centi- 
number in meters  inmeters layer, in meters | meters per hour 
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815 

(438 
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anisotropic or layered soil. Thus, the discharge with thin transition 


: be as high as that with thick transition layers. This is true at least over a 
_ sufficiently long time interval during which the wetting front penetrates beyond — ? 
& thin transition layer. In a thick transition layer the lateral discharge -_ = 
an opportunity to build up over a a period. _ Therefore, thin layers may “Tequire *] 
’ : less rain in order to attain saturation point, but if the rate of rain is ane enough 
they may neverreach saturation, = 
4 Derans oF Frow Recime 
_ Table 2 is a partial list of problems whose solutions will be used in od 
present analysis. Each plot of the results has an identification label including __ 
the problem number at the upper right | corner, the time 7, the rain . intensity — 
-- the half-cycle length i the rain depth D, the n number of time “steps in the 
_ program 7S, the amplitude A, and the thickness of transition layer 8. conil ody 4 
4 The results will be presented together with the development of the analysis. — 
They will include moisture distribution, hydraulic head distribution , integral 
- moisture accumulation, etc. In order to show the details of the head distribution, 


"three vertical sections transition layer have been n enlarged (Fig. 
—— 


— atCla UCDCHUC U Utd | 
‘TABLE 2.—List of Problems Analyzed in Detail 
a 
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At the top Q) in the region of maximum and (3) at the 


‘some iW mm of rain have entered the soil. It is clear that at the wetting front 
the flow is normal to the transition layer. Behind the wetting front, streamlines i. 


; (that would be - approximately normal to the lines of equal head) clearly — 


downstream. ‘This i is most pronounced i in high slope (Section B). a ns : 


dayer, 8 = 0.45 m (problem 522). The lower diagram shows an almost vertical 


Legend 


to 


6 8 10 12 ‘4 
Horizontal Distance X (in meters) 


ina Vertical Crossection 
_ No. Scale Distortion 


= 
N 
3 


flow i in a soil lacking a transition layer, 5 = -0 (problem 438). Both soils received * 
_ rainfall of the same, very low intensity, equivalent to about one third of the 7 
subsoil hydraulic conductivity, and both diagrams represent sections of a " 
after precipitation of about mm of rain. 


interesting to note the ‘elevations marked on the left side and right side of 
the diagrams. If the hydraulic head registered on the 150 g is lower than the 
elevation, the flow is under suction and may be unsaturated. On studying Fig. © a 
‘Wand. 11, one can one can observe: the high suction suction at 


pe nee Af head high rate At raintall and ati 
q 
4 
= 
| 
ent 
teds 
9.—Se 


suction. 

there i is no ny rene or perched water formation anywhere. Nevertheless, there : 
is a lateral flow component both when the rains greatly exceed the hydraulic — 
conductivity and when they are much lower than it. In Fig. 12 one can see 
! that at the onset of the rain the flow tends to be normal to the soil surface 


. since suction gradients are predominant. Later, normal flow i is maintained = 
OW CO 


Within the wetting front. Lateral flow continues for a long, time after 
the rain has (the lowest part of Fig. 12). 


he ) Transition Layer {in ony 


vidi No Seale Distortion 7 


Horizontal Distance = 


T=43 P=4 L=20 
D=n2 TS=40 Az 1045 


otby, 
FIG. 10.—Segments A, B, and c: Lines of Equal Head (Problem 815); “Head in : 


comparing Fig. (13 with previous ones, one notes how a soil 
lacking a transition layer at its surface differs from those with such a layer. 
:. both, the flow starts more or less normal to the soil surface. After prolonged _ ; 

b~ it tends to become vertical in the uniform soil. During drainage time after 7 

_ the rain has ceased, lateral flow seems to develop, but it is of negligible magnitude. : 


the horizontal flow component is to the slope. 
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eheeee avast now ask whether such a simple rule also applies applies i in the case “ae 


to ae the moisture content © over a vertical column | in order to obtai j 
= § O(xzt)dz. Then by conservation of mass 


which P = the rate of rainfall. We assume the lower boundary of z over 


‘whi ch © ‘is ‘integrated t to exclude any deep percolation of water. The values 
of w were obtained by computation, as were In this 
— of Equal Head in the Transition Layer cm) 


T=2620P=0.01 L=20/522 
(D=262TS=270 A=2 5:04 


Horizontal Distance X Stim meters) 


— ou nes of Equal Head in the Transition Layer e cm 
a No Scale Distortion, 


FIG. 1". _—Middle Segment B of a Slope: Lines of Equal Head in Uniform 
Soil (Problem 438) and Soil with Transition Layer — 522) oe 

_ possible the determination of Q(x) by integration over x from the point x = 
0 in which Q(x) = 0. One can now divide this calculated discharge at every 
point by the local slope. ‘Table 3 is such a table for problem 522. Instead of — 
the discharge we have used the weighted average of the horizontal hydraulic 

7 gradient, the latter being exactly proportional tothe former, 


It is quite evident that in the aforementioned case the horizontal gradient © 
very the local slope. This conclusion is far from general. 


more water accumulates, the atere decreases. At the top of the hill, moisture 
depletion will also tend to limit the 


‘al 
| 
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tgrizontal Distance X 
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_ TRANSIENT FLOW 

_ From the numerous results we have chosen two sets to show he ain = 

7 In these figures, four sequences of diagrams (from top to bottom) show four : 
‘Stages of flow within the top soil layer. The vertical scale of this layer has 
been enlarged in 20: I in order to show Ons renee Lines of equal moisture 


Lines of Equal in the er (in cm) 
No Scale Distortion, 


| Horizontal | 


“Sub Seoment > 


T=55) P=0 L=20| 815 | 


7045 


Segment B = <4 “25 


Ft 


> with Transition Layer and High- Intensity Rain at Three Different Times Ft, ey 
are depicted, the moisture content being indicated on a volumetric ‘out, raed 
of saturation are shaded. Owing to the enlargement of the vertical scale, flow 
directions should not be inferred from the lines of equal moisture. Wil Lite ie ; 


_ The first sequence (Problem 522) in Fig. 14 shows early precipitation stages 
; with the onset of moisture build up at the valley bottom. At some stage saturation 
_ occurs. As the rain continues the saturated area spreads upwards and sideways 


thez zone of higher moisture is ; maintained for : a a long period despite respreading 
_ due to downward and lateral flow of the excess water. Saturation and Seepage = 


| 
igh 
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occur despite the fact dia the rain is saliitieas to hie about 1/3 of t 


hydraulic conductivity of the subsoil. The moisture depth accumulation in the r 


concave parts greatly exceeds the average rainfalldepth, 
438 shown in 14 involves the same rate of rainfall but 


1S (problems 388 and 533) the comparisod between soi ils with ‘and 
hout t transition layers is repeated. The previous two problems differ from > 
in their -rain intensity times larger than the 


wit 


of Head in the Transition (in cm) 
No Scale Distortion, 


|Sublayers 


at 


ang] 


(inm 
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Ak 
high and low rates of rainfall. If a transition layer is present, in both cases » 4 
“Z the rain first enters the soil. It then produces horizontal flow components, wa a 
followed by moisture accumulation in concave spots. Saturation initially appears 
__ in such concave spots within the transition layer. Later, it spreads and reaches * 
_ the surface in one place only, where runoff, in the usual sense of a well-defined __ 
leads and conclusion. It has been n argued previously 


= 
| 
& 


can be produced by rains that first enter the ground 

- and then accumulate laterally. ‘Here it is shown that i in the pr presence of a surface _ 
transition layer this is the only possible mechanism. Every drop of rain, be 7 
it of high or low intensity, first penetrates into the soil. The classical model 


. had Details of the flow at the valley bottom c cannot ot be se seen in Figs. 14-15. oe ; 
it has been found that an upward hydraulic gradient is actually formed and — 
= that there is an outflow from n the soil. Moreover, additional rain is able to 


penetrate the soil at this point. » 
_ The ‘seepage i is s capable of ‘Starting an erosive piping 


According to ‘the traditional models, the high rate of rain seseed 


TABLE 3.—Ratio of Horizontal Gradient to Local 
Horizontal distance from - 
_ the hilltop, in meters _ gradient to local slope 


Hilltop 


greater t than the ¢ hydraulic scilleenitte should have produced a v a water excess” 
= ‘soon after its inception, and shortly after the rain has stopped the runoff 


. should also have stopped. Note the type of _— ed obtained. As = 


Let a | to that in all the aforemen-. 

tioned solutions but with one difference. There is an eroded vertical cut at 

the valley bottom, (e.g., on the right-hand side of all the profiles shown in 

: Figs. 14 and 15. The solution remains unchanged up to the time when saturation | 
4 occurs. The exposed boundary at the cut acts $ as an impermeable one when = 
‘under : suction. It is then a streamline and does not differ from the symmetry - 
line assumed in the original solution, oft wom 
_ Seepage flow from the cut will start as soon as saturation is attained. This oe 
_ occurs much earlier than in the case where saturation has to reach an uneroded ~ : 
soil surface. Recall the saturation times or the total rain ‘in necessary 


4 
; gs 7 fails not only in the case of low rates of rain. It is incorrect under all circumstances 
| 
110 times 
= 
Se 
> 
4 
| 
(1) 
| | 
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= the onset of : saturation aie seepage ge (in Table 1). For the p present analysis 


and those necessary for seepage from a transition layer. 10 
sy The conclusion i is that gullies considerably mitra the time e and reduce the 


T8270 


q FIG. 1 4. _ of Equal Moisture » Content St Showing | Stages of Moisture Aewstindetion > 
deep | plough layers, a greater amount of rain is able to penetrate into the soil _ 
before the onset of seepage. From Table 4 one also notes the following: (1) 
The more gentle the slope (and the concavity), the more rain required to bring 
about saturation and seepage; and (2) the thicker the transition layer 5, the — 
more rain necessary for the onset of saturation and seepage. 8 —t*t~*s 
sf seepage forces are the principal erosive mechanism, erosion can commence 


ms 
= 
: 


TRANSIENT wel | : 


‘teen and greatly increased seepage forces. Once erosion has been | 
initiated, the seepage gradients and outflow at the gully sidewall increase 


considerably and cause elongation, widening, and branching of the gully. 
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FIG. 15.—Lines of Equal Moisture Content Showing Stages of Moisture Accumulation 
_ The foregoing offers an explanation for the formation of tunnel erosion so 
_ typical i in many areas. In Israel it is pronounced, especially in the wind-blown 
- Joess soils of the south. It is a well-established fact that the initial bulk density 
of these soils is low. The delicate cohesive structure of the loess can collapse 
_ On saturation. The manifestation of such saturation in the profile is well 


demonstrated in the foregoing. The sites of its occurrence coincide exactly | 


=f 
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SOON as seepage starts olution, the upward 
and outward gradient may appear to be small. However, it should — 
_ | Rein q 
— 
= —— 
= 
Wi Ose where underground tunnels were discovered in tne loess area, Wher 


the soil structure collapses a cavity is produced. Water seeps from the upper — 
end of this cavity and reenters the soil at its lower end. When such — 

- soil is transported the cavity is enlarged (the difference between the initial 
soil and the hydraulically redeposited soil). The tunnel is thus elongated and © 
widened until its top caves in. Often it finds an outlet to a larger truncation. = 
Evidently erosion will progress at an increasing rate once it has been initiated. ™ 
Moisture concentration will start erosion an e erosion in turn 
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FIG. 16. 
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Horizonrat OF Rainwater 
_ Vertical integration of the moisture content, as in Eq. 12, was performed 7 
at 11 grid columns of the numerical solution. The initial value w, at ¢ = o- 
was subtracted from this integral to give Aw (t) resulting from the rain. The © 
average Aw(t) was calculated and finally the ratio Aw/Aw was expressed. In © 
“cases where this ratio is unity, the local increment equals the average added» 
ae. If the accumulation of moisture occurs, the ratio is higher than one. Figs. sl 
17-20 show the results for five different cases. Let us first consider Fig. 17 
for problem 388 (24 cm of rain/day). The curves in the lower part show moisture 
g., at 2.5 h and 2.5 cm of rain there is some > 
Si 
, 


dic 
q ge by Seepage due to Concentration of Moisture in . 
Concave Landscape (Problem 812) | 


(87 


hill accumulation as the is to and the force 
of gravity is still negligible compared to the pressure gradient. Later, hen oo F 


i locally added water in the concave part (on the right) increases re 


Total Rain for 
Seepage 
in | Truncated | No trun- 
| Transition | i layer (smooth 
layer =| sity, P, (gullies), surface), 
plitude, thickness, i i- | 
‘ Lin | A,in 
number meters meters | meters 


522 
A 045 

| 0.15 


enough 4 


even reaches twice the average. _ After the rain has stopped (upper part of the 
diagram, i.e., drainage), the build up of moisture continues in the lower, concave 


and continues in the other parts. relates to a low 
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rate of rainfall (0. 24 cm /day). The relative excess is very significant in the y 
concave part, i.e., more than three times the average rain. Further concentration © . 
continues for | some time after the rain has stopped. Despite deep drainage it : 
is maintained for a relatively long period. After more than a year it is still 
twice the averae added rain. In all the previously mentioned problems, tie caren 


Problem 438 | (Fig. 19) does not heettan a transition layer and serves ion the — 
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FIG. 17.—Horizontal Distribution of Relative Added Moisture Depth at Different Times: 


purpose of comparison. ‘There is very very oe, moisture accumulation ent it is 
largely uphill. Fig. 20 exhibits the behavior at a very high rate of rainfall (40 

mm/h-96 cm/day). It becomes evident that during the rain itself there is less — 


time for lateral movement of moisture and the accumulation is less than that — 


in problem 388. In both of the aforementioned problems the accumulation is 


less than that in problem 522. Later on, after the rain has stopped, there is 
still time for accumulation to develop. wee 
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FIG. 19. aemeaaition Distribution of Relative Added | ere: Depth at Different Times: 
Problem 438 tah he 
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One is s tentatively led to the following important conclusions: () Relative 
moisture accumulation depends, first and foremost, on the total rainfall tative 

not on its intensity; (2) the lower the rate of rainfall, the greater the lateral 
accumulation; (3) accumulation due to lateral flow continues long after the rain 

has stopped; and (4) intermittent rain can generally be regarded as rain at a 

low, continuous rate. 
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20.— Distribution of Relative Added Depth at Ditterent Times: 


We We should ‘recall that lateral migration of rainwater and its accumulation in 
concave parts of the landscape may but do not necessarily produce saturation 
seepage and runoff. At a very low rate of rainfall, moisture accumulation can 4 
occur over a long column of soil. The moisture is then so distributed as to | 
have low local values. High-intensity rains, while causing somewhat less moisture — 

accumulation, have less time to D penetrate into the Thus, they 


2 t=20 
— 


ov jer r the surface. Note that j in reality splashing rain may 
ause higher local concentrations, | on which effect of the surface transition — 
layer is superimposed. The latter can be by far the most important process, 
especially if it is bound up with erosion via seepage. Furthermore, this effect 
has been shown to apply over a wide range of slopes and curvatures. 
The distribution of Aw(xt) in Figs. 17-20 is reminiscent of a negative image 
. of a ‘surface sinusoid showing the higher moisture accumulation in the more 
_ concave parts of the landscape. The moisture will show a high, though not 
: correlation with the concavity, as in the case of the Beer-Sheba experiment 
_ (6). It should also be noted that in the foregoing a reasonable proportionality 7 
but not a perfect one) was found between the horizontal discharge and the © 


ene In view of the relatively limited range of | our | - calculations, they cansiot 


serve as a full proof to this effect. 
_ Many of the conclusions have already been mentioned in the preceding analysis. 
_ he us summarize them here in concluding the entire report, paying special 
Sen to the answers they provide to the questions raised in its first part 


g Lateral flow does occur in the transition layer. 
2. This lateral flow i is roughly proportional to the slope. 
3 It is also roughly proportional to the rain intensity itself. “ye 
4. It leads to the concentration of water in concave parts of the landscape. 
% 5. Saturation will initially appear within the transition layer, in its most concave - 
% part, almost regardless of the rain intensity. 
During the early ‘Stage of a rainfall suction prevails ‘over the force 
gravity y and the flow is normal to the \ "wetting front. It can thus give rise to ake 4 


7. No moisture accumulates in concave spots in a uniform soil during ‘the £ 


ae The accumulation of ‘moisture in concave ‘Spots c continues long : after the f- 

ae. The total excess moisture in a concave - depends more on the pai 


proach 
the upper part of the transition layer as seepage. 
Seepage can form runoff along the soil dip. 
> 13. The fact that all of the rain first ‘enters the soil and _— then may partly 
‘seep out explains certain leaching phenomena as well as erosion. gio! = Ng 
14. The preceding explains why in many areas runoff actually commences 
. only after a certain amount of rain has fallen, irrespective of its intensity. _ = 
15. It also explains why a thick and level plough layer prevents early 1 runoff. 
16. If a gully develops at the soil dip, seepage starts much earlier. This is 
significant in the case of low rates of rainfall. = 


17. For uniform water distribution in a field, the latter has to be moothed 


. Much work relating to all aspects of the new approach to surface hydrology | 
_ is still necessary. We have raised sufficient doubts about the current formulation | 
_ of certain hydrologic phenomena and furnished sufficient evidence in support 
of the new approach to justify further research. It could be directed at measuring 
soil properties at a point, characterization of a watershed as a whole, transforma- 
tion of moisture accumulations into runoff and groundwater recharge, improved — 
algorithms that will make possible the solution of three 
flow problems, correlation wi with present watershed data, etc. oes ble blood! 


c., 
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The following sy mbols are used in this | 

= amplitude of surface elevation; = | off 


ih = length of half-cycle in soil surface elevation; differential operator; hee 
P = rainintensity,incentimetersperhour; 
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total rain time; 
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time; 
volume; 


- = volumetric moisture content; 
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HYDRAULICS DIVISION 
ULTIPORT SLEEVE VALVE DEVELOPMENT 


APPLICATION 


water have presented the need to develop a control valve and energy dissipator _ 
= would be compatible with the aqueduct system. A valve design is needed 
that will: (1) Adequately dissipate high energy flow at small discharges; and i 
: (2) pass the maximum design flow with a minimum of energy loss. The aqueduct 7 
_ would dissipate the majority of available energy in line losses when operating - 
under full design discharge flow conditions and the control valve would dissipate — 
- the majority of excess energy when operating under a low-flow condition. = 
_ Flow control stations on many long aqueducts have generally been limited . 
to pressure head differentials of approximately 15 m. This head restriction resulted _ 
from the cavitation damage associated with the use of butterfly valves when 
used for throttling at high head differentials. If a valve and associated energy 
dissipator could be developed to accommodate pressure head differentials 
exceeding 15 m, the number of * ‘control stations’’ ” required along an aqueduct 
- could be reduced, resulting in significant cost savings. With these needs in 
mind, the Division of Research and the Division of Design at the Bureau of 
s (now ‘Water and Power Resources Service) Engineering and 


The design and construction of long aqueducts to supply municipal and industrial | y 


develop a more versatile control valve and ¢ energy -dissipator. 


_ 'Hydr. Engr., United States Dept. of Interior, Water and Power Resources Service, A 


Engrg. and Research Center, P.O. Box 25007, Building 67, Denver Federal Center, Denver, i 
al ? Retired; formerly Supervisory Mech. Engr., United States Bureau of Reclamation, 


7 * Supervisory Civ. Engr., 


_ Note.—Discussion open until June 1, 1981. To extend the closing date one month, 


a written request must be filed with the Manager of Technical and Professional Publica. 
tions, ASCE. Manuscript was submitted for review for possible publication on March 12, 
1980. This paper is part of the Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, © ASCE, Vol. 107, No. HY1, January, 1981. ISSN 
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Nov. 6, 1979, ne Bureau of Reclamation was renamed the Water : and 
_ Power Resources Service of the U.S. Department of the Interior. The new 
name more closely identifies the agency with its principal 
water and power. Some of the references listed in this paper were prepared 
prior to adoption of the new name; all references to the Bureau of Reclamation _ 
or any derivative thereof are to be considered synonymous with the Water q 
ti _ This paper presents the laboratory results of tests conducted to determine 4 
the discharge coefficients and overall valve performance of two 200-mm sleeve — 
Fi valves. The essential components of a sleeve valve are two concentric cylinders. 
cylinder is stationary and the other serves as a a control sleeve traveling 
ah 

Hydraulic Gradient ovisy evesle ofT 
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- ontal Multijet Sleeve Valve ebay 

| 

across the ports of the stationary cylinder. The first valve tested was 
_ the Bailey Polyjet valve. As can be seen in Fig. | the control for the Bailey . 

valve consists ofa cylindrical s sleeve that travels o overa a a smaller diameter stationary 


: of this valve is fixed and contains the regulating ports; the inner, smaller diameter 

APPLICATION OF Mucriport Sueeve VALVES the quired 


aoe sleeve valve in ‘the laboratory, Fig. 2. The outer cylinder 


_ Previous studies (2) on a 50-mm model of a multiport sleeve valve” and a 
_-vertical stilling well indicated that such a valve could meet the control requirements 


| 
q 
te O3mm _| [Pe 
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Ltd., une a » eubmonged valve similar to ‘the Wanship sleeve valve (4). 


= developed a new y device which could be attached | to their standard sleeve 


dissipation resulting from the email individual jets leaving the valve. 
- The MWD (Metropolitan Water District) of Southern California (5) conducted 
tests on an improved submerged discharge valve (without ports), ‘but found, | 7 
that at pressure heads in excess of 30 m, cavitation damage occurred on the = 


bottom plate of the valve and edges of the control sleeve. To develop a valve 
that could control high energy flows, the concept of flow nozzles was used. 
MWD engineers developed an outer sleeve containing a large number of small 
nozzles and attached it to a sleeve valve under study. The nozzles produce > 
a flow acceleration through | the outer sleeve; thus cavitation does not occur 
the metal flow passages. As the jets exit the valve, cavitation occurs 
_ the water surrounding the valve and not against the flow surfaces. The develop- 
mental stages of the MWD ‘“‘multijet’’ sleeve valve design were summarized 
The multiport sleeve valve is quite versatile from the designer’s point of 
view. The ports can be designed to meet almost any hydraulic requirement. 
It can regulate flows from fine regulation of small discharges to maximum flow — 
with minimum pressure head loss. _ The physical arrangement of the valve « can mf 
be varied to meet a wide range of installation requirements. It is adaptable — 
to various types of operating mechanisms with different mounting positions. — “4 
The sleeve valve can operate over a wide ide range adage heads including = 
a wide variation at any given installation. __ 
“the multiport sleeve valve is designed for use on water pipeline pressure 
systems where the differential head across the valve at low and medals 
_ flows usually produces cavitation damage in other types of valves. What this, 
in effect, means is that the back pressure on the downstream side of the sleeve = 
valve can be relatively small, just enough to keep the valve submerged while 
the pressure on the upstream side of the valve can be quite high. =e. - a. 
“h addition to eliminating valve damage from cavitation, the multiport — 
a valve also offers another major advantage in that, at full flow, head loss through 7 
the entire valve structure can be as small as 1.5 m. The low head loss features” 
can yield significant savings in costs of pipe. The ability to dissipate high- “pressure 
heads enables the designer to use fewer rate-of-flow control stations along an 


controlling high- -pressure head oui ranging fr rom 150 to 300m at one installation. j 
Previously, such high-pressure head flow was controlled using several steps 
of energy reduction to prevent cavitation damage to the control valves. = 


. The flow and energy ase quite upon leaving the valve, 
thus requiring a relatively small energy dissipation structure. ew 


The = of and of cavitation bubbles, 


SLEEVE VALVES mh 


‘which oc occurs dating ‘throttling of high energy flow, can be occur 
in the water surrounding the valve and not against the flow surfaces of the I 
Jet velocity deceleration rate is related to port size. The location of the 
on 
_ “cavitation bubble collapse zone”’ is related to the port exit shape and its relative _ 
. proximity to other ports. Fig. 3 shows the flow pattern developed by a submerged 
jet (7). Albertson (1) and Yevdjevich (11) examined the characteristics of a 
submerged jet. . Albertson described the phenomena as two stages: zone of flow 
_ establishment and zone of established flow. In the zone of flow dae 
the core of the submerged jet is penetrated by viscous shear until the oneal 
velocity begins to decrease. Thus the jet decelerates while the fluid surrounding — 
the jet gradually accelerates. The zone of established flow is defined as that 
this Was the cage. Test | bern 


ubmerged Jet Flow Patterns 
zone where the entire jet becomes turbulent and the center-line velocity begins - 


7 to decelerate. The center- line velocity V,, is defined by Albertsonas y 


(slotted port) 


fn exit port anne jet center line; C, = discharge coefficient of port; and 
AH = pressure head differential across port. ireyents The Gala 

_ From Eg. 2 it is evident that for a circular port, the distance X from the 
port, needed to reduce the jet velocity V,, to a “certain fraction of the jet exit 
q we V,, is directly related to the port diameter D,. . Thus a a reduction = 5 
the port diameter by one-half would reduce the distance X required to produce a 

the same center-line velocity V, it = Oy one-half. Ports as small as 3.2 mm have — 


| ee could be reduced from 51 m/s at the exit port to 4.0 m/s in a distance — 


SANUARY 
_ Hydraulic Laboratory Testing.—The Water and Power Resources Service 
Pom Laboratory’s high head test facility was used for laboratory tests _ 
of two multiport valves. The facility consists of a seven stage vertical turbine _ 
_ pump driven by a 186-kW d-c motor. The rectifying unit and motor speed control 
_ converts a-c into d-c. needed for the motor and provides speed selection from _ 
— 200 to 1,800 rpm. Rate of flow is measured with a 200-mm venturi meter a 
permanently installed 3 m downstream from the pump outlet. A 200-mm, aa 
_ motor-operated valve was used to control the test 
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PERCENT VALVE OPENING « (jest (100) 
BAILEY POLYJET VALVE DISCHARGE COEFFICIENT 


FIG. 4.—Bailey Polyjet Valve Discharge Coefficient bebost 
Ss Bailey Polyjet Valve.—Tests were conducted on a 200-mm test valve provided 4 
_ by the Chas. M. Bailey Co., Inc., 1301 59th Street, Emeryville, Calif. 94608 


= 


| 
(Fig. 1). The purpose of the test was to determine the performance characteristics _ 
of the valve under the two operating conditions previously 
a energy dissipation at low-flow conditions and minimal energy loss at design & 


MULTIPORT SLEEVE VALVES» 
a | ‘The control for the Bailey valve consists of a cylindrical sleeve located inside q 
the annular chamber of the valve that travels over the ports. ‘Movement of 


the sleeve controls t the open port area and, thus, the valve discharge. The flow 


perforated ‘sleeve into. the pipe, which is the same diameter as 
? _ Tests were conducted to determine the discharge coefficient, C,, for val 
a _ Bailey valve. The discharge coefficient curves based on the 203-mm-diam pipe © 
inlet area’ and the 192-mm-diam perforated sleeve area are shown in Fig. 4. 
> Since all 1,835 ports were the same size and configuration, it would appear ze 
_ that the overall valve flow characteristics would be the same as the local flow 
characteristics of each port. This would result in a linear relationship between 
the valve coefficient of discharge C, and the area ratio: A ,.,,/A eae As indicated 9 
on Fig. 4, this was not the case. Test results conducted by the MWD of Southern 
, _ California ia (10) on a similar 200-mm polyjet valve are included in Fig. 4. bischa 
\ a. the control sleeve is ; opened, exposing n more ports, the flow near the 
TABLE 1 —Bailey Polyjet Valve Test Data aie “a 


«> 


; 
4 opening, i ae i rs| in meters Cavitation 


in percent ] rs per second 


“ao = (H, — H,)/(H, — H,); H, = —8.47 m (atmospheric pressure at 
_ and Power Resources Service test facility); and H, = downstream pressure head. 


upstream end of aii perforated sleeve passes the upstream ports in an effort =f 
to satisfy the downstream flow demand. The resultant approach velocity and — 
pressure drop near the upstream ports causes a reduction in the flow rere cal 
these ports. This phenomenon is similar to that which occurs in manifold pipes. _ 
To evaluate potential cavitation damage to the 200-mm polyjet test valve, 
the inside flow surfaces of the p ported sleeve were painted with a white concrete 
curing compound. In earlier tests concrete curing compound proved to be an 
excellent indicator of the location and degree of pitting of a metal surface 7 
resulting from cavitation bubble collapse. Table | presents the test data for a 
the polyjet valve for sleeve openings ¢ of 5, 10, 15, 20, and 30%. Pressure head qi a 
on | the valve was measured 1.83 m upstream (H,,) and 1.73 m downstream 
=. (H,) from the valve. The vapor pressure H, at the ‘Water and Power Resources 
- Service laboratory elevation is equivalent to —8. 47 m of water. Seeds ened 
: __ The high head vertical turbine pump head was designed to deliver a maximum 
flow of approximately 0.17 m° */s. Therefore, the cavitation damage tests were — 


limited to valve openings between 5 and 30% and cavitation index, = (Hy 
shows 


| (219 | 0.066 206 008 | 20 
7 


the resulting paint: removal (small pockmarks) © on ‘the internal flow surfaces 

of the ported sleeve under the 15% valve opening test condition described in 

Table 1. The photograph shows cumulative damage from the 5% opening _ 
to and including the 15% open level. (The scratch lines on the paint in the 4 : 
damaged area were caused by insertion of a mirror for photographs.) a 

_ The cavitation damage shown on Fig. 5 is typical of locations where individual 

cavities implode near the flow surface. The damage » was located in a zone 

which extended 152 mm downstream from the ports ‘to 76 mm into the port 

Se zone, Fig. 1. Although the valve was operated at an extremely low sigma value ~ 


(o = 0. 08) for the 5% test, there was no apparent cavitation ag So fl 


indicate that at the cavitation « damage in ‘the test of this | particular va valve is related 

to the quantity of vapor bubbles produced as well as the pressure head siedanihis: 

_ The amount of paint removal would increase for greater valve openings at low 

. values. However, on long pipelines, the majority of the energy head — 

would be dissipated in | upstream pipe losses at larger valve Openings, resulting — 

in higher sigma values with less cavitation yn potential. It is most likely that the 

- critical point, with respect to cavitation damage for polyjet valves on long — 

- aqueducts with high friction losses, will be in the range of 10 to 15% open. 
= In general, the polyjet valve performed well when operated within the manufac- 


— 
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a sleeve valve designs could not be found that | fully satisfied the design b 
criteria desired by the Service, that is, a valve that will dissipate high energy ; a 
= flows i in a low-flow condition and deliver design flows with a minimum energy 
ve loss at the valve. Although previous ; designs function quite well as pressure- 


reducing valves, most do not emphasize minimal energy me when delivering | 


2 s _ The 200 mm laboratory test valve discharged into a 1 ,370-mm- diam, 1,220-mm- 
long stilling chamber. The basic concept of bud valve and stilling chamber is 
“1. illustrated in ‘Fig. 2. Flow enters the valve | from the high-pressure side and 
_ is discharged outward through the perforated body of the valve into the stilling 
chamber. A cylindrical sleeve located inside the valve body travels over the 


oer section of the valve, pire zg the port area and thus the valve 


6 —Ideal Multijet Sleeve Valve Characteristics 
as a function of sleeve travel for a 200-mm pipeline where the static upstream 


_ pressure head is 137 m. A 200-mm-horizontal multiport sleeve valve with nozzles et 


and slots for ports is proposed. As the valve is opened and the po port area slowly : 

increases, the valve discharge increases linearly with sleeve travel. This increase 
in valve discharge results in a reduced upstream pressure head due to friction - 

_ losses in the long aqueduct. When the valve sleeve has been opened 230 mm, — 
it has completed the function of a pressure reducing valve and assumes the 


due to ‘the low available pressure head across the valve. When the valve port 
- area equals the 200-mm pipe area, the pressure head differential, AH, —_ 
decreased to 1.1 m across the valve, resulting in a discharge coefficient C, 4 


| 
| ¥ 
end of the stilling chamber, 
The concept of a linear relationship between sleeve travel and valve discharge 
_ was also sought as a desired characteristic in the new design. Such a valve — 
would provide better control characteristics on long aqueducts where water 
wy ‘hammer presents a potential problem. Fig. 6 shows ideal valve characteristics _ 
| 
| 
. role of a low-head-loss control valve. At this point the valve port area increases ¥ 
idlv with sleev j iminishine in se in discharge 


— 


4 


— 
in in which Q= valve area; and AH = head dif 
, across the valve. The key to the multiport design is the proper placement of 
the ports along the several turns of a spiral that extends along the length of = 
valve to produce a nearly linear ‘Telationship between the re 


“corrected for a pressure differential, AH, between ‘the: upstream of the 

‘sleeve valve, H,, and the downstream 200-mm pipe flange, H,. Eqs. 4 and 

were used to the through | the and slots, respectively 


in whi ich QO, = nozzle discharge; Q, = slot discharge; A = open nozzle area; 
_A, = open slot area; Ah = pressure head difference between valve and stilling © 
chamber; A, “= pipe inlet area; v= Vin - Q, and H = 
_ When the nozzles alone are exposed, ‘oe velocity term on the right side 
of Eq. 6 is negligible since the velocity head is small with respect to the total 
available energy, H. As the control sleeve opens, the total energy term decreases _ 
due to friction losses in the upstream aqueduct, and the velocity term increases. 
Thus the nozzle and slot coefficients, C, and C,, decrease in value to reflect — 
the phenomena of a manifold where the majority of the valve discharge is 
released through the downstream ports of the valve. arte aoe 
_ Eqs. 6 and 7 are similar to the equations presented by Vigander (8) and ‘ 
Enger (3) dealing with large diffusers and manifolds. They are empirical equations — 
_ based on the results of laboratory studies conducted on the 200-mm test valves. = 
4 The investigation included studies conducted on the two-ported (nozzle and 
a - slots) valve sleeve shown in Fig. 7. The sleeve travel versus port area relationship 
; : curves for the two configurations and the polyjet valve are shown in Fig. 8. a 
The Bailey polyjet valve has a linear relationship between port area and sleeve 
_ travel. The area for the slotted port configuration increases at a somewhat 
‘slower rate initially. The area for the valve with nozzles and slots increases = 


very at the start of the sleeve travel and then once 


= 

[ 
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the slots are exposed, yielding a near ‘Telationship between sleeve travel 
The us use of nozzles or orifices for ‘the: discharge ports should be carefully 
_ considered. The present spiral arrangement of the ports results in the pa 
, "blockage of some of the ports when the valve is used to control the flow. 


design be considered for ports 19 mm in diameter or larger oc 
at -~ = for the head loss coefficient K are plotted with respect to percent: valve 


that « an orifice design [similar to ore shown in Fig. 7(5)] instead of the ‘ct 


« \2.7mm siots 
6-7 


—Reterence Pere 


2.7mm 
half) 


gar sia SLOTS 
Sketches of Ported Valve Sleeves 


ad 
opening (100 x port area) 9 for both port 
shown in Fig. 7. The head loss coefficient K is based on the pressure head 
_ differential, AH, between the upstream valve flange and the 200-mm pipe flange 
es of the stilling chamber. The loss coefficient, therefore, includes _ 
the total system loss for the control structure. The Bailey polyjet valve data 
Computer Model.— Presently, multiport sleeve valve designs include a large 
number of apertures in the valve body. To reduce manufacturing costs and 
to obtain a better flow control characteristic, the Service’s sleeve valve design 


utilizes the fact that pipe friction head losses can be epgeenimated for various 


Stipe Uld De Conducted LO dclermine tne Pressure 


flows by using the equation H on H, = in which K = constant for a 
pipeline and Q = CA V2gH_ Using this approach, a computer program (2) | 
(U.S. customary units) has" been developed to locate and space nozzles along = 


alll 


Bailey polyjet 


O-Sleeve with slots (USBR) 
_4-Sleeve with nozzles and siots 


AG. 
FIG. —Head Loss Coefficient 


a - spiral around the valve body ouch that a linear change i in flow is maintained 
- through the zone of the nozzles. As the control sleeve is retracted to expose — 
_ more apertures, the ee and/or diameter of the nozzles is increased until : 


q 
= 
toe 
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‘the largest diameter nozzle is equal to the wall thickness of the _— body. — 

When a point is reached where the computed number of nozzles cannot meet a 

_ the minimum spacing ‘Tequirements along the spiral i ina quadrant that is being 

incremented, then the program m automatically changes its output from nozzles | 

to a preselected number of longitudinal slots. The length of the slots is computed _ 

_ by short intervals such that at the final length, full flow is achieved: that is, Ml 
- when the control sleeve is retracted past the end of the last slot [Fig. ‘1a)). 


The area, » Aq, of each set of ina quadrant the spiral 


V 2gAH. 
4 The linear change in flow is prorat by setting a uniform mete of went between — 
_ spiral peaks; hence the change for a quadrant can = calculated. d. Rewriting 
1 in which Q,, = the theoretical flow (in cubic feet per second) through the nozzles 
q based on the control sleeve travel from the reference index line [Fig. Ha); 
— C,, = 0.94, a value which Temains constant until slots ar are ¢ included i in n the Output; — 
2g = 64.32 ft/s’; 
quadrant being (For ‘each “quadrant, AH is reduced by the 
tot pipe friction for the reach being considered. Since Q is a function of | 
distance from the reference index line and H, = KQ’, then H = Hyinic — 


KQ’); and A, = = the total amount of nozzle area in square feet et required to 4 


iad From this relationship, A, for - quadrant being incremented can be obtained 
by subtracting the total area, A,, of the nozzles in the preceding eaten 
frm A,. Because A, will not exactly match with the sum of the areas of the 
 nonales that can be used in a quadrant, the small excess area is added to the - 
calculated area for the next quadrant before its number of nozzles is } computed. = 
_ The physical model test program was used to evaluate the empirical co coefficients od 
used and has produced good agreement with discharges indicated in the computer — 

: - model. The linear change in Q is maintained until the control sleeve is retracted | 
to the zone of the slots. At this point, flow is approximately 82% of the design 
flow. The end result is a valve which can achieve design flow with a much 
smaller number of apertures. It is hoped that physical ‘model testing can be — 


continued to evaluate larger design flows where two or more multiport sleeve — 
valves would be required in a flow control structure. queens 
_ Mechanical Considerations.—The heart of the sleeve valve is two concentric ; 

_ cylinders. The outer cylinder or body is fixed and contains the regulating ports; 


the inner, unperforated movable cylinder | is the control sleeve. fits closely 


= "providing regulation of flow. Minimal practical clearances a are provided between = 
_ the body and sleeve to control leakage. No seals or packings have been used 
7 on the sleeve. A guard valve upstream of the sleeve valve provides tight shutoff 
when no flow is desired and the sleeve valve is closed. If minimal leakage | 
is important, piston rings are provided on the sleeve. The Metropolitan weter 
- District of | Southern California designed their valves for tight closure (5, > 
Pockings are proviies on the end of the sleeve which — not move across — 


j 
| 
q 
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1e ‘regulating ports and a fixed se seal is mounted in tee te body to seal the other 
dof the sleeve when the valve isclosed. 
The valve ports can be varied in size, type, and spacing to meet a wide | 
ange of requirements. . According to MWD round nozzles are limited in size 
as their diameter cannot be larger than the body ‘shell thickness because of — 
changes which occur in jet characteristics. Because of the limitation on size, 
a large number of nozzles is required. Ports can be machined as slots with 
the same taper as the round nozzles. The width of the slots has the same | ~— 
limitations as the diameter of the nozzles with respect to cavitation damage, 
and the length is limited by the physical strength of the body shell. Because 
of the taper, these ports must be reamed from the inside of the cylindrical 
body. This is a difficult and expensive machining process. A 200-mm valve | 
; s considered the smallest that can be machined from the inside. 
In river or reservoir outlet work installations where the water may carry 

ebris large enough to plug the nozzle-type ports, ports that are designed as — 
sharp-edged orifices | are utilized. These ports have parallel sides for the first : 
1.5 to 7 mm of depth from the inside of the body [Fig. 7(b)]. To eliminate 
any possibility of cavitation on the sides of the port, the remainder of the 
az The size and shape of the ports can be specified to suit specific requirements. 


48 


size e of a port is limited by the physical ‘strength. of the body. shell. However, . 
the larger the port, and in turn the jet, the greater the outward distance — 
(perpendicular to the valve body) that is required to dissipate the energy of 


< the jet. For this reason, the ports should be kept as small as practicable consistent 7 
with the s size of expected debris that must | be passed. 


As noted ‘previously, valves. used in 1 long municipal- -industrial pipelines gain 
_ benefits from linear discharge characteristics even though the effective energy | 
head at the valve drops rapidly as the discharge increases. It has already been 7 
_ shown that this requirement can be easily met by varying spacing and size 
of the nozzles and, when effective head drops sufficiently, changing to slots. 
This. general principle can be used for other discharge characteristics. As an 
example, valves requiring fine regulation of small off-season discharge can have 
a few small ports as required for the off-season releases. The remainder of | 
- the ports can be longer and spaced closer together to give needed equate 
of normal releases with minimum sleeve travel. As port size increases, it is 
- necessary to place ports diametrically opposite to maintain hydraulic balance 
_ The sleeve valve lends itself to many y physical arrangements. The most efficient oa 
: and straightforward is the horizontal valve of the general type used in the Service — 
laboratory test and shown in Fig. 2. In this type, the stem attached to the | 
control sleeve is extended downstream through packings to the valve operator. __ 
Thus all moving parts of the sleeve are downstream of and unaffected by the __ 
water flow. The sleeve is hydraulically balanced. The only unbalance is the 
hydraulic load on the valve stem. Since this stem load and friction are the 
only loads to be overcome, the operating mechanisms can be relatively small. — 
Sleeve valves are adaptable to many types of operating mechanisms. The 
Service h has used a type v with torque- switches, hydraulic 


= 


| 


and even a en a hand- gear lift for a a remote installation without p power. svie 

A vertical arrangement has been used for some valves with the flow being 

discharged downward into a square stilling well (Fig. 10). This is a neal 

a "arrangement for valves | with large ports as the distance required to dissipate 

the jet energy is ‘much greater. The bottom of the stilling well and the — 

a walls should be steel lined and the well sized to allow jet velocities as high 

as 12 m/s to impinge on the steel liner. The stilling well will very effectively 

“the ino the remaining energy of the jets. However, with this arrangement 


__ the incoming water flow must pass through a 90° elbow and the valve operator 
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FIG. 10. Sleeve Valve in Vertical Stilling Well 
is mounted on top of the elbow. Thus the control stem is in the water passage 
and the water must flow through the inside of the control sleeve. The valve 
_ interior must be streamlined t to impede water flow as little as practical and 
to > provide adequate support to: the stem. The Service’ S experience has shown | 
that vibrations can cause fatigue and breaking of a stem. Flow velocity through — ~f 


the elbow and sleeve must be kept low and, when at maximum valve opening, . 


_ From the mechanical standpoint, sleeve valves are fairly simple devices. 


| 
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‘sive. The body, sleeve, and: stem must tbe made fre rom corrosion-resistant mustache. 
The sleeve and body must be carefully and closely fitted. As noted before, ~! 
the ports are difficult and expensive to machine. The stainless steel body also also 
‘Prototyre INSTALLATIONS AND ExpeRiIENCE hee od 
fiw Hove oo ogning ot Shae 
D> A 356-mm-diam prototype horizontal multiport sleeve valve has been — 
_ installed in the Deep Red Run rate-of-flow control station near the end of the i 
_ Frederick Aqueduct of the Service’ s Altus Project in Oklahoma. The valve _ 
is designed for Q = 0.23 m’/s with a static head at no flow of 70 m. From — 
the Frederick regulating tank which is located upstream of the valve, there 
2 a gravity flow-pressure pipeline consisting of 5.52 km of 610-mm pipe, 6.4 © 
km of 458-mm pipe, and 3.8 km of 406-mm pipe. The multiport sleeve valve 
has allowed the designers to use one rate- -of-flow control station in lieu of 
four which would incorporate butterfly valve type rate-of-flow controllers. It — 
2 estimated that savings resulting from elimination of three flow control stations hl 
are well above $200, 000. Field testing of the prototype valve is spunuel for 
Several large (1,372 mm) multiport sleeve valves (Fig. 10) have designed 
for control of flow in vertical stilling wells. These valves have circular or 
elliptical-shaped ports with diameters up to 305 mm. Installations are ed : 
for the Mt. Elbert Aqueduct in Colorado and Sugar Pine Dam River a 


Conciusions oF 


a 


1, Laboratory tests were conducted on a 200-mm polyjet valve. The valve 
_ performed well when operated within the manufacturer’s suggested pressure 
head ranges. However, the multiport sleeve valve design will be used by the | 


Water and Power Resources Service for those applications where operation 


a Physical and mathematical models have ni used to develop a horizontal ; 
multiport sleeve valve design for use asa flow controller on long aqueduct © J 


es The design of the multiport sleeve valve is _ versatile and can be 


in a wide variety of flow regulation situations. 
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= port area; 
rd = pipe inlet area; 
quadrant port area; 
= 
slot port area; 
total port area; art 


nozzle discharge coef 
= slot discharge coefficient; 
circular port diameter; ‘ 3 
pressure head differential across pens, 
pressure head differential between valve and stilling chamber; I 
nozzle discharge; 
slot — 
Q/A,; 
= ‘jet exit 
jet center line velocity at 
distance from exit port along jet center line; and 


= cav x. 


“al Colorado State University, Fort Collins, Colo., Aug., 9-14, 1970, (published 1971). "an 
Miller ‘‘The Submerged Discharge Valve.”’ Glenfield Gazette, No 229 
following symbols are used in this paper; 


i. 


| 
|= 
| 
| = 


9810 


NS 


Note. .—Discussion open until June 1, 1981. To extend the closing date one month, 


a written request must be filed with the Manager of Technical and Professional Publications, — 
ASCE. This paper is part of the Journal of the tye Division, ao of 
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Bint sate a place within ASCE for - publication of technical ideas that have not advanced, > pe 


as yet, to the point where they warrant publication as a Proceedings paper in a Journal, the © 
of Technical Notes | was the Board of Direction on October 16-18, 


 ¢ wl Professional Publications, ASCE, 345 East 47th Street, New York, N.Y., 10017, ad 
bee a request by the author that it be considered as a | Technical Note. ms ea ee 


| If the Division or Council approves the contribution for ern i 

to Society Headquarters with appropriate comments. = 7 
4. The technical publications staff will prepare the material for use in the omnes _— ’ 
of the Journal, aftez proper coordination with the author, 


_§. Each Technical Note is not to exceed 4 pages in the Journal. As an approximatio: 


d an full manuscript page of text, tables, or figures is the equivalent of one-half a Journal 
6. The Technical Notes will be grouped in a special section of each Journal. a_i a. 
4 £ ‘Information retrieval abstracts and key words will be unnecessary for Technical Notes. — 
* ‘The final date on which a Discussion should reach the Society is given as a footnote 
Technical Notes will not be included in Transactions. 
10. Technical Notes will be included in ASCE’s annual ont cumulative subject and author 


The manuscripts for Technical Notes must meet the following requir 
1. Titles must have a length not exceeding 50 characters and spaces. 2 
_ 2. The author’s full name, Society membership grade, and a footnote reference stating present 


7 ‘employment must — on the first page of the manuscript. Authors need not be Society 
3. The manuscript is to be cubmined as an 1 original c« copy (with two ~~ that 
double-spaced on one side of 8-1 /2-in. (220-mm) by 11-in. (280-mm) white bond paper. 


.. All mathematics must | be typewritten and special symbols | must be Properly identified. 


5. Standard definitions and symbols must be used. Reference must be made to the lists 


“published by the American National Standards Institute and to the Authors’ Guide to the 
% 6. Tables must be typed double-spaced (an original ribbon copy and two duplicate copies) 
7 on one side of 8-1/2-in. (220-mm) by ll-in. (280-mm) paper. An  engtenation of each table 
7. Figures must be drawn in black ink on one side of 8-1 /2-in. (220-mm) by L-in. | (280-mm) 
paper. Because figures will be reproduced with a width of between 3 in. (76 mm) to 4-1/2 
in. (110 mm), the lettering must be large enough to be legible at this width. Photographs must 


be submitted as glossy prints. Explanations and descriptions must be made — the text 
8. References cited in text must be typed at the end of the Technical Note in alphabetical 


9. Dual units, i.e., U.S. Customary followed by SI (International System) un units in ance 


should uld be used be used throughout the the wager. 
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RE-EXAMINATION OF NIKURADSE ROUGHNESS DATA © 


William R. Brownlie," A.M. ASCE 
ot gigh dP (11 @ 23%] 41) tolq ows ods 

ob b tud lapigito, ad) ni awode ton 2 ai 
Two sets of flow resistance data are “commonly u used in the evaluation of — _ 

friction factors for pipes and open channels. The data compiled by Colebrook _ 
and White for commercial pipes were used by Moody to construct his well 
known friction factor diagram (3, Fig. 5. 32). A similar diagram based on the 
data of Nikuradse (1) for sand-roughened pipes 2 appears in most texts of fluid ; 
mechanics (2, Fig. 108 and 3, Fig. 5.31), however, with a much more limited — 
+ range of relative roughness and Reynolds number than the Moody diagram. bas 
‘While the Colebrook and White data are appropriate for commercial pipe 
applications, the Nikuradse data, with its sand roughness, may be more applicable __ 
for problems involving open channels with uniform- “sand beds for which grain 
‘friction factor is required. This note describes an inconsistency in the ; original — 
presentation of some of the Nikuradse data and provides a Moody-type ‘diagram 
with some engineering applications for a range of the data believed to be valid. 
data are reviewed here because they appear in classical texts of 


- mechanics for engineers (¢.g., 2, 3). 
4 


‘The experiments reported by Nikuradse were conducted using pipes with 
diameters of 2.474 cm, 4. 94 cm, and 9.94 cm. Roughness was created by gluing 
uniform sands to the pipes. In all, five sands were used, with mean diameters — 
ranging from 0.01 cm-0.16 cm, to give six values of relative roughness (grain 
_ diameter over pipe diameter). Uniformity of sand grains was created by sieving, — 
resulting in a typical geometric standard deviation of 1.02 for the grain-size 
distributions. Measurements in the pipes were taken using an approach length — 
of approximately 40 pipediam. = 
‘The data has traditionally been presented graphically in two different forms — 


following the original presentation of Nikuradse (1). In the Moody-type form, _ 


_ friction factor is plotted against Reynolds number on a log-log scale with a 


b different curve and set of data points for each of the six values of relative . 
roughness. In the alternate form, by transforming the plotting coordinates, the — 


“— Grad. Research Asst., W. M. Keck Lab. of Hydr. and Water Resources, Calif. Inst. 
Note. —Discussion open until June 1, 1981. To extend the closing date one month, 

a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on April 25, 1980. 
This paper is part of the Journal of the Hydraulics Division, Proceedings of the American | 
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q _ Fig. 1 shows data from 90 runs randomly selected from the 362 that are a 
‘ative. The figure also shows the Colebrook transition function upon which 
- the Moody diagram is based. Since the e equivalent sand roughness of the Colebrook 
and White data was calibrated to the Nikuradse data in the fully rough — 
_ the two curves converge to the same asymptote on the right side of Fig. 1. od 
An inconsistency in the original data presentation can be seen by comparing 
the two plot types (1, , Figs. 9 and 11) with the data tables. The data in the + 
tables cover the range of parameters shown in Fig. 1; however, all points plotted 2B 


; on the original diagram do not appear in the tables. Conversely, all of the : 
data in the tables are not shown in the original diagram, but they do conform a 
closely to the curve in Fig. 1. On the other hand, the Moody-type diagram — 
Shows de data wi with h Reynolds t numbers as low low as as 500 whereas the lowe lowest Reynolds _ me 


2 


Nikuradse Data and Colebrook ond white 
Transition Function (about 25 percent of Published Data are Shown) 


g for Reynolds numbers greater than. 10, 000. Finally, the “unpublished 
are somewhat suspect because they show a smooth transition from turbulent _ 
- to laminar flow occurring at a Reynolds number of about 2,000, for all given 
- values of relative roughness. Such a condition seems walikely & due to the nature 


: ‘The Moody- type flow resistance. chart shown in Fig. 2 was derived from | 


the curve fitted to the data points in Fig. 1. Although there are inconsistencies 
in the original diagrams, the experiments appear to have | been carefully conducted 
and the data in the tables are reasonable. Reynolds r numbers lower than 10, 000 “4 
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«ct is hoped that Fig. 2 will be a useful and accurate tool for engineers. The 
- hart can be used for side-wall corrections as well as for separating total resistance: 
_ into grain resistance and form resistance. For open channel flow calculations, 
Pipe diameter D should be replaced by 4r in which r r= hydraulic radius. ci 


_ Fig. 2 is based on three equations which apply to different domains along “oo 


+2log~ 


4, (ios 0.5 < log ~—< 2. 0 
— -2los .. for > 2.0 


= k, 
sand g grain ‘to grain R = Reynolds 


# 
i, NIKURADSE SAND GRAIN N ROUGHNESS _ 


FIG. 2. —Frietion Factor Diagram, | for Pipes of Diameter, D, or or Channels c of Hydraulic 
and A, = empirical constants. . Eq. 1 is for smooth pipes, and relative roughness - 
can be ‘removed by factoring both sides of the equation. Eq. 2 was fitted by > 
Md the writer to the transition data from the smooth to the rough regime, with 
the coefficients A, through A, defined as 1.3376, —4.3218, 19.454, —26.480, 


(16. 509, — 4. 9407, 0. 57864, respectively. Eq. 3 describes the - rough regime — 
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_ Fig. 2 can be used to eoctname a side-wall correction for flow at a given 

z. in flumes with a known friction factor, f, and roughness, k,, using a procedure 

- analogous to the smooth-wall procedure described by Vanoni and Brooks (4). 
the derivation in ‘Ref. 4, the following equations can obtained: 

| 


) 


‘in which p P = wetted | perimeter the subscript | b denotes bed, and the subscript 


The procedue fc for usii using 2 to to calculate and r, is as follows: 

Plot Eq. 4 on Fig. 2 as a b anniate Tes with a | slope of ‘1 in log units, 

“a an intercept of 0.01R/f at hoes 0.01. The desired values of f,, and R,, 

will lie on this line. | 


2. Pick a trial value of r, and compute 4r and from 


should converge after two or three interations. 


4. The quantities f, and r, can now be calculated from 
Grain Resistance 


some open channel flow it is is often desirable to separate 
_ resistance from bed-form resistance. Two procedures are possible for separating | 
_ the bed shear stress into its two components. Either the slope may be broken | 
into components or the hydraulic radius of the bed may be broken into components. 7? 
_ Vanoni and Brooks (4) have presented a graphical solution of the Einstein-Bar- 
barosa approach which divides the hydraulic radius into two components. Fig. 
_ 2 could also be used to carry out this procedure by applying a technique si similar 
to that of the side-wall correction procedure just described. However, a more 


convenient and perhaps more conceptually reasonable approach is to ie 


the energy slope intotwocomponents, 
‘The ey equations ca can be used with Fig. 2 to perform this procedure: — 
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in which S’ and f, = the energy slope and bed friction factor, respectively, 
‘resulting from grain resistance; and S” and fi = those quantities resulting from — 
form drag, ‘for a flow with a given velocity and bed hydraulic radius. The 
"quantity can be determined directly from Fig. 2, given R, and 4r,/k,. The 
‘remaining quantities can be calculated from Eqs. 8, 9, and 10. 


Special thanks to A. Massengale. neo ai dot vit 
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BED Accra 


The sediment ‘in natural streams evolved over 
time are in “‘equilibrium’’ and produce no objectionable scour or deposition. - 
_ However, many factors, both man made and natural, can contribute to the — 
_ imbalance between the supply of sediment to and the sediment transport capacity — 
of a river, which in turn can lead to degradation /aggradation along the river a 
reaches. The channel bed along some reaches of the Missouri River has been 7 ‘ 
progressively degrading and aggrading since about 1960 (3). Reliable, quantitative _ 
estimates of the future rate and extent of long- -term men and aggradation — 


“and experimentally and recently reported by Soni et al. (4). When rib! Re 
results failed to agree with their experimental data, they empirically modified — 
the value of the aggradation coefficient to achieve conformity between predicted 

and measured bed profiles. The nonconformity between the analytical and 
experimental 1 results were attributed to the several simplifications made in the — 

_ analysis. _ The writer, in a discussion of their paper (2), pointed out an error 

in the boundary condition of their analytical model. ‘This paper presents ; the 
analytical solution of the aggradation problem for a more appropriate boundary : 

condition. The analytical results are with the laboratory experimental 

THEORETICAL ANALYSIS 

~ Consider a uniform flow with a mean flow ‘velocity U, and mean flow depth | 
h, in a wide rectangular alluvial channel as shown in Fig. 1. The equilibrium 
sediment transport ‘‘rate’’ under the uniform flow condition is equal to G, and 

_ the bed slope is S,. Let the sediment supply rate at a section x = 0 be increased 

bya constant rate | AG. - The bed level for x = 0 will rise with time due to 

; increase in the sediment load. The differential equation for the bed level obtained 
es Prof. in ees Div., and Research Engr. in Inst. of Hyar. Research, ae 
University of lowa, lowa City, lowa 52242, 
 Note.—Discussion open until June 1, 1981. To extend the closing date one month, bi 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on May 8, 1980. @ 
Tis paper is part of the Journal of the Hydraulics Division, Proceedings of the American 

a Society of Civil Engineers, © ASCE, Vol. 107, No. HYI, January, 1981. ISSN cone ma 
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a» 
= ed 


depth of deposition at time at a distance x x from ei 
- origin; N = the Porosity of the sediment; and b = the exponent in the sediment — 


in which G = total sediment transport rate at any cross section in absolute P 
_ volume per unit width of the channel; V= mean flow nes and a = constant. 
constant. 


[= of sediment ‘injection = 
| 


initial a and for the problem under considera- 


which AG = the volumetric sediment ‘supply per r unit 
The explanation for the boundary condition (Eq. 5) is given in Appendix I. 7 


The forms of the governing equation (Eq. 1) and the initialand boundary conditions _ 
(Eas. 4 and 5) are exactly similar to those of heat conduction i in a semi-infinite SB 


“equation for the prescribed initial and boundary enudidens is(1): 
q 


» 


_ The bed elevation at x = = 0is is given by sat - ta norhag 


wall 

q 
a 

| 
| 
| 


a 


ch con be rearranged i in nondimensional form as a 
Ma 


bed profile i in n the ‘nondimensiona form i is 


The analytical r are the experimental data of Soni 
: al. (4). Since the large changes in the bed elevations which occur for ae 
values of AG/G, can be more accurately measured, the experimental data for 
_ AG/G, greater than one is used in this study. The comparison of Eqs. 9 and 
ae with the experimental data is shown in Figs. 2 and 3, respectively. The -% 
¥ a between the analytical and experimental results is fairly satisfactory. “ 


The relatively large ‘scatter = of in Fig. 3 is due to 


For explanation of 
symbols, see 


| 
‘Furthermore, the length of a 
J, on the assumption that aggradation 
ere z/z, = 0.01 (from Eq. 10) is given by 


Hing water quality” hes 
pletion in meters Sx 10° per second 
3.00 0.385 


8: 
ver 


«FAG. 3. —Comparison of Eq. 10 with Experimental tor Bed Profiles 
_ possible errors in the measurement of small depths of deposit near the end 
by ibe and Bor which uve sedunent yighd 
The aggradation in a stream due to increase in sediment load can be analytically ' . 
determined from Eqs. 9 and 10. The analytical ‘results are in agreement with 


ons 


s The volume of sediment between the initial bed profile (t¢ = 0) and the cat 


"profile at time ¢ (see Fig. 1) must be banal to the volume of additional sediment 


time derivative of E 

ied 
AG= ) 


ox 
Eq. 14 reduces to Eq. 5 waive: pol benimdue sew J 


| 
Sub (at) from Eq. land integration gives 
| 
| | 
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- Section 208 of the Federal Water Pollution Control Act of 1972 has focused: 
s attention on meeting water quality orga through management and control 


proce 
subsequent effects « on water quality. These include the ‘prediction of 


_ yield from agricultural lands, construction sites, silvicultural operations, mining 
activities, and from other existing and changing land use practices. 
7 Efforts to quantify sediment in runoff have centered largely on statistical i 
analyses of measured sediment yields from experimental plots and watersheds. 
_ from these studies, which reflect the combined linkages between land — 
use, rainfall-runoff, and soil detachment and transport aaline used to determine — 
‘ soil loss and sediment yield equations applicable only under long-term, reasonably 
=" hydrologic conditions. Among the better known equations of this type 
7 is the Universal Soil-Loss Equation developed by the United States Department 
of Agriculture, Agricultural Research Service for cultivated soils" cast « of the 


Rocky Mountains. In the basic form of this equation, the long term average 


cover factor, erosion control practice factor, and soil erodibility factor. Typically, : 
when used to predict nonpoint loadings, the Universal Soil Loss Equation times 
a sediment delivery ratio is used to estimate the sediment load or yield, and 
source loadings for various s pollutants are then generally assumed ‘Proportional 


or histographs of various water quality constituents are ‘determined. ‘However 
broad or detailed the output from such nonpoint load generators, the estimates _ 


- remain limited by the assumptions and conditions for which the sediment yield = po 
equation was developed. However, recent progress in sediment yield modeling 
has improved predictions and has added some ‘deterministic predictors to the — 


It is recognized that a variety of constraints for both designated an and non~ desig- 
= areawide waste treatment management studies often precluded more than. 7 
a cursory analysis of sediment runoff from agricultural and cultivated lands. 
Still, the environmental effect of erosion and sediment runoff is potentially 
significant. Fortunately, Section 208 planning is a Continuing process and 
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_ from which improved predictors of sediment yield from agricultural lands can 
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NOTES, 


FLOw INVESTIGATIONS IN Two PuMPING Pits 
_ Model studies with pumping pits te predict the possible flow conditions in— 


a prototype have become increasingly important due to use of more ‘pumps 


detail the associated hydraulic problems in model studies, and stressed the 
need for systematic studies and publication of such results. The hydraulic problems 
_ that affect the performance of a pump are surface and submerged vortices, © 
_ air entrainment in the inflow, and swirl and asymmetric flow in pump intakes. 
In this study, relationships between \ hydraulic problems and various dimensie- 
less groupings, and Reynolds number effect on the determination of swirl angle - 
and the axial velocity distribution in the pump intake are investigated with 
two simple pumping pits. Also, a quantitative method is emengee to qranins 
“The: perfomance of a pump represented by i its 's efficiency (a) ii is affected by 4 
the adverse effects in a pumping pit (2) and this efficiency will be different a 
¥ from the design efficiency (n,) of a pump. If an ideal pumping pit is designed, 
then both efficiencies would be same and the ratio (yn/7,) leads to 


identified to affect the problems: in a pumping pit, the 
pre. 


be = Sei dimensions of pumping pit; yo = - velocity of inflow; p = kinematic 
viscosity; p = density; and t = surface tension of fluid. 
_ Therefore, for complete similarity of a geometrically similar ‘mote all three 
numbers have to be satisfied. The Froude number can be satisfied quite easily, — 
g Sr. Lect., Faculty of Engrg., Univ. of Singapore, Kent Ridge, Singapore 0511, seat 
_ Note.—Discussion open until June 1, 1981. To extend the closing date one month, - 
a written request must be filed with the Manager of Technical and Professional Publications, 
_ ASCE. Manuscript was submitted for review for possible publication on April 8, 1980. 
= paper is part of the Journal of the vol WO Division, yey: of the American a 
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=(wV/t),, “equal- velocity’ simulation is obtained if the 


same fluid is used for the model (m) and the prototype (p) (1). It is rather 
difficult to have Reynolds number simulation and thus a critical value of 5x 10° 
is chosen to simulate flow conditions at the pump intake. - 
PumpinG Pits anp EXPERIMENTAL Setup owe | 
‘ee pumping pits are given in Fig. 1 where continuous lines show 4 S| 
4 the original pumping pit and broken lines the modifications. The dry dock pumping 
_ pit (DDP-pit) has one vertical axis and a mixed flow pump of capacity 474.9 © 
: gal/s (1,800 L/s). The return activated sludge pumping pit (RASP-pit) has two- 
identical vertical axis, axial flow pumps of Tag 105.5 gal/s (400 L/s) each, 
i Model scale ratios of 1:8.750 and 1:4875 were uantia for DDP-pit and RASP-pit 
| respectively, taking into consideration the availability of standard pipes and 4 
4921 
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‘BENCH ‘BAFFLE WALL 


FIG. 1.—Line of Pits, Broken Line Showing Modifications to Orginal 


Pd Pumpin Pits = 
suitability fam a minimum 1m model size for of at at dierent 
modelling laws. Geometrically similar models were made with plywood and 
transparent material (perspex) was used for walls and pipe of 3. 
diam t to visualize the flow within the pumping pits. tf 

A standard experimental setup (2) was used. It had a pair of antiswirt vanes 
and a constant speed pump. The swirl near the pump intake was measured 

: using a vortometer (3) and it can be shown that tan @ = U/ V, in which $ = swirl — 
angle; U = tip velocity; and V, = mean axial velocity (4). Axial al velocity profiles 
att the intake were measured using a log linear velocity rake. rake. 

4 Experiments | were conducted with models of both pumping pits at flow rates 


g 
be 


corresponding to Froude, “equal- velocity,” critical Reynolds number simula-_ 
_ tions and at intermediate values. Part of the investigation was flow visualization és 
within the models using dyes and tufts, and the other part was determination __ 
_of swirl angle with and without modifications. The acceptable configurations 
; which had the best flow conditions were determined by the usual trial and = “i 
error approach and in the case of the DDP-pit model. The best configuration — Pr 
W was also found by determining the efficiency of the pump used in the model Be 
For the best configuration of pumping pits the velocity profiles were 
4 _ also determined. Different problems were identified with both models, the DDP-pit 7 
ie had air entrainment, surface undulation, and surface and submerged vortices, _ 
while the RASP-pit had high swirl angle and submerged vortices. 


Flow Conditions and Swirl Angle. —In the DDP-pit, introduction of a bench | 
reduced the strength of and a baffle wall eliminated ‘the vortices formed 
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“SWIRL ANGLE IN DEGREES 


‘ 


R = 4.46 * 10°) 


— RATE, 


vortex persisted which v was eliminated by an inlet with a split 
The RASP-pit had a high swirl angle which could have keane’ reduced to a 
_ lower value using a cross vane within the bellmouth, but it was not acceptable a 
= to the user of pumps. Therefore different types of modifications such as sback 
wall, cross vane at the pit invert, etc. were attempted. Finally an inlet cone ea 
with a split vane and a back wall reduced the swirl angle to within an acceptable 
_ limit. The inlet cone also eliminated the submerged vortices between the eee Z 
_ and the pit invert. As the pumps were symmetrically placed, the same values _ 
_ for swirl angle and velocity aw were obtained but the queen was | mirror a. 
In both pumping pits the swirl angle increased with ‘the nen ta ‘number 
? 


| — 


ot 
FIG. Variation 


have very little effect on swirl angle, but ‘substantial difference is seen 
with the RASP-pit. This shows that the cross vane within the bellmouth is — 
very effective in controlling the swirl angle in the pump intake. The swirl angle a 
obtained for the final configuration of models is less than 2°, which was a he 
stipulated requirement of the pump manufacturer. bos 
® Velocity Profiles at Intake.—The design performance of a pump could not > 
be met and would be affected by the type of velocity profiles at the inlet 
of the impeller. The velocity profiles were determined for the best configuration 
of both pumping pits at flow rates corresponding to critical Reynolds number > 
4 


EQUAL 
OR 6.46 108 


that the critical Reynolds number satisfactorily simulates the velocity profiles 
of the prototype. Also, velocity profiles would depend on the size of the intake 
a Efficiency as a Parameter.—The best configuration is usually determined ; 
subjectively by observing the flow conditions. In this investigation for various 
i ‘positions of the baffle wall, the efficiency of the pump used in the ae 


study of the DDP- P-pit was determined at ‘‘equal-velocity’’ simulation. Fig. 4 


Dut Varial sma ip chow: ha he modi ions on a DDP-ni 
j 
d . 3.—Normalized Velocity Protiles at Ury Pump Intake Determined at Ditterent 
‘‘equal-velocity”’ simulation. Fig. 3 shows the normalized velocity distribution 
with respect to the center velocity of the DDP-intake which is almost the same. 


TECHNICAL. NOTES 


shows that the relative efficiency (n/1,) has an optimum position { for the baffle 
‘ - wall, which incidently coincides with the position determined by the trial and 
_ error method. The variation of efficiency of the pump used in the model study 4 
_is mainly caused by the airentrainment inthe inflow. 
‘Cates and Submerged Vortices.—At a flow rate corresponding to Froude ms 


simulation, surface circulation around the suction were aoticed | in both 


those in the DDP-pit became strong intermittent air entraining vortices. Senin 
that when the velocity in the model pit approaches that in the a 


; ». Velocity profile measurement also leads to detection of he 
the Gh The ir the ‘The 


b 


a = 
‘FIG. 4.—Variation of Relative Efficiency of Pump Determined at Different Positions 
of Baffle Wallin Dry Dock Pumping Pit 
ofiles in 


‘ umneget vortex which causes localized distortion of the velocity pr 


~ the Reynolds number, and the critical Reynolds number of 5x 10° gives 

> 
A satisfactory simulation of the velocity profiles at the pump intake. The model 
_ pump efficiency may be used as a parameter to determine the best modification 
- aeeerey. From the model investigation it appears that it is preferable to 

= conduct the study at different flow rates corresponding to Froude “equal- -veloci- 
and critical Reynolds number simulations which would lead to a better 

of the of pump | pit. 


q f Submerged vortices are detectable at Froude simulation using a dye, but J 
a high flow rat 
Tree 
q 
{ 
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URFACE TENSION EFFECT ON PROFILE 


OF A FREE 

7 Nevzat Yildirim’ and Subhash C. Jain,’ M 3 


call A aumber of experimental studies (1,3,4) have been reported on the similitude 
in free-surface vortex formation. The knowledge of these similitude relations 

is of importance in designing intake models to study the prototype  eations 

_ that lead to vortex formation. Based on the experimental results (1,3,5) it has» 

~ been generally accepted that the surface tension effects on the formation of 
_ the free surface vortex are not important in most test systems. In this technical — 

note, the effect of the surface tension | on the surface _ of a free a 


in = z and r directions, respectively, canbe writtenas = | 


Note. —Discussion open until June 1, 1981. To extend the closing date one month, - 
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ASCE. Manuscript was submitted for review for possible publication on July 8, 1980. 
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- Consider a differential element HEFG of a rotating fluid with the a 
i in Fis j j i ial element 


TECHNICAL NOTES 


1 


‘ free surface and the bottom of the element i in 1 the vertical The 
first and second terms on hand = of 2 are the resultants 


ae FIG. 1.—Forces on Element of Rotating 
the pressure forces acting on the tangential faces a and radial faces of the clement, 
- gespectively. The third term is the pressure force due to atmospheric pressure 
i acting on the free surface of the element, and the fourth term is the resultant 
Ss due to surface tension in the radial direction. The term on the right hand — 
side is the product of the mass and the acceleration of the element. For the 
- free vortex, the tangential velocity is assumed to be inversely ‘Proportional t to 
a The elimination of of (p — - Po) from rom Eqs. 1 and 2 and some trigometric manipula- 


sin + = 2 A. 


“which after dividing by « or cos’ sand using the transformation z = C — reduces 


¢ 
&g in which o = suriace tension, p = pressure at the bottom of the element; 
: Po = atmospheric pressure; p = fluid density; g = acceleration due to gravity; _ 
@g _ The various terms in Eq. | are the resultant force due to the surface tension, —— 
il 


‘in which C=a constant; € = the drop in the free surface; =r 2/4n%g 

primes represent derivatives with respect to r. on EN 
a Eq. 4 is an ordinary, nonlinear, differential equation. An analytical solution — 
- of this equation is difficult to obtain. Therefore it was solved numerically wsing 
‘the Runga-Kutta method. For a third order equation, three boundary conditions _ 
are needed. At large r, the curvature of the free surface is small, and the 
oad? 


(ANWAR 1965) 
 Be1.04x107 
&=0,000 Ib/ft - 
Ib/ft 
—-— EXPERIMENT - 
(ANWAR 1965) 
4.271073 #3: 
Ib/ft 


effect of the surface tension can be neglected. Thus the boundary conditions 


at large r rcan be assumed to be t the same as those without the surface tension. 


4 
| 
Fue. 
| 
q 
4 
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at any “point, | r, can be obtained by "differentiating Eq. 5 with respect 


used to solve 


Sa The comparisons of the solution of Eqs. 4 and 6 with the free-vortex surface i 
- profiles ‘measured by Anwar Q) are shown in Fig. 2. There is ‘no effect o of 


O16 


ft? and 0.00427 ft*), 
_ The variation of the normalized difference between the surface drops without 


surface tension (€,) and with surface tension (é, ) for several values of circulation — 
is shown in Fig. 3. The effect of the surface tension for the low values of 


the circulation becomes significant, particularly near the core. The effect all 


“ significant at low values of the circulation pe may give rise to some = 
effects i in model studies on vortex 
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name should be indicated below the title (see Discussions herein as an enampie) together with 
ASCE membership grade (if applicable). 
The discusser’s title, company affiliation, and business address should appear on the first | 
page of the manuscript, along with the Proceedings paper number of the original paper / technical — 
note, the date and name of the Journal in which it appeared, and the original author’s name. 
Note that the discusser’s identification footnote should follow consecutively from the original 
paper / technical note. If the paper/technical note under discussion contained footnote numbers 


_ | and 2, the first Discussion would begin with footnote 3, and subsequent a would 
continue in sequence, 
z Figures supplied by the discusser should be designated by letters, starting with A. A. This also 


_ applies separately to tables and references. In referring to a figure, table, or reference that | 
iE appeared in the original paper/technical note use the same number used in the original. _ sae, " 


It is suggested that potential discussers request a copy of the ASCE Authors’ " Guide to 
_ the Publications of 4 


EC 
q 
7 
4 
7 | 
of manuscript including figures and tables); the editors will delete matter extraneous to the > tn 
a _ subject under discussion. If a Discussion is over two pages long it will be returned for shortening. __ 
_ All Discussions will be reviewed by the editors and the Division’s or Council’s Publications _ in 
“a 
4 
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A VerticaL Rounp Buoyant Jet" 
Closure by Nikolas E. Kotsovinos 


Fine 


vs The writer would like to thank the discussers for their interest in the paper 
and for their valuable comments. The writer thinks that high values of o are 

_ probably caused by experimental errors due to: (1) Stratification of the ambient , J 
y water; and (2) boundary or free surface effects. The writer’s opinion is that v7 
- accurate measurements in concentration, velocity and their turbulence charac- - 
“teristics are necessary to. resolve questions about the details of the buoyant 
jet behavior in the plume region. However, it is worthy to point out that Eg. 

6 (or Eq. 11) is based on various assumptions and simplifications and therefore 


even the most accurate experimental data would not have to follow the curve | 


Errata.—The corrections should be made to the paper: 
Page 796, Eq. 4: Should read =aTy (x)i instead of =a Ty(x) é 
Two-Dimensiona BusBLe Puumes® 
_ The writer is grateful to W. Hall C. Maxwell for drawing his attention to ot 
: “errors i in the expressions for the Reynold’s Number, R, and the Weber Number, 7 z 


_ Errata.—The following. corrections should be made in the original paper: . 


Page 142, Table 1, Col. 6: Should read “300, 49, 10. 3, 8.2, and 21” instead 


a Page 146, paragraph 1 line 5: Should read “‘R = B,/gv is a form of Reynolds 7 


*Prof., School of Engrg., Democritus Univ. of Thrace, Xanthi, bute 


“February, 1979, by David L. Wilkinson (Proc. Paper 14385). aie 
* Visiting Lect. in ‘Civ. Engrg, School of Engrg., ‘Univ. of 


number; and W = op’? /og i is a form of Weber nu mb 
BY? /evi is a form of Reynolds number; and W = pB2/og 


_ (These changes do not affect other calculations or conclusion in the paper. 


Hy DRAULIC RESISTANCE OF ARTIFICIAL Strip RouGuness* ey 

bas 

by Donald w. Knight® and J. Alasdair Macdonald” 

of 

4 The writers are grateful to Eric Valentine, R. re and J. J. Garcia, | 


- for their comments. . The topic of nonuniform roughness is still clearly © one e of 


ns 
‘the shape of the roughness elements. In the writers’ " experiments, the roughness 


elements were carefully machined to a square cross section such that b “4 
_ k = 3 mm, ie., b/k = 1, in which 6 = width of element. The additional 
* presented by Valentine, which confirm the general usefulness of Eq. 11 q 

have b/k = 4 and A/k = 8. The data presented by Fuentes and Garcia have — 

(b/k = 0.52, E: 0, and 10.0 for varying \/k. In both sets, the values of the 

coefficient, C,, differ somewhat from those suggested by the writers. Evidence _ 

ee experimental work in pipes with different transverse two-dimensional ribs 
y (26) strongly suggests that the parameter b/k has an effect upon the roughness © 
by as much as 25% or more. Generally the friction factor increases with decreasing _ 

_ b/k for given \/k and h/k values. The variation of C, with \/k should therefore | 


reflect this as pointed out by Fuentes and Garcia. The writers’ curve is of 
course only valid for b/k = 1. However, the trend indicated in Fig. 11, with 
‘the discussers’ additional data, does not appear to be the 
data points for b/k = 0.52 are below those for b/k = 1, the — points for 


pao = 10 appear to lie on the s same curve as as those for bike . The reasons 
et results 


Sayre and Albertson is subject to a certain amount of error, but it is — 
that all the discrepancy may be attributed to this cause. Until more data are 
- published i in a similar format it is unlikely that this problem will be resolved. ca 
The comments about the inadequacies of current sidewall correction n procedures ‘ 
are thoroughly endorsed by the writers. The boundary shear results undertaken | 
by the writers (25) clearly demonstrate some shortcomings. The first writer 
is continuing his studies in nonuniform roughness by examining wyor-ee 
transfer in symmetric and asymmetric sections using a wind tunnel and a nell “= 
The boundary will be differentially ‘Toughened and it is hoped that _ this will 
shed some light on these basic correction procedures. It is a source of : astonishment 
“June, 1979, by Donald W. Knight and J. Alasdair Macdonald (Proc. Paper 14635). 
*Leet., Civ. Engrg., Univ. of Birmingham, P.O. Box 363, Birmingham B15 2TT, England. 7 
Asst. Engr., Soar Div., Severn Trent Water Authority, Leicester, England. 


= 
instead of “‘R = 
is a form of Weber 
i ] 
i 
| 
| 


- that we know so little abo about shear distribution i in asia cross sections s considering — 
the importance of open channel flow to so many problems. | ‘The writers are 


grateful to have Ref. 22 drawn to their attention. As regards viscous effects, 

* = from presenting the data in the form of Figs. 2 and 3, no further analysis 

was undertaken. Until a more rational sidewall correction procedure is formulated 
= seems to be little point i in examining this issue due to the complex i interaction =) 


“Open Channel Flow with ‘Varying Bed 
" Roughness, ft Journal of the Hydraulics Division, ASCE, Vol. 105, No. HY9, Proc. 
Paper 14639, Sept.. 1979. pp. 1167-1183. 
26. “Losses Caused by Friction in | Straight Pipes with Systematic Roughness Elements,”’ am 
Data Item No. 79014, Engineering Sciences Data Unit, London, England, 
ASCE, of My Val. 106, Prot Peper 


Aswar, G., Flow Vertes,” Jounal af Resverch, Vel 


_ The writers thank Hooke for emphasizing the role of sediment transport rate 


_ over the entire bend in shaping the bed topography. The contention that the 
a shape. of any stable channel i is a consequence ¢ of the balancing of th the prevailing 
shear stress and the ‘sediment transport rate is correct. snail asad 


_ The important fact, however, is that for flow through | bends, secondary 


circulation is an intrinsic feature and is accompanied by and concurrent with 
asymmetry of velocity and shear distribution. Interestingly for flow through 
rigid boundary laboratory channel bends (3), high shear-stress zones along the 
inner bank : at entry shifting | to the o outer r bank a at exit occur r with the development — 


7 bank, the suggestion that the velocity gradient (and consequently bed shear- aco) 

is higher or lower consequent to smaller or greater « depth of flow is not valid. — 

“= _ The locus of the bed shear maxima, however, is associated with the path 

“of of the longstream high velocity filament which has a characteristic swing from 
the inner bank to the outer bank 3, 11). This has also been observed for highly ; 
curvilinear and accelerating flow through end-dump constrictions (12,13). — ee 

_ The writers agree with the discusser that more detailed studies to assess 

_ the role of the secondary flow in a bend in shaping the local channel geometry _ 

¥ (or not?) must be undertaken. In this context it is perhaps worth noting that 


“October, 1979, by Mamdouh A. Nouh and Ronald D. Townsend (Proc. Paper 14898). _ - 
_ “Post-doctoral Student, Dept. of Civ. Engrg., Univ. of Calgary, Calgary, Alberta,Canada. 
of Engrg., Univ. of ‘Ottawa, Ontario, Canada KIN 9B4. 
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Il. Bathurst, J. C., Thorne, C. R., and Hey, R. D., “Secondary Flow and Shear Stress 
at River Bends,” * Journal of the Hydraulics Division, , ASCE, Vol. ™ , No. HY10, — 
Proc. Paper 14906, Oct., 1979, pp. 1277-1295, 


(12. Das, B. P., ‘Stability of Rockfill in End- River “Closures,” of the 


Division, ASCE, Vol. 98, No. HY11, Proc. Paper 9381, Nov., 1972, 


13. Das, B. P., ‘Bed Scour at End-Dump Channel Constrictions,’’ Journal of the Hydraulics _ 


ee ASCE, Vol. 99, No. HY12, Proc. Paper 10228, Dec., 1973, pp. —— 2291. 
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Cavcuarion OF STRONGLY 
an 
Discussion by Habib O. Anwar* and James G. Rodger 


when the flow i is isotropic for which the coefficient C,, C,, C,, 0, anda, remains — 
universally constant, otherwise the model must be * ‘tuned”’ by adjusting the 
_ numerical values of the constants, which, usually differ from experimental values. 

The writers have, however, evaluated the constant, from measurements 


data obtained in the Great Ouse Estuary (Ref. 36, United Kingdom). The global 
a Reynolds numbers for these flows were 9x 10° and 4 x 10°, respectively. 
, It was found that the laboratory C,, value, depending | on y, varied between 
0.008 ar and 0. 08 (C, = 0. 008, near the e free surface) and the field va value o of Cc, 4 
= measured at y= : 1.14 m varied between 0.08 and 0.012. The energy spectra 
showed the existence of a -—5/3 distribution (35, 36), but further examination 
_ Showed that the flows were anisotropic (36). The C, value in a homogeneous © 
_ isotropic turbulent shear flow (38) is about 0.11; al a turbulent boundary layer 
_ with zero pressure gradient (39), depending on y/, 5 being the boundary layer 
thickness » C, was 0. 075 and 0.059 for y/8 = 0. 1 and 0.9, Tespectively. The 
“4 Cc, value used by the authors was 0.09. These results show clearly that the | 
‘constant, C,,, depending on the flow, can vary considerably, = 


‘ag 


_ The K-e model, although successful in describing isotropic turbulence with 


October, 1979, by Michael A. Leschziner and Wolfgang Rodi (Proc. Paper 14927). 
Principal Scientific Officer, Dept. of the Environment, Hydr. Research Station, 
*Sr. Scientific Officer, Dept. of the — Research 
Oxon OX10 8BA, England. ¥ 
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recent analytical approaches to the problem (1,8) employ the ra 

+ 
7 
4 

= 

tal 
results and those of their numerical K-e model. Although it was not mentioned 4 

f 

i 


DISCUSSION 


universal constants, ‘could not be applied, as without some 
tuning to more general flow problems. Therefore, it is not obvious that the 
results from a K-e model compared with aoe from other simpler ee 


computational effort required to solve the K and € equations. ‘were 
_ The authors in Eq. 4 have assumed, by analogy with laminar flow, that the : 
_ Reynolds stress is proportional to the rate of strain. Prandtl (41), on the other | 
hand, postulated by analogy with Rayleigh’s on the stability of 
carved flow that the Reynolds stress, being positive, is proportional to the i 

vorticity (37); what changes would the authors expect had Prandtl’s . 2 


iy Fcuree be tn 


3 35. Anwar, H. O., and Atkins, R., “‘Turbulence Measurements in a a Simulated Tidal Flow,” — 
mel ASCE, Journal of the Hydraulics Division, Vol. 106, No. HY8, Proc. Paper 15609, 
36. Anwar, H. O., ‘‘The Study of the Turbulent Structure in a Tidal Flow,” in press. 
37. Anwar, H. O., ‘‘Turbulent Flow ina Vortex,” Journal of Hydraulic Research, Vol. 
7, No. 1, 1969, pp. 1-29. 
_ 38. Champagne, F. H., Herris, V.G., and Corrsin, s.,“ Experiments on Nearly Turbulent 
Flow,”’ F. Fluid Mech., Vol. 41, Part 1, 1970, pp. 81-139 (see also ‘‘Turbulence,”’ — 
Bradshaw, ed., section 6, Springer-Verlag, Berlin, Germany, | 
39. Klebanoff, P. S.., ‘*Characteristics of Turbulence in a Boundary Layer with Zero | 
" Pressure Gradient,’ ,”’ Report No. 3178, National Advisory Committee for Aeronautics, 
i ‘M0. Odd, N. V. M., and Rodger, J. G., ‘*Vertical Mixing in Stratified Tidal Flows,” 
ASCE, Journal of the Hydraulics Division, Vol. 104, No. HY3, Proc. Paper 13599, : 
Prandtl, L. .. “Einfluss stabilisierender Krafte auf die Turbulenz,’’ Gesammelte Abhand- 
Zweiter Teil, Springer- Verlag, Bertin, Germany, 1961, pp. 


#) Recent investigations i, 32) have shown that the dispersion process in natural 
rivers cannot be accurately predicted by the well-known Fickian dispersion theory. 
The concentration-time distribution curve is generally found to be skewed -* 
a sharper rise and more prolonged fall in the initial period. The author has 
d proposed a dispersion model with two coefficients to reduce the discrepancies _ 
between field data and the Fickian theory. Application of this model to natural — 


"January, 1980, by Spyridon Beltaos (Proc. Paper 4 

Research Asst., Iowa Inst. of Hydr. Research, The Univ. of Iowa, lowa | City, | lowa ] 
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_ streams shows good agreement between the data and the model ~The pencticnd 
; advantage of the method using several C versus ¢ curves at fixed cross sections — ‘i 
_ to determine AT values with a higher accuracy rather than &? should be — 
_ emphasized. However, the writer would like to discuss a few points which 
do not seem to be clear in the development of the model and presentation 
_ The value of a is given as 5.55 and assumed to be a universal constant. 
_ However, it is not explained how this constant quantity was obtained. ‘The 
coefficient a can be calculated either numerically or using the definition of 5 
 @ as given in Eq. 30. For a numerical solution, Eq. 22 can be rearranged for 


to ¢, and ay - the C versus ¢ curve, respectively. Since AT = T, — T, and 

AT= (a) » one can am calculate a. According to the writer’s calculation 

TABLE 1.—Application of Model to Fischer's Flume Data 


3100 . 10.2 | 30.0 23.1 0.800 3.734 3.559 


“Taylor length, L,, is calculated from Eq.7. 


was found th that for 0 0.0070 00700 < < 0. 01400, a= 5, 55, a increases slightly 


this Tange. Use of t to ) define a requires known values of 


a Fischer’s experiments (15) conducted i in laboratory flumes should nearly satisfy 


the condition of one-dimensionality required for the application of the author’ 7. a 
“model. Values of a were computed using Eq. 30 for series 2800, 2900, fT 7 
3100 of Fischer’ s experiments, and are presented in Table 1. The a wae 
obtained are seen to be much smaller than 5.55, and to increase with | increasing — 


‘x. It ible that Id Il h the val f ’ 
is e that a ‘wou vould eventual reach the of 5. 


444 
X,, | Xz |G, |G, [ters |AT,, |AT,, | 
— in | in | in per in in | in | 
ischer's| me- | me- | sec- | sec- | sec- | sec- | sec- | me- _|me-| — a ; 
data | ters | ters ondsjonds* ond |onds jonds| ters | a, tes | B 
| (1) (2) | (3) | (4) 7) | | | | 13) 
2800 | 16.07) 26.11} 369 651 0.251 12.3 | — | 30.9 0.410 2.954 0.300 
2900 |16.08| 26.12} 111 179 |0.454| 6.0 | 18.0 | 28.1 |0.324 |1.810| <0.0 
25800 
= 
i accurately, determination of accurate &? from experimental data is difficult 
ine since &? is very sensitive to the tail region of C versus ¢ curves. It is possibly bs ; 
q 


le. the aforementioned expesinente, AT is not to x; therefore, 

_ it is unlikely that Fischer’s measurements were conducted in the linear | Tange. 

a Therefore, values of B and L were e calculated from Eq. | 15 using the data: measured 2. 
at two locations, X,, and X,, and the Taylor length, L, was calculated using __ 
Eq. 7. The results of these computations are also included in Table 1. The 
following comments can be made: vemmegative ordinsiss oF the 

1. The values of B are much larger than those suggested by the author. _ ; 
2. The length L is less than one third of re that was found by the oe 
* 3. Since X, and X, are greater than 3L in Fischer’s experiments considered 
herein, the corresponding C versus f curve should be in the Fickian range accord- bn “4 
ing to the author’s criterion. On the other hand, “ee and X, are still significantly 
Smaller | than ' and, thus, they are within the intermediate range. This is 


17.50 m,v=0444 ome 


FIG. 10.— omparison of Fischer’ s Data with Proposed Model P 


ed 


also confirmed by the fact that the Cc versus ¢ curves for the aforemention 
riments are very skewed (more s skewed at X= than at =X, ). 


- _ It was then assumed ‘that the low values for a as calculated above were 


due to inaccurate determination of ?. The method proposed by the author - 
for the determination of B and L ‘from Eq. 31 was | thus attempted, using the 
3 AT values at the t two. ‘measuring, stations in Series 2900 and 3100, and a = 


= 


3.39. However, for these two particular s series no finite value for L can be 
obtained. If indeed the L values were so large as to be much greater than 
the flume length, Fischer’s experiments then would have been all conducted = 
in the linear range, which is in contradiction with the observation that 4T a 7 
a not proportional to X. Furthermore, when the data of series 3100 at X¥ = X, fj 
and X = X, are compared with the prediction of Eq. 25 with a = 5.55 and 
: obtained heen’ Eq. 27, it can be seen from Fig. 10 that the agreement is 


very poor at X = X,, although somewhat betterat X= X,, 


| 
i 
| 
4 
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data and the author’ s model prediction can be offered since on the one hand, 

§ flaw was discovered in the author’s analysis, and on the the other hand, Fischer’s i. 
data are widely accepted as being reliable. The author’s comments on the possible P 
source for the disagreement discussed here would be extremely welcomed. a 


31. McQuivey, R. S., and Keefer, T.N., “Convective Model of Longitudinal Dispersion,” bs 
Journal of the Hydraulics. Division, ASCE, Vol. 102, No. HY10, Proc. Paper 12478, 
32. McQuivey, R. S., and Keefer, T. N., ‘ ‘Dispersion- Mississippi River Below — 5 
Rouge, La.,’” Journal of the Hydraulics Division, ASCE, Vol. 102, No. HY10, il : 


4 Discussed by Dale D. Meredith,° M. 
= Timothy Byrnes’ 


authors. have presented an interesting paper on the of a 
composite unit hydrograph derived simultaneously from several runoff hydro- 
graphs. There are, however, items that need further explanation. 
_ The authors used Eq. 2 to compute the number of observed direct runoff 
- bydrogaph ordinates from each storm to use in the analysis. However, the _ 
number of unit hydrograph ordinates, M, must be determined before Eq. 2 
= be used. From the data it appears that M could be 14, 18, and 16 for 
storms A, B, and C, respectively. The authors chose M = 16 without explaining — 
how or why this value was chosen, 
relatively easy to formulate the problem such ‘that the optimal num 
— of unit hydrograph ordinates can be obtained simultaneously with the values 
of the ordinates. Choose M somewhat larger than the largest value indicated 
by the data, e.g., here choose M = 20. Then Eq. 2 can be used to determine 
_ the number of direct runoff hydrograph ordinates required. If the number of 
required ordinates is greater than the number of observed ordinates then the 
additional ordinates are assumed to” be zero, €.g., here the last six runoff 
hydrograph ordinates for storm A would be 0. This formulation increases the 
; number of decision variables and the number of constraints. However, Eqs. 
16 and 17 are still valid. The number of nonzero unit hydrograph ordinates 
in the linear programming solution will be the optimal number of unit nycrograpa 


wu * January, 1980, by Larry W. Mays and Lynn Coles (Proc. Paper 15129), 

*Prof., Dept. of Civ. Engrg., Parker Engrg. ‘State Univ. of New Now York 

Water Resources Planner, U. Ss. a — of Engrs., Buffalo D Dist. ‘Buffalo N.Y 


_ Using the formulation presented here the optimal value of M is 17 for storms 7 7 
A, B, and C. The : resulting unit hydrograph i is shown in Col. 2 of Table 6. —— 
.* _ The authors stated that: “‘least squares methods cannot guarantee nonnegative — re 
ordinates of the derived unit hydrographs. The linear programming method has 
nonnegativity constraints that eliminate any occurrence negative ordinates.” 
Least squares methods that guarantee nonnegative ordinates of the derived unit 
__ hydrographs can be used. The procedure is simply to formulate the problem — 
asa nonlinear programming problem with nonnegative constraints. Statisticians 4 


problem except that ‘the objective function is now nonlinear. 
Nonlinear programming are now available at at most commercial com- 


ra Collins, M. As ase 


ALB, C, in inches per hour* a 


Sumof | $Sumof least 

differences differences | squares 


00279 


a 0.01580 
0.01137 


0,00083 


was ‘stated as a least 

with M seen as indicated in the foregoing. The solution was obtained using 

sg available nonlinear programming algorithm (23) and is shown in Col. 3 of 

Table 6. Again the optimal value of M is 17. Tina 

‘The nonlinear model function was then stated as a 

weighted least squares 


least squares ee 


January 1981 


the mean of the ordinates for the ith direct runoff hydrograph. 

— The optimal value of M for this model is 15 and the > resulting” unit hydrograph ‘ 
is shown in Col. 4 of Table 6. The author’s unit hydrograph is shown in Col. : 

The authors state that: ‘‘An important factor in deriving unit hydrographs 


to define the “peak and time to peak.’’ In order to evaluate the 


7. of Runoff Hydrographs for Various Unit Hydrographs, ir in inches 
per hour” 


te 


Sum of differences (M = 16) sale 08890 = | 
Sum of differences (M = 17) 0.08888 

‘Least squares (M = 17) 09025 


4 Observed 0.02560 


results from the programming models the four - unit hydrographs in. in Table 6 
were then used to compute direct runoff hydrographs for storms A, B, and 
C. The peaks for these computed hydrographs are shown in Table 7 along 
with the observed peaks for comparison. The linear programming models (sum 
of differences) resulted in a peak greater and two time periods later than the — 
observed peak for storm C. The nonlinear programming models (least squares — 
; and weighted least squares) resulted in a more accurate estimate of the <_ 
for storm C but in slightly larger peaks for storms Aand 
_ Therefore, from this preliminary analysis it appears that the nonlinear program- 
ming models may be more accurate than the linear programming m models because 
_-peak weighting can be preformed. In addition the models appear ir to be i 
accurate when the number of unit hydrographs ordinates are determined by 


Pesther study should be done to evaluate these programming models. It would 


ond total runoff hydrograph ordinates. Ref. 24 « can be used to determine a an 
optimal unit hydrograph for a single storm event with a weighted least squares 


= 
148 
Mind > d 
in which d.. = a weighting factor used by the U.S. Army Corps of Engineers _ 
&g 
(4) 
9.09228" 
0.09207 
0.090815 
0.090465 
i 
d 


Ne. 320, Version 3, Northwestern University, Evanston, Ill., 1976. ia @anagae 
, 24. HEC-I, Flood Hydrograph Package, User’s Manual, Corps. of | of Engineers, Hydrologic — 
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EXTENDED Ser OF CoMPONENTS IN NETWORKS© 


these will by Michael A. Collins, M.ASCE 


as check valves and pressure ssoacins valves, and providing a systematic 


_ description of the proper network topology to be used when such elements 


are presented ina network, 


‘The author, in his conclusion, raises the important is issue ¢ of solution uniqueness, 


that a unique solution exists (9). The seouite of Ref. 9 can also be teadily — 
_ extended to networks in which the only control elements are simple check © 
valves since, as has been demonstrated (10, 11), a network containing both — 
‘ unidirectional elements, such as check valves, and bidirectional elements, such _ ty 
; 2 as pipes with simple frictional losses, can be modeled by an equivalent network 
consisting of solely unidirectional elements and the results of Ref. 9 apply to | 
_ La Other types of control elements, which for modeling certain modes of operation 
_ require a changing of the network topology beyond the simple replacement 


of bidirectional elements by two equivalent but oppositely directed unidirectional © 


; } elements, indeed appear to offer the opportunity to yield nonunique steady 


state solutions (8). This potential nonuniqueness arises, of course, because when 
the network topology is altered one essentially has another network problem, - 
and there appears to be no reason yet proposed why one solution to a steady-state. 
network should rule out a solution to another different, = similar, steady-state 
3 In networks wherein the control elements are. pressure ond we valves, 
when a unique solution does exist, the criteria as described by the author to 
identify the proper solution do not appear to this writer to be Cones. ~The 
remainder of this discussion addresses this apparent deficiency. __ | ae eo} 


January, 1980, by M.  Chandrashekar (Proc. Paper 15130). 
Prof. of Civ. Engrg., Civ. and Mech. Engrg. Dept., School of Engrg. and Applied 


Sei. , Southern Methodist Univ., Dallas, Tex. 75275. 
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this issue. In networks without control elements and in which the head loss- = 
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criteria are given, of which only are concern (1) If 
: < H,,, and H, < H,, . then the pressure reducing valve is inactive; and (2) 
if H,> H,,.. ont H, < H,, then the pressure reducing valve is active and the Pt 

admittance is changed from y, to y,. Because one does not necessarily know 

if a unique solution does exist, one should consider all the possibilities of 

alternative operating modes of the pressure reducing valves. One might argue — 

that if one initially finds that H,< H,.,, one needs to go no further. This 


Seems ‘Plausible, but if one admits to the of roe 


one does accept the premise that if the first criteria (H,< H aa is satisfied 
one has the correct —. then the author’s criteria would imply that if 
it is found that H, > H,,., then the proper solution corresponds to an active 


= = operation for the pressure reducing valve. But what happens if there 4 


12.—Example Network with Two Pressure Reducing Valves 


in the author's ‘criteria, consider the four node network of 
= 12 consisting of three pipes, of which two have pressure reducing valves. : a 
Nodes 1, 2, and 4 are reservoirs with fixed heads of 100, 100, and 0 units 
of head, respectively. For simplicity and without loss of generality, flows, heads, - 
and admittances will be described without reference to any specific set of units re 
with the understanding a consistent set of units is being used. Pipe ane 


_ umber 3 between nodes 3 and 4 is assumed to be described by Eq. 3 with 
_ The pressure reducing valves in elements 1 and 2 prevent arseabid pressure 
at node 3. The control set head in the valve in pipe 1 is H,., = Aes ass 
90.7 while that for the pressure reducing valve in pipe 2 is H.. = Hy. = 
90. Note that the different set heads for the two valves do not make the problem — 
unrealistic since the network of Fig. 12 could be considered to be a simplification — ] 
“of some more complex network with more pipes and nodes between the = 


| 
—— 
4 


y, = 0.200 and y .= 0. 250. Both pressure reducing valves are located midway 
along | the pipe on which they occur, ie., £,/L, =L »/L, = 0.5, in which 
i” ‘iL, = the length of pipe k; and Z, = the distance along pipe k from its downstream - 
end to the pressure reducing valve [see Fig. 4(a)]. Since the author implies 
that admittance for an active pressure reducing valve is modeled wg a simple 
4 Pipe (se (see Eq. 11) in which the effective admittance 7, is based on L,. Pe ‘it follows 
that = y,(L,/L,)°™*. Thus = 0.291 and 7, = 0.363. 
® A tentative solution (in the sense that it is ‘not known a priori whether the 
= will satisfy the author’s criteria) can be obtained from a contianty 


z reducing valves and node 3. When these pressure sedecing valves are inactive, 


equation for node 3 expressed in terms of head losses: 


in which the heads H * and H* have been introduced to denote the fact that 
these heads will have different valves according to the operation of the pressure 


valves. The variable, yt, is merely Ve or y, depending « on 
the pressure reducing valve i is ina inactive (or active, 


1.—Both reducing valves are inactive; thus Hf = H, 

and H} = 100. 


2 —The valve in pipe 1 is active but the other valve 


inactive; thus Hf = H,.,,, 
Case 3.—The valve in is active but other 
inactive; thus H? = =! = 100 and H H 


= 


tentative for case is from E Eq. using a simple 


81.7840 
81.7840 
81.7840 
81.6872 
81.6872 
81.6872 
81.2095} 1 | 90. 9749 | 0.200 | — | 0.2908 
81.2095 1740 0.363 0.3386. 
81.2095 a 
81.2518 
81.2518 
81.2518 


ent 


4 
4 
TABLE 7.—Summary of Alternative Tentative Solutions for Example Network of 
Fi 
: 
| 
| 
= 
| “If the y, column is blank, pressure reducing valve is activated; if the 7, column is 
blank, pressure reducing valve is not active. 


ah 


M ariable, H,, is the computed head at node 3 while ¢ qk in| Col. 5 is the discharge 
computed from this head the ac admittances- and the “upstream (or 
downstream in the case of pipe 3) heads appropriate to the assumed operating — 
: — of the pressure reducing valve. The variable H} is the head at the 
downstream side of the pressure reducing valve. In the case of inactive valve 
operation it is merely the average of H, and the reservoir head since the valve 
a. midway between the reservoir and node 3. In the. case of active operation 
ofa a valve H?= The variable yi is the apparent admittance that mt aang 
; ee upstream of the pressure reducing valves in view of the computed diachasge 
a the tabulated valve of H*. The variable, H,, on the other hand, is the 
: A, = at the upstream side of the pressure reducing valve as computed from a 
‘ the discharge, the upstream reservoir head, and admittance y, adjusted for 
the shorter length of pipe rd the valve; these latter admittances 


_ It is seen from Table 7 that « case 1, ae which both valves are inactive, yields al 
heads at the pressure reducing valve which exceed H,,.,. Thus, by the author’s _ 

- criteria, this case is not the proper solution. But examination of the other cases — 2 
Shows, as well, that H . does not simultaneously satisfy H, < H,,., for both — 


e author’s criteria as given would ‘imply no solution 


valves i in any case. Thus, the auth 


-* would appear another physical condition, not apparently explicitly cited by 
author, is necessary. Physically, a pressure reducing valve becomes active 
| £ ‘ by increasing the effective head loss coefficient across the valve, i.e., the valve 
a admittance decreases. Thus, for a physically acceptable solution, it would be Jaa 
necessary that for each active pressure reducing valve y; be less than the yy; 
FE for the inactive operation of the valve. Table 7 shows that case 3 | satisfies ; 
this criterion and is thus apparently the correct solution. 
This writer would suggest author’ s second criteria be replaced 
If criterion number | is not satisfied, asi a solution in which for each 
active pressure reducing valve the apparent upstream admittance is less 
than the corresponding inactive value of the mer admittance a an 


acceptable solution, provided H, “9 
Of course, these criteria still do not indicate which nihil is the proper solution — 


Collins, M. in Pipe Network Techniques,” 
No. TES, Proc. Paper 15675, ‘Sept. 1980, 
9. Collins, M. A., Cooper, L., and Kennington, J. L., “Multiple Saw ey Points in 
Complex Pump Networks,’’ Journal of the Hydraulics” Division, ASCE, Vol. 105, 
No. HY3, Proc. Paper 14452, Mar., 1979, pp. 229-244. 
(10. Collins, M. A., Cooper, L., Helgason, R. V., and Kennington, J. L., “Solution of — 
Large Scale Pipe Networks by Improved Mathematical Approaches,”’ Technical Report — 
IEOR 77016-WR 77001, Southern Methodist University, Dallas, Texas, Aug., 1978. 
. Collins, M. A., Cooper, L., Helgason, R., Kennington, J., and LeBlanc, L., ‘Solving . 
the Pipe Network Analysis Problem Using 
Science, Vol. 24, No. 7, Mar., 1978, pp. 747-760. 
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Hyprautics tn Estuarine CHANNELS" 
‘The author has presented an interesting and novel contribution to tidal 


hydraulics. The writer would like to raise several points. that - his opinion — 


_ The author began his paper with the statement that ‘Intuition tells us that 


- alternating tidal flow in an estuarine channel has similarities to alternating current © 
g in an electrical circuit.’ ’ As intuition is not an objective research tool, 
analogy can be proved by comparing the mathematical describing 
certain electrical and hydraulic systems (10,11), on | 

x Furthermore, the dimensional analogates for DC flow (Eqs. 3 = 4), suggested 
b by the author, are confirmed in the analysis of similarities between the mathematics — 
of basic electrical and hydraulic laws. As has been shown i in Ref. 10, sp pal 

_ first law generally corresponds to the continuity equation (conservation of mass) 
in open channel hydraulics, and the Kirchhoff’s second law corresponds to 
the momentum equation for steady, and uniform, open channel flow. 

_ The author introduced a lumped, electrical analogy model to describe t 7 
‘ge " velocity-head relationship in each of two gauged cross sections of the —— 

_ Waterway. The physical phenomena occurring in each of these « cross uli 
"were represented by means of ‘an electrical resonant circuit f orced bya a sinusoidal ee 

However it is difficult to why this particular series circuit strectere 

_ wasa priori assumed to model the hydraulic phenomena in question. The lumped 

- circuit, suggested by the author, cannot account for the transformation of a 


flow along the reach. ’ The applicability of this approach is limited to modeling 
the -amplitude- phase relationship between ve velocity and head of tidal flow at 
- agiven section. Indeed, even at a single section, the proposed analogy is doubtful 
i between measured tidal head and an ideal voltage source, independent of 
i conditions in the electrical system. The sinusoidal forcing functions (analogous ] 
_to harmonic sources of electricity) appearing in the system of flow in the analyzed 
ay are physically ; located at its two boundaries. Velocities and i 
‘- of flow in intermediate cross sections are the distributed system responses to a 
_ the ocean inlet conditions and to the tributary inflows. The process of flow 
transformation along the waterway can be modeled by a distributed electrical 4 
_ Based « on the harmonic linearization of the resistance ce term (procedure identical 2 
to the author’ 's development. of the * “equivalent damping coefficient’ Ippen 
(9) modeled the tidal flow in estuaries by means of the telegraphers Ss 


“February, 1980, by Robert M. Snyder (Proc. Paper 15202), 
-? Alexander von Humboldt Foundation Fellow, Inst. Wasserbau III, Univ. of ie, 7 
Federal Republic of Germany; also Research Hydro., Inst. of Geophysics, Polish rye eel 


of Sci., Warszawa, Poland. 


— 
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| 
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line, and the linear parabolic flow model (diffusion is mathematically 

amplification of the tidal head, in ‘the area measurements: were 

re can probably be explained by variations in the channel geometry (ie e., * "] 
aa of channel width), and accounted for in a distributed electrical model. 
“Moreover, the meaning of the new term “‘DC flow,”’ and its relation to the — 
well established terms like ‘ “steady” and ‘‘uniform’’ have not been precisely _ 


As seen nin Fig. 2, the agreement of calculations with the recorded data, + 


for the case of 57 acres, is questionable. In fact, the tidal current | leads the ol 


tidal head in the first se ment of sinusoid presented in Fig. 2 and so contradicts 


-Ippen, A. T., ed., Estuary and Coastline e Hydrodynamics, _McGraw-H Hill Bo Book 


esis “Inc., New York, N. Y., 1966, Pp. 505- 507. th 
10. Kundzewicz, Z., ‘‘Elektrische Analogien zur Modellierung hydrologischer Systeme,” 
submitted for publication in Archiv fur Elektrotechnik, 
. Shen, J., ‘Use of Analog Models in the Analysis of Flood Runof: on, S. — 
Pro, — Paper 506-A, U.S. Geolo ical Survey, 1965. 


Ocean Outratt Diwution: EFFECTS OF 


if Discussion by Henry H. C. Lin,” . ASCE 


poet ieee the dynamic effect of prevailing currents for diffuser design. 
This concept was utilized in the design of an ocean outfall em, for th he 
and county of San Francisco, California. ad 4 q 
_ The San Francisco Ocean outfall system, “currently being designed to have 

a peak capacity of 670 mgd (2.55 x 10° L/day), will include three separate — 
outfalls. The dry weather pipe will discharge all flows up to a capacity of | 
190 mgd (722 x 10° L/day), a distance of approx 4 miles (6.4 km) a fowe 
at a depth e 80 ft t (24. 4 m). Two wet weather pipes will discharge all on 


May, 1980, by Philip J. W. Roberts (Proc. Paper 15429). 


CH2M Hill, Inc., , 450 mad Sen Calif. 941 il. 


‘introduced in electrical engineering for a transmission line. The electrical-hydraulic 
i analogies go even further; as shown by the writer in Ref. 10, the frictionless or 
’ 
4 
in 
% 
4 
J 
<7 


a nominal diffuser length of 3, 060 ft 


design to determine diffuser length and location. The first model, “‘finite line 
plume”’ (5), as developed by the author, considers the discharge to be a finite — 
TABLE —Dilution of Plume | Model Test with Ma 


g Two mathematical models were used in the San Francisco ocean outfall system 


OF, age 
Density ine ini- | Av dilu- er 
profile tion 


ent to determine average dilution. 


red minimum dilution — predicted average initial dilution) /measured 


__ Note: The initial dilution predicted by Roberts’ model is not as sensitive as Koh’s s 
- model for a change of the port configuration i in a fixed diffuser length. as 


ocean current in both speed and direction, but ignores ‘the momentum of the ; 
discharge and individuality of the multiple ports. It is equivalent to using the . 4 
_— plume formula (3) at no current condition, which does not take into me 


. account } the age of discharge, jet momentum, and merging of plumes. The 


a of interference from the multiple ports, as well as the effect of | the 
_ effluent field’s finite thickness. The ambient ocean is assumed motionless during P 


the buoyant rise of the — Yow 


Koh’s Model | Roberts Model 
tien | Error, 
gallons asa 
(1) | (2) | @& | @ | | T (10) 
1,300h) 0 | 32 | 34 27 15.5 | 31 
1,300 0.25 41 | 60 | — | 4 
DwW-l |July |e 75 | 270 | 9% 7.0 
DW-2 107 | Uniform | ‘192 | 220 | 210 | -95 183 | 45 
DW3 107 July 0.25 280 | 300 126 | 55.0 
bw-4| 107 [Uniform | 0.25 | 400 | * | — | 
DW-5 | 190 |July oO 80 64) 60) 8 | -23.5 
DW-5S | 190 |July 90 64 | 23.0 84 
pw | 190 |Uniform| | 139 | 170 | 156 | -120 125 | 100 


ae the project design, a physical model to study plume behavior was 
: conducted at the California Institute of Technology (1,6). This model study was 
_ used to investigate the near-field dilution, height of rise, and other plume behavior __ 
_ expected from alternative diffuser designs under various flow aos ocean condi- 
tions, and to verify and to calibrate the mathematical models. 
physical model tests determined the initial dilution to be i 


=a 


| 
DESIGN PROFILE 


: 
iy a2 2B 


« 


AIG. 14 Donat Profiles Simulated in Laboratory Tests 

the horizontally spreading Portions of the effluent plumes from the diffuser 


‘The ocean conditions ‘included. typical and worst -density- y-stratified cases as 
_ measured at the proposed outfall diffuser sites during the design period. Ambient 5 
ocean cross currents of 0 knots (0 m/sec) and 0.25 knots (0.12 m/s) were 

modeled. The wet vet weather diffuser was tes tested at its design ms rate of 240 


| 
| 
| 


DISCUSSION 


(912 x 10° L/day). 1 ‘The dry diffuser was tested at its average 
_ discharge rate of 107 mgd (350 x we L/day) and peak discharge rate of 190 a 
‘The test results of the physical ‘compared favorably with the predictions 
of the mathematical model developed by Koh (7, 8) for the case of no current. 
A This writer further utilized the test results to compare with the predictions — 
} of initial dilution by using the author’s “‘finite line plume” model. 


Table 1 shows the comparison of the physical model test results with the 


PREDICTED FLOW | WEIGHTED AVERAGE DILUTION ody sud 


‘ 
OMPARED WITH ROBERTS * MODEL 


FIG. 15. 15. —Comparison of Test with Math Models’ Predictions oD 


‘mathematical m model’s pre predictions. The test results c clearly indicate that the effect 
_ Of ocean currents is to increase dilution. Fig. 14 shows the selected density — 
- profiles to be simulated in the tests. The comparison between the flow-weighted 
average dilution predicted by Koh’s model and the author’s model and the 

h | minimum dilution measured in the laboratory (minimum in the vertical profile, : 
temporally and laterally averaged) i is shown on Fig. 15. T The comparisons indicate 
both Koh’s model and the author’ s model are conservative, and are useful 
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8. Koh, R. C. Y., and Fan, L. N., ‘‘Mathematical Models for the Prediction of Temperature ibe 
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DOXES IN THE His sToRY OF was 


ohammad A. Gi, 2M. ASCE 


hydraulics written: in his usual masterly style. History described in the form 
f paradoxes is ve very interesting indeed. Attention is drawn here toa few parado xes 
which are not perpetuated yet (hopefully?). 
_ The hydraulicians are well acquainted with the Shields cirtical shear stress, 
st Shields parameter, and Shields diagram, thanks to an increase in research od 
; in sediment transport in the last few decades. However, Shields assuredly was 
me not the first investigator who introduced the parameter of | critical shear stress 


os author has presented another most interesting paper on the history of 


textbooks which describe Shields critical shear stress also ‘describe | in . detail: 
_ the Du Boys bed load formula, after a chapter or so. In the derivation of 
. Du Boys formula, the critical shear stress not moody makes a formal appearance 


a any standard, is the functional correlation of what is now ; called Shields 
parameter with the grain Reynolds number on the basis of experimental data z 

Shields collected for his doctoral thesis. Paradoxically, Shields, a mechanical 

engineer, did not realize the importance of his contribution (unwittingly?). 

The story was told by the author of how Ippen’s name came to be associated 

a with the equation of suspended sediment in one-dimensional, steady, » uniform, — 

“open channel flow, although the work had been developed independently and 
iuhed by the author. On the other hand, Ippen never published his work. s 

- Fortunately, this paradox appears to have been laid to rest. Ee ie 

What of the Einstein-Brown (or the other way round) and Kalinske-Brown | 

sediment transport formulas? The writer had the occasion to use these formulas 
im some of his publications and the formulas were described, as in the foregoing, q 


= “June, 1980, by Hunter Rouse (Proc. Paper 15475), 


_ *Prof. of Hydr. E , Civ. E Dept., Ahmadu Bello Ur _ Zaria, Ni 
of Hydr. iv. Engrg. ello Univ., Zaria, Nigeria. 


> 


| 
| 
| 


"ment note appearing at the end of the chapter on ‘‘Sediment Transportation” we 
+ in Ref. _ 1 came to the writer’s attention . Although the chapter is credited to 
eS C. Brown, the sections in which the previously mentioned formulas were given 
a for the first time were written by no other f¢ person than the author of the original e., 
ts paper. As far as the writer is concerned, the paradox had been born and probably Q 
# 
2 perpetuated. Just before writing this discussion, it was discovered that at least rt 
one book (2), if not more, ceepee the foregoing formulas as Rouse-Einstein 
and Rouse-Kalinske formulas, , respectively . The author has the right and the 
privilege to | provide the background intormetion to at least set the facts i in 


wets 


Brown, C. B., Sediment Engineering Hydraulics, H. 
John Wiley and Sons, Inc., New York, N.Y., 1950, pp. 769-857. a ae = 
2. Morris, H. M., Applied Hydraulics in Engineering, Press Co., New York 
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